Wydział Matematyki i Nauk Informacyjnych 


Politechnika Warszawska 


Algorytmy 
i podstawy programowania 


- wykłady 


Marek Gągolewski 


Y KAPITAL LUDZKI (E| PROGRAM ROZWOJOWY "noes EF] 
NARODOWA STRATEGIA SPÓJNOŚCI (5) POLITECHNIKI WARSZAWSKIEJ FUNDUSZ SPOŁECZNY 


Projekt współfinansowany przez Unię Europejską w ramach Europejskiego Funduszu Społecznego 


| PROG RAM ROZWOJOWY 
U] POLITECHNIKI WARSZAWSKIEJ 


ALGORYTMY 
I PODSTAWY 
PROGRAMOWANIA 


Marek Gągolewski 


I. Etapy tworzenia oprogramowania. Algorytm 

II. Podstawy organizacji i działania komputerów 
III. Deklaracja zmiennych w języku C++. Operatory 
IV. Instrukcja warunkowa i pętle 

V. Tablice jednowymiarowe. Sortowanie 

VI. Funkcje cz. I 

VII. Funkcje cz. II. Rekurencja 

VIII. Wskaźniki. Dynamiczna alokacja pamięci 

IX. Macierze 

X. Podstawowe abstrakcyjne struktury danych 


UNIA EUROPEJSKA 
KAPITAŁ LUDZKI EUROPEJSKI 
CZŁOWIEK — NAJLEPSZA INWESTYCJA! FUNDUSZ SPOŁECZNY 


Projekt współfinansowany przez Unię Europejską w ramach Europejskiego Funduszu Społecznego 


"PROG RAM ROZWOJOWY 
la] POLITECHNIKI WARSZAWSKIEJ 


ALGORYTMY 
I PODSTAWY 
PROGRAMOWANIA 


Marek Gągolewski 


I. Etapy tworzenia oprogramowania. Algorytm 


KAPITAŁ LUDZKI UNIA EUROPEJSKA 
CZŁOWIEK — NAJLEPSZA INWESTYCJA! EUROPEJSKI 
FUNDUSZ SPOŁECZNY 


Projekt współfinansowany przez Unię Europejską w ramach Europejskiego Funduszu Społecznego 


Spis treści 


1 Etapy tworzenia oprogramowania 


1.1 Sformułowanie 1 analiza problemu] ............. o... 11 
1.1.1 Komputery 


1.2 Projektowanie 


pu 


1.2.1 Przykład: Problem młodego Gaussa| ..................-. 
1.3 Implementacja 
ESA +05 6Ro GO GE ay ay PY $E BóĘ WÓZ GR a 


EA «owsa ii RW 4 AYO kW KSW kW KE 


16 Podsumowanie 


O © © © JU W W NM 


2 Cwiczenia 


p 
= 


— 
ww 
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1 Etapy tworzenia oprogramowania 


1.1 Sformułowanie i analiza problemu 


Punktem wyjścia naszych rozważań jest pewne zagadnienie. Problem ten może być bardzo 
ogólnej natury, np. numerycznej, logicznej, społecznej czy nawet egzystencjalnej. Oto kilka 
przykładów: 


— obliczenie wartości jakiegoś wyrażenia arytmetycznego, 

— znalezienie rozwiązania układu 10000 równań różniczkowych, 

— zaplanowanie najkrótszej trasy podróży od miasta X do miasta Y, 

— postawienie diagnozy medycznej na podstawie listy objawów chorobowych, 

— rozpoznanie twarzy przyjaciół na zdjęciu, 

— dokonanie predykcji wartości indeksu giełdowego, 

— zaliczenie przedmiotu Algorytmy i Podstawy Programowania, 

— rozwiązanie dylematu filozoficznego, np. czym jest szczęście i jak być szczęśliwym? 


Znalezienie rozwiązania danego problemu jest często dla nas bardzo istotne. Rzecz jasna, 
kluczowym pytaniem jest: jak to zrobić? 


Sytuację tę obrazuje rys. 
PROBLEM 


297 


ROZWIAZANIE 


Rysunek 1: Punkt wyj$cia naszych rozwazañ. 


Po dosé abstrakcyjnym sformutowaniu problemu przechodzimy do jego analizy. Tutaj do- 
precyzowujemy, o co nam naprawde chodzi, co rozumiemy pod pewnymi pojeciami, jakich 
wyników się spodziewamy i do czego ewentualnie mogą one się nam przydać. 

W przypadku pewnych zagadnień (np. matematycznych) zadanie czasem wydaje się względ- 
nie proste. Wszystkie pojęcia mają swoją definicję formalną, można udowodnić, że pewne kroki 
prowadzą do spodziewanych wyników, które są jednoznaczne itd. 

Jednakże inne zagadnienia mogą przytłaczać swoją złożonością. Na przykład „zaliczenie 
przedmiotu AiPP” wymaga określenia, jaki stan końcowy jest pożądany (ocena bardzo dobra?), 
jakie czynności są kluczowe do osiągnięcia rozwiązania (udział ćwiczeniach i laboratoriach, 
słuchanie wykładu, zadawanie pytań, dyskusje na konsultacjach), jakie czynniki mogą wpływać 
na powodzenie na poszczególnych etapach nauki (czytanie książek, wspólna nauka?), a jakie je 
wręcz uniemożliwiać (codzienne imprezy? brak prądu w akademiku? popsuty komputer?). 
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1.1.1 Komputery 


Istotną cechą niektórych problemów jest to, że nadają się do rozwiązania za pomocą kompu- 
tera. Wbrew pozorom, jak pokazuje rozwój współczesnej informatyki, takich zagadnień jest 
wcale niemało. Często też mamy do czynienia z sytuacją, w której komputery mogą pomóc 
uzyskać rozwiązanie częściowe lub zgrubne przybliżenie rozwiązania, które może być lepsze 
niż zupełny brak rozwiązania. 

Często słyszymy o niesamowitych „osiągnięciach komputerów, np. wygraniu w szachy 
z mistrzem Świata, uczestnictwie w pełnym podchwytliwych pytań teleturnieju Jeopardy! (pier- 
wowzór niegdyś emitowanego w Polsce Va Banque), samodzielnym sterowaniu pojazdem ko- 
smicznym, wirtualną obsługą klienta w banku itp. Wyobraźnię o ich potędze podsycają od kil- 
kudziesięciu lat opowiadania science-fiction, których autorzy przepowiadają, że te maszyny — 
kiedyś — będą potrafiły zrobić prawie wszystko, co jest do pomyślenia. 

Niestety, oprócz tego, komputery podlegają trzem następującym ograniczeniom. By łatwiej 
można było je sobie uzmysłowić, posłużymy się analogią z dziedziny motoryzacji. 


— Komputer ma ograniczoną moc obliczeniową. Każda instrukcja wykonuje się przez pe- 
wien czas. Im bardziej złożone zadanie, tym jego rozwiązywanie trwa dłużej. Wraz z roz- 
wojem techniki, sytuacja ta jednak się poprawia. (Samochody mają np. ograniczoną pręd- 
kość, ograniczone przyspieszenie). 


— Komputer ,,rozumie” określony język (języki), do którego syntaktyki (składni) trzeba się 
dostosować, którego konstrukcję trzeba poznać, by móc się z nim „dogadać ”. Języki sa 
zdefiniowane formalnie za pomocą ściśle określonej gramatyki. Nie toleruje on najczę- 
Ściej żadnych odstępstw lub toleruje tylko nieliczne. (W samochodach, aby zwiększyć 
obroty silnika, należy wykonać jedną ściśle określoną czynność — wcisnąć mocniej pe- 
dał gazu. Nic nie da uśmiechnięcie się bądź próby sympatycznej konwersacji na temat 
zalet jazdy z inną prędkością). 


— Komputer ograniczony jest ponadto przez tzw. czynnik ludzki. potrafi zrobić tylko tyle 
(i aż tyle), ile mu sami dokładnie powiemy, co ma zrobić, krok po kroku, instrukcja po 
instrukcji. Nie domyśli się, co tak naprawdę nam chodzi po głowie. Każdy rozkaz mu wy- 
dawany ma bowiem określoną semantykę (znaczenie). On jedynie potrafi go posłusznie 
wykonać. (Samochód zawiezie nas gdzie chcemy, jeśli będziemy odpowiednio posługi- 
wać się kierownicą i innymi przyrządami. Nie będzie protestował, gdy skręcimy nie na 
tym skrzyżowaniu, co trzeba). 


Celem głównym naszych wszystkich rozważań jest więc takie poznanie natury i języka kom- 
puterów, by mogły zrobić dokładnie to, co my chcemy. 


Na początek, dla porządku, przedyskutujemy definicję obiektu naszych zainteresowań. Otóż 
słowo komputer pochodzi od łacińskiego czasownika computare, który po prostu oznacza obli- 
czać. Jednak powiedzenie, że komputer zajmuje się tylko obliczaniem, to stanowczo za wąskie 
spojrzenie. 


Definicja 1. Komputer to programowalne urządzenie elektroniczne służące do przetwarzania 
informacji. 
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Aż cztery słowa w tej definicji wymagają wyjaśnienia. Programowalność to omówiona po- 
wyżej zdolność do przyjmowania, interpretowania i wykonywania wydawanych poleceń zgod- 
nie z zasadami syntaktyki i semantyki używanego języka. 

Po drugie, pomimo wykonywania licznych prób m.in. przez biologów molekularnych i fi- 
zyków kwantowych, współczesne komputery to w znakomitej większości maszyny elektro- 
niczne. O implikacjach tego faktu dowiemy się więcej z drugiego wykładu poświęconego or- 
ganizacji i działaniu tych urządzeń. 

Dalej, jednostkę informacji rozumiemy jako ciąg liczb lub symboli. Oto kilka przykładów: 
6 (liczba naturalna), PRAWDA (wartość logiczna), "STUDENT MINT” (napis, czyli ciąg zna- 
ków drukowanych), 3,14159 (liczba rzeczywista), (4,-3,-6,34) (ciąg liczb naturalnych), 011100 
(ciąg liczb binarnych), (' WARSZAWA”, 52*13"56”N, 21700'307E) (współrzędne GPS pew- 
nego miejsca na mapie). 

Wartym zanotowania faktem jest to, iż wiele obiektów spotykanych na co dzień może mieć 
swoje reprezentacje w postaci ciągów liczb lub symboli. Często te reprezentacje nie muszą 
odzwierciedlać ich wszystkich cech, lecz tylko te, które są potrzebne w danym zagadnieniu. 
Ciąg ( Jola Kowalska”, "Matematyka II rok”, 4,92, 21142) może być wystarczającą informacją 
dla pani z dziekanatu, by przyznać pewnej studentce stypendium naukowe. W tym przypadku 
kolor oczu Joli jest zbędnym szczegółem. 

I wreszcie, przetwarzanie informacji to wykonywanie różnych działań na liczbach (np. 
operacje arytmetyczne, porównania) lub symbolach (np. zamiana elementów, łączenie ciągów). 
Rodzaje wykonywanych operacji są Ściśle określone, co nie znaczy, że nie można próbować 
samodzielnie tworzyć nowych na podstawie tych, które już są dla nas dostępne. 


2+2 => 4 
T<e > FAŁSZ 
”KATA”O9 © STREFA”(E>0)) — "KATASTROFA" 


są przykładowymi operacjami przetwarzającymi, odpowiednio: liczby naturalne, wartości lo- 
giczne i napisy. 


Ciąg operacji, które potrafi wykonać komputer, nazywamy programem komputerowym. 


Informacje wejściowe 
PROBLEM 


PROGRAM 
ROZWIĄZANIE 
NĄ, Informacje wyjściowe 


Rysunek 2: Rozwiązanie problemu za pomocą programu komputerowego. 
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Przyjrzyjmy się diagramowi na rys. |2| Otóż problem, który da się rozwiązać za pomocą 
komputera posiada dwie istotne cechy: 


a) daje się wyrazić w postaci pewnych danych (informacji) wejściowych, 

b) istnieje dla niego program komputerowy, przetwarzający dane wejściowe w taki sposób, 
że informacje wyjściowe mogą zostać przyjęte jako reprezentacja poszukiwanego roz- 
wiązania. 


1.2 Projektowanie 
Po etapie analizy wstępnej następuje etap projektowania algorytmów. 


Definicja 2. Algorytm to abstrakcyjny przepis (proces, metoda, „instrukcja obsługi”) pozwa- 
lający na uzyskanie, za pomocą skończonej liczby działań, oczekiwanych danych wyjściowych 
na podstawie poprawnych danych wejściowych. 


Przykładowym algorytmem ,,z życia” jest przepis na ciasto marchewkowe] przedstawiony 
w tab. [I] 

Widzimy, jak wygląda zapisany algorytm. Ważną umiejętnością, którą będziemy ćwiczyć, 
jest odpowiednie jego przeczytanie. Dobrze to opisał D. E. Knuth (2002, s. 4): 


Na wstępie trzeba jasno powiedzieć, że algorytmy to nie beletrystyka. Nie należy 
czytać ich ciurkiem. Z algorytmem jest tak, że jak nie zobaczysz, to nie uwierzysz. 
Najlepszą metodą poznania algorytmu jest wypróbowanie go. 


A zatem — do kuchni! 
Oto najistotniejsze cechy algorytmu (por. Knuth, 2002, s. 4): 
— skończoność — wykonanie algorytmu musi zatrzymać się po skończonej liczbie kroków; 


— dobre zdefiniowanie — każdy krok algorytmu musi być opisany precyzyjne, Ściśle i jed- 
noznacznie, tj. być sformułowanym na takim poziomie ogólności, by każdy, kto będzie 
go czytał, był w stanie zrozumieć, jak go wykonać; 


— ściśle określone dane wejściowe pochodzące z odpowiednio określonych zbiorów; 


— dane wyjściowe, czyli wartości powiązane z danymi wejściowymi, powinny odpowiadać 
specyfikacji oczekiwanego poprawnego rozwiązania; 


— efektywność — w algorytmie powinno unikać się operacji niepotrzebnie wydłużających 
czas wykonania, w miarę możliwości należy zastępować je prostszymi bądź w ogóle je 
usuwać. 


!Przepis ten pochodzi z serwisu Kwestia Smaku: 
http://www.kwestiasmaku.com/desery/ciasta/ciasto_marchewkowe/przepis.html. 
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Tablica 1: Przepis na ciasto marchewkowe. 


Dane wejściowe: Ciasto: 2 jaja, 200 g brązowego cukru, 150 ml oleju roślinnego, 200 g 
drobno startej marchewki, 50 g posiekanych orzechów włoskich, 75 g drobno pokrojo- 
nego ananasa (świeżego lub z puszki), 50 g wiórków kokosowych, 200 g mąki, po 1 ły- 
żeczce: cynamonu, sody, soli, 1/2 łyżeczki proszku do pieczenia. Polewa: 400 g cukru 
pudru (można mniej), 100 g kremowego serka, np. Philadelphia, 50 g masła. 


Dane wyjściowe: Ciasto marchewkowe. 


Wykonanie: 
Polewa: 


a) Mikserem ucieram serek wraz z masłem. 

b) Dodaję cukier puder, w 3 częściach, cały czas ucierając pomiędzy kolejnymi par- 
tiami. 

c) Wstawiam do lodówki, aby polewa lekko stężała. 


Ciasto: 


a) Ubijam jajka do podwojenia objętości. Dodaję cukier i dalej ubijam aż masa bę- 
dzie gładka i puszysta. Wciąż ubijając na wysokich obrotach, dolewam ciągłym, 
cieniutkim strumieniem olej. 

b) Do powstałej masy dodaję marchewkę, ananasa, orzechy, wiórki kokosowe i deli- 
katnie mieszam. 

c) Dodaję przesianą mąkę, cynamon, sodę, proszek do pieczenia i sól, delikatnie łączę 
wszystkie składniki. Przekładam do małej formy 21 x 21 cm, wyłożonej papierem 
do pieczenia. 

d) Piekę przez 1 godzinę, w piekarniku nagrzanym do 150 stopni. 

e) Wystudzone ciasto przekrawam poziomo na 2 części. Spód smaruję 1/3 ilości po- 
lewy. Przykrywam górą ciasta i smaruję resztą polewy. Wstawiam do lodówki. 
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Wygląda na to, że nasz przepis na ciasto marchewkowe posiada wszystkie te cechy. Wykonując 
powyższe czynności ciasto to uda nam się kiedyś zjeść (a więc jest skończony). Wszystkie 
czynności są zrozumiałe nawet dla mało wprawionej gospodyni. Dane wejściowe i wyjściowe 
są dobrze wyspecyfikowane. Sam proces przygotowania jest efektywny (nie każe np. kucharce 
umalować się w trakcie mieszania składników bądź pojechać na zagraniczną wycieczkę). 


Jeden program rozwiązujący rozpatrywany problem może zawierać realizację wielu algo- 
rytmów, np. gdy złożoność zagadnienia wymaga podzielenia go na kilka podproblemów. I tak 
zdolna kucharka wykonująca program „obiad” powinna podzielić swą pracę na podprogramy 
„rosół”, „kotlet schabowy z frytkami” oraz „ciasto marchewkowe”, 

Ważne jest, że algorytm nie musi być wyrażony za pomocą języka zrozumiałego przez kom- 
puter. Taki sposób opisu algorytmów nazywa się często pseudokodem. Stanowi on etap po- 
średni między analizą problemu a implementacją, opisaną w kolejnym paragrafie. Ma on nam 
po prostu pomóc w bardziej formalnym podejściu do tworzenia programu. 

Dla przykładu, w przytoczonym powyżej przepisie, nie jest dokładnie wytłumaczone — 
krok po kroku — co oznacza ,,pieke przez jedną godzinę”. Wszak wymaga ono wielu czynno- 
$ci (włączenie piekarnika, rozgrzanie, pilnowanie ustawionego zegarka itp.). Jednakże, jak już 
wspomnieliśmy, czynność ta jest dobrze zdefiniowana. 


Jakby tego było mało, może istnieć wiele algorytmów służących do rozwiązania tego sa- 
mego problemu! Przyjrzyjmy się następującemu przykładowi. 


1.2.1 Przykład: Problem młodego Gaussa 


Szeroko znana jest historia Karola Gaussa, z którym miał problemy jego nauczyciel mate- 
matyki. Aby zająć czymś młodego chłopca, profesor kazał mu wyznaczyć sumę liczb a, a + 
1,...,b, gdzie a,b € N (oryginalnie było to a = 1 i b = 100). Zapewne myślał on, że tamten 
użyje algorytmu I i spędzi nad tym trochę czasu. 


Algorytm I wyznaczania sumy kolejnych liczb naturalnych. 


// Wejście: a,bEŃ (a<b) 
// Wyjście: a+a+1+..+b0EN 


niech suma,iEN; 


suma = 0; 
dla (i = a,a +1 ,... ,b) 
( 

suma = suma + i; 


) 


zwróć suma jako wynik; 


Jednakże sprytny Gauss zauważył, żea+a+1+:::+b= 2 (b — a + 1). Dzięki temu 
znalazł on bardzo szybko rozwiązanie korzystając z algorytmu II. 
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Algorytm II wyznaczania sumy kolejnych liczb naturalnych. 


// Wejście: a,bEŃ (a <b) 
// Wyjście: a+a+1+..+bEN 


niech suma € N; 
suma = B(b=a+1); 


zwróć suma jako wynik; 


Zauważmy, że obydwa algorytmy rozwiązują poprawnie to samo zagadnienie. Jednakże 
inna jest liczba operacji arytmetycznych (+, —, *, /) potrzebna do uzyskania oczekiwanego wy- 
niku. Poniższa tabelka zestawia tę miarę efektywności obydwu rozwiązań dla różnych danych 
wejściowych. Zauważmy, że w przypadku pierwszego algorytmu liczba wykonywanych ope- 
racji dodawania jest równa (b — a + 1) [instrukcja suma = suma + i] + (b — a) [dodawanie 
występuje także w pętli dla... ]. 


Tablica 2: Liczba operacji arytmetycznych potrzebna do znalezienia rozwiązania problemu mło- 
dego Gaussa. 


a b  Alg.l Alg.Il 
1 10 19 5 
1 
1 


100 199 5 
1000 1999 5 


Charakteryzacją efektywności algorytmów zajmuje się dziedzina zwana analizą algoryt- 
mów. Czerpie ona szeroko z wyników takich działów jak matematyka dyskretna czy rachunek 
prawdopodobieństwa. Z jej elementami zapoznamy się podczas innych wykładów. 


1.3 Implementacja 


Na kolejnym etapie, abstrakcyjne algorytmy zapisane często w postaci pseudokodu (np. za po- 
mocą języka polskiego) należy zapisać w formie, która jest zrozumiała nie tylko przez nas, ale 
i przez komputer. Nadto, należy dokonać scalenia (powiązania) podrozwiązań w jeden spójny 
projekt. 

Intuicyjnie, tutaj wreszcie tłumaczymy komputerowi, co dokładnie chcemy, by zrobił, tzn. 
go programujemy. Efektem naszej pracy będzie kod źródłowy programu. 

Formalnie rzecz ujmując, zbiór zasad określających, jaki ciąg symboli tworzy kod źró- 
dłowy programu, nazywamy językiem programowania. Reguły składniowe (ang. syntax) Ści- 
śle określają, które (i tylko które) wyrażenia są poprawne. Reguły znaczeniowe (ang. seman- 
tics) określają precyzyjnie, jak komputer ma rozumieć dane wyrażenia. Dziedziną zajmującą 
się analizą języków programowania jest lingwistyka matematyczna. 


Podczas tego wykładu będziemy poznawać język C++, zaprojektowany ok. 1983 r. przez 
B. Stroustrupa (aktualny standard: ISO/IEC 14882:2003). Język ten powstał jako rozwinięcie 
języka C i jest przedstawicielem klasy języków cechujących się podobną doń składnią. Wśród 
innych podobnych można wymienić inne popularne narzędzia, takie jak Java, CH, PHP czy 
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JavaScript. Widzimy zatem, że C++ jest tylko jedną z wielu możliwości wydawania poleceń 
komputerowi, jednakże spośród innych wyróżnia się on zwięzłością, efektywnością i względną 
łatwością nauki. 

Kod źródłowy zapisujemy zwykle w postaci kilku plików tekstowych (tzw. plików Zródto- 
wych, ang. source files). Pliki te można edytować za pomocą dowolnego edytora tekstowego, 
np. Notatnik systemu Windows. Mimo to wygodniejsze jest korzystanie z całych środowisk 
programistycznych. My podczas laboratoriów będziemy używać Microsoft Visual C++] 


I wreszcie narzędziem, które pozwala przetworzyć pliki źródłowe (zrozumiałe przez czło- 
wieka i komputer) na kod maszynowy programu komputerowego (zrozumiały tylko przez kom- 
puter) nazywamy kompilatorem (ang. compiler). 


1.4 Testowanie 


Gotowy program należy przetestować, to znaczy sprawdzić, czy dokładnie robi to, o co nam 
chodziło. Jest to, niestety, często najbardziej żmudny etap tworzenia oprogramowania. Jeśli coś 
nie działa, jak powinno, należy wrócić do któregoś z poprzednich etapów i naprawić błędy. 

Warto zwrócić uwagę, że przyczyn niepoprawnego działania należy zawsze szukać w swo- 
jej pracy, a nie liczyć na to, że jakiś chochlik robi nam na złość. Jak powiedzieliśmy, komputer 
robi tylko to, co my mu każemy. Zatem jeśli mu nakazaliśmy wykonać instrukcję, której skut- 
ków ubocznych nie jesteśmy do końca pewni, jest to nasza odpowiedzialność, żeby te skutki 
opanować. 


1.5 Eksploatacja 


Gdy program jest przetestowany, może służyć do rozwiązywania wyjściowego problemu (wyj- 
ściowych problemów). Czasem zdarza się, że na tym etapie dochodzimy do wniosku, iż czegoś 
nam brakuje lub że nie jest to do końca, o co nam na początku chodziło. Wtedy oczywiście 
pozostaje znów powrót do wcześniejszych etapów pracy. 

Jak widać, nauka programowania komputerów może nam pomóc rozwiązać wiele proble- 
mów, których rozwiązanie bez nich byłoby często niedostępne. Poza tym jest wspaniałą roz- 


rywką! 


1.6 Podsumowanie 


Omówiliśmy następujące etapy tworzenia oprogramowania, które są niezbędne do rozwiązania 
zagadnień za pomocą komputera: 


— sformułowanie i analiza problemu, 
— projektowanie, 

implementacja, 

— testowanie, 

eksploatacja. 


| 


| 


2Do pobrania bezpłatnie ze strony http://www.microsoft.com/express/Downloads/. 
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Efekty pracy po każdym z etapów podsumowuje rys. 


Godne polecenia rozwinięcie omawianego tematu można znaleźć w książce Harela (2001, 
s. 9-31). 


PROBLEM 


Analiza IDEE 


| RSE 


Rysunek 3: Efekty pracy po każdym z etapów tworzenia oprogramowania. 


Projektowanie 


Implementacja 


Testowanie 


Eksploatacja 
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2 Cwiczenia 


Zadanie 1.1. Dany jest algorytm Euklidesa znajdywania największego wspólnego dzielnika 
(NWD) dwóch liczb a,b € Z, 0 < a < b. 


// Wejście: O<a<b 
// Wyjście: NWD(a,b) 


niech cEN; 

dopóki (a #0) 

{ 
c = reszta z dzielenia b przez a; 
b =a; 
G=ć; 


zwróć b jako wynik; 


Prześledź działanie algorytmu Euklidesa (jakie wartości przyjmują zmienne a,b, c w każ- 
dym kroku) znajdywania największego wspólnego dzielnika dla następujących par liczb 


a) 42, 56, c) 199, 544, 
b) 192, 348, d) 2166, 6099. 


Zadanie 1.2. Pokaż, w jaki sposób za pomocą ciągu przypisań można przestawić wartości 
dwóch zmiennych (a, b), by otrzymać (b, a). 


Zadanie 1.3. Pokaż, w jaki sposób za pomocą ciągu przypisań można przestawić wartości 
trzech zmiennych (a, b, c), by otrzymać (c, a,b). 


Zadanie 1.4. Pokaż, w jaki sposób za pomocą ciągu przypisań można przestawić wartości 
czterech zmiennych (a, b, c, d), by otrzymać (c, d, b, a). 


Zadanie 1.5. Dany jest ciąg n liczb rzeczywistych x = (x[0], x[1],...,u[n— 1]) (umawiamy 
się, że elementy ciągów numerujemy od 0). Rozważmy ich średnią arytmetyczną, określoną 
jako 

1 
-X zi] = „ lO] + z[1] +: +ajn— 1)). 


Rozważmy następujący algorytm służący do jej wyznaczania. 


// Wejście: n>0 oraz zf0j,z|1),...,zjn-1]ER 
// Wyjście: średnia arytmetyczna elementów a|0),zf1],...,z|n — 1]) 
niech sredniaarytm€ R; 
niech ¡EN; 


sredniaarytm = 0; 
dla (i=0,1,...,n—1) 

sredniaarytm = sredniaarytm + x[i]; 
sredniaarytm = sredniaarytm / n; 


zwróć sredniaarytm jako wynik; 
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Wyznacz za pomocą powyższego algorytmu wartość średniej arytmetycznej dla ciągów (1, —1, 
2, 0, —2) oraz (34, 2, —3, 4, 3,5). 


Zadanie 1.6. Dany jest ciąg n dodatnich liczb rzeczywistych x = (x/0], x[1],...,z[n — 1)) 
(różnych od 0). Napisz algorytm, który wyznaczy ich średnią harmoniczną, określoną jako 


EG 1/40) + 1/20] +---+1/2[ín- 1] 


Wyznacz za pomocą tego algorytmu wartość średniej harmonicznej dla ciągów (1,4, 2,3, 1) 
oraz (10, 2, 3, 4). 


x Zadanie 1.7. Dany jest ciąg n liczb rzeczywistych x = (z(0|, x[1],...,u[n—1]). Rozważmy 
tzw. sumę kwadratów odchyleń elementów od ich średniej arytmetycznej, określoną jako 


SKO(x) = > (s — (ESen) 


Rozważmy następujący algorytm służący do wyznaczania SKO. 


// Wejście: n>0 oraz zf0j,z|1),...,zjn-1]EeR 
|// Wyjście: SKO(x[0], x[1],...,x[n— 1]) 

3| niech sko, sredniaarytm ER; 

niech i,j EN; 


sko = 0; 
dla (i=0,1,...,n—1) 
( 
sredniaarytm = 0; 
dla (j=0,1,...,n—l1) 
( 
sredniaarytm = sredniaarytm + x[/]l; 


) 


sredniaarytm = sredniaarytm Í n, 


sko = sko + (x[i] — sredniaarytm)*(x[i] — sredniaarytm); 


) 


zwróć sko jako wynik; 


a) Wyznacz za pomocą powyższego algorytmu wartość SKO dla x = (5,3, —1, 7, —2). 

b) Policz, ile łącznie operacji arytmetycznych (+-, —, *, /) potrzebnych jest do wyznaczenia 
SKO dla ciągów wejściowych o n = 5, 10, 100, 1000, 10000 elementach. Wyraź tę liczbę 
jako funkcję długości ciągu wejściowego n. 

c) Zastanów się, jak usprawnić powyższy algorytm, by nie wykonywać wielokrotnie zbęd- 
nych obliczeń. Ile teraz będzie potrzebnych operacji arytmetycznych? 
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N 


3 Wskazówki i odpowiedzi do ćwiczeń 


Odpowiedź do zadania 1.1] 
NWD(42,56) = 14. 


5: a=42, b=56, c=0 
7: a=42, b=56, c=14 
8: a=42, b=42, c=14 
9: a=14, b=42, c=14 
5: a=14, b=42, c=14 
7: a=14, b=42, c=0 
8: a=14, b=14, c=0 
9: a=0, b=14, c=0 

53 14, c=0 


NWD(192,348) = 12. 


5: a=192, b=348, c=0 
7: a=192, b=348, c=156 
8: a=192, b=192, c=156 
9: a=156, b=192, c=156 
5: a=156, b=192, c=156 
7: a=156, b=192, c=36 
8: a=156, b=156, c=36 
9: a=36, b=156, c=36 
5: a=36, b=156, c=36 
7: a=36, b=156, c=12 
8: a=36, b=36, c=12 

9: a=12, b=36, c=12 

5: a=12, b=36, c=12 

7: a=12, b=36, c=0 

8: a=12, b=12, c=0 

9: a=0, b=12, c=0 

5: a=0, b=12, c=0 
Wynik: 12 


NWD(199, 544) = 1. 
NWD(2166, 6099) = 57. 


Odpowiedź do zadania|1.2 


// Wejście: (a,b) 

// Wyjście: (a ',b )=(b,a) 
niech x — zmienna pomocnicza; 
XS g az $ => 
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N 


w 


N 


w 


Odpowiedź do zadania |1.3 


// Wejście: (a,b,c) 
// Wyjście: (a ,b',c )=(c,a,b) 


niech x — zmienna pomocnicza; 
+ b; 
„RZE: 
a= g: 
CSa 


Odpowiedź do zadania 1.4] 


// Wejście: (a,b,c,d) 


// Wyjście: (a',b',c ',d')=(c,d,b,a) 


niech x — zmienna pomocnicza; 


aSa R > 
Il 
= QA 0 


Odpowiedź do zadania 1.5] 
Wynik dla (1, —1, 2,0, —2): 0. 
Wynik dla (34, 2, —3, 4, 3,5): 8,1. 


Odpowiedź do zadania 1.6] 
H . 60 
Wynik dla (1, 4, 2,3, 1): 37. 
Wynik dla (10,2, 3,4): #2. 


Odpowiedź do zadania 1.7] 
SKO(5, 3, —1, 7, —2) = 59,2. 


Podany algorytm wymaga n? + 5n operacji arytmetycznych. Nie jest on efektywny, gdyż 


można go łatwo usprawnić tak, by potrzebnych było 5n + 1 działań (+-, —, 
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1 Historia i organizacja współczesnych komputerów 


1.1 Zarys historii informatyki 


Historia informatyki nieodłącznie związana jest z historią matematyki, a w szczególności za- 
gadnieniem zapisu liczb i rachunkami. Najważniejszym powodem powstania komputerów było, 
jak łatwo się domyślić, wspomaganie wykonywania żmudnych obliczeń. 

Oto wybór istotniejszych wydarzeń, które pomogą nam rzucić światło na przedmiot naszych 
zainteresowań. 


— Ok. 2400 r. p.n.e. używany jest już w Babilonii kamienny pierwowzór liczydła, na którym 
rysowano piaskiem (specjaliści umiejący korzystać z tego urządzenia nie byli liczni). 
Podobne przyrządy w Grecji i Rzymie pojawiły się dopiero ok. V-IV w. p.n.e. Tzw. abaki, 
żłobione w drewnie, pozwalały dokonywać obliczeń w systemie dziesiętnym. 


— Ok. V w. p.n.e. indyjski uczony Pánini sformalizował teoretyczne reguły gramatyki san- 
skrytu. Można ten fakt uważać za pierwsze badanie teoretyczne w dziedzinie lingwistyki. 


— Pierwszy znany nam algorytm przypisywany jest Euklidesowi (ok. 400-300 r. p.n.e.). 
Opisał on operacje, których wykonanie krok po kroku pozwala wyznaczyć największy 
wspólny dzielnik dwóch liczb. 


— Matematyk arabski Al Kwarizmi (IX w.) określa reguły podstawowych operacji arytme- 
tycznych dla liczb dziesiętnych. Od jego nazwiska powstaje potem pojęcie algorytmu. 


— W 1614 r. John Napier (szkocki teolog i matematyk) znalazł zastosowanie logarytmów 
do wykonywania szybkich operacji mnożenia (zastąpił je dodawaniem). W roku 1622 
William Oughtred stworzył suwak logarytmiczny, który jeszcze bardziej ułatwił wykony- 
wanie obliczeń. 


— W. Schickard (1623 r. ) oraz B. Pascal (1645 r.) tworzą pierwsze mechaniczne sumatory. 


— Jacques Jacquard (ok. 1801 r.) skonstruował krosno tkackie sterowane dziurkowanymi 
kartami. Było to pierwsze sterowane programowo urządzenie w dziejach techniki. Po- 
dobny ideowo sposób działania miały pianole (poł. XIX w.), czyli automatycznie stero- 
wane pianina. 


— Do czasów Charlesa Babbage'a jednak żadne urządzenie służące do wykonywania obli- 
czeń nie było programowalne. Ok. 1837-1839 r. opisał on więc wymyśloną przez siebie 
„maszynę analityczną” (parową!), która to umożliwiała. Była to jednak tylko koncepcja 
teoretyczna, nigdy nie udało się jej zbudować. Pomagająca mu Ada Lovelace może być 
uznana za pierwszą programistkę. 


— H. Hollerith (1890) — maszyna wspomagająca spis powszechny w USA. 


— Niemiecki inżynier w 1918 r. patentuje maszynę szyfrującą Enigma. (Nota bene polski 
matematyk Marian Rejewski w 1934 r. łamie jej kod jako pierwszy). 
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— W 1936r. Alan Turing i Alonzo Church definiują formalnie algorytm jako ciąg instrukcji 
matematycznych. Określają dzięki temu to, co daje się policzyć. Kleene stawia później 
tzw. hipotezę Churcha-Turinga, która mówi o tym, że jeśli dla jakiegoś problemu istnieje 
efektywny algorytm korzystający z nieograniczonych zasobów, to da się go wykonać 
na stworzonym przez nich prostym modelu komputera: Każdy problem, który może być 
intuicyjnie uznany za obliczalny, jest rozwiązywalny przez maszynę Turinga (tutaj jednak 
pojawia się problem, co może być uznane za „intuicyjnie obliczalne '). 


— Claude E. Shannon w swojej pracy magisterskiej w 1937 r. opisuje możliwość użycia 
elektronicznych przełączników do wykonywania operacji logicznych (algebra Boole'a). 
Jego praca teoretyczna staje się podstawą konstrukcji wszystkich współczesnych kompu- 
terów elektronicznych. 


— W latach 40” XX w. w Wielkiej Brytanii, USA i Niemczech powstają pierwsze kompu- 
tery elektroniczne, np. Z1, Z3, Colossus, Mark I, ENIAC (zob. rys. [D, EDVAC, EDSAC. 
Pobieraja bardzo duże ilości energii elektrycznej i potrzebują dużych przestrzeni. Np. 
ENIAC miał masę ponad 27 ton, zawierał około 18 000 lamp elektronowych i zajmował 
powierzchnię ok. 140 m”. Były bardzo wolne, np. Z3 wykonywał jedno mnożenie w 3 
sekundy. Pierwsze zastosowanie: łamanie szyfrów, obliczanie trajektorii lotów balistycz- 
nych. 


— W 1947 r. Grace Hopper odkrywa pierwszego robaka (ang. bug) komputerowego w kom- 
puterze Harvard Mark II — dosłownie! 


— W międzyczasie okazuje się, że komputery nie muszą być wykorzystywane tylko w ce- 
lach przyspieszania obliczeń. Powstają teoretyczne podstawy informatyki (np. Turing, 
von Neumann). Zapoczątkowywane są nowe kierunki badawcze, m.in. sztuczna inteli- 
gencja (por. np. słynny test Turinga sprawdzający ,,inteligencje” maszyny). 


— Pierwszy kompilator języka Fortran zostaje uruchomiony w 1957 r. Kompilator języka 
LISP powstaje rok później. 


— Informatyka staje się dyscypliną akademicką dopiero w latach sześćdziesiątych XX w. 
(wytyczne programowe ACM). Pierwszy wydział informatyki powstał na Uniwersytecie 
Purdue w 1962 r. Pierwszy doktorat w tej dziedzinie został obroniony już w 1965 r. 


— W 1969 r. następuje pierwsze połączenie sieciowe pomiędzy komputerami w ramach pro- 
jektu ARPAnet (prekursora Internetu). Początkowo ma ono mieć głównie zastosowanie 
wojskowe. 


— Dalej rozwój następuje bardzo szybko: powstają bardzo ważne algorytmy (np. wyszu- 
kiwanie najkrótszych Ścieżek w grafie Dijkstry, sortowanie szybkie Hoara itp.), teoria 
relacyjnych baz danych, okienkowe wielozadaniowe systemy operacyjne, nowe języki 
programowania, systemy rozproszone i wiele, wiele innych... 


— Współcześnie komputer osobisty (o mocy kilka miliardów razy większej niż ENIAC) 
podłączony do sieci Internet znajdziemy w prawie każdym domu, a elementy informatyki 
wykładane są już w szkołach podstawowych! 
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Rysunek 1: Komputer ENIAC, źródło: http://www.ifj.edu.pl/str/psk/img/eniac4.jpg. 


Wielu znaczących informatyków XX wieku było z wykształcenia matematykami. Informa- 
tyka teoretyczna była początkowo poddziałem matematyki. Z czasem dopiero stała sig samo- 
stanowiącą dziedziną akademicką. Nie oznacza to oczywiście, że obszary badań obu dyscyplin 
są rozłączne. 

Można pokusić się o następującą klasyfikację głównych kierunków badań w informatycd!ł 


a) Matematyczne podstawy informatyki 


i. Kryptografia (kodowanie informacji niejawnej) 
ii. Teoria grafów (np. algorytmy znajdywania najkrótszych Ścieżek, algorytmy koloro- 
wania grafów) 
iii. Logika (w tym logika wielowartościowa, logika rozmyta) 
iv. Teoria typów (analiza formalna typów danych, m.in. badane są efekty związane 
z bezpieczeństwem programów) 


b) Teoria obliczeń 


i. Teoria automatów (analiza abstrakcyjnych maszyn i problemów, które można z ich 
pomocą rozwiązać) 
ii. Teoria obliczalności (co jest obliczalne?) 
iii. Teoria złożoności obliczeniowej (które problemy są „łatwo” rozwiązywalne?) 


c) Algorytmy i struktury danych 


i. Analiza algorytmów (ze względu na czas działania i potrzebne zasoby) 
ii. Projektowanie algorytmów 
iii. Struktury danych (organizacja informacji) 
iv. Algorytmy genetyczne (znajdywanie rozwiązań przybliżonych problemów optyma- 
lizacyjnych) 


!Por. http://www.newworldencyclopedia.org/entry/Computer_science 
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d) Języki programowania i kompilatory 


i. Kompilatory (budowa i optymalizacja programów przetwarzających kod w danym 
języku na kod maszynowy) 
ii. Języki programowania (formalne paradygmaty języków, własności języków) 


e) Bazy danych 


i. Teoria baz danych (m.in. systemy relacyjne, obiektowe) 
ii. Data mining (wydobywanie wiedzy z baz danych) 


f) Systemy równoległe i rozproszone 


i. Systemy równoległe (bada wykonywanie jednoczesnych obliczeń przez wiele pro- 
cesorów) 
ii. Systemy rozproszone (rozwiązywanie tego samego problemu z użyciem wielu kom- 
puterów połączonych w sieć) 
iii. Sieci komputerowe (algorytmy i protokoły zapewniające niezawodny przesył da- 
nych pomiędzy komputerami) 


g) Architektura komputerów 


i. Architektura komputerów (projektowanie, organizacja komputerów, także z punktu 
widzenia elektroniki) 

ii. Systemy operacyjne (systemy zarządzające zasobami komputera i pozwalające uru- 
chamiać inne programy) 


h) Inżynieria oprogramowania 


i. Metody formalne (np. automatyczne testowanie programów) 
ii. Inżynieria (zasady tworzenia dobrych programów, zasady organizacji procesu ana- 
lizy, projektowania, implementacji i testowania programów) 


i) Sztuczna inteligencja 


i. Sztuczna inteligencja (systemy, które wydają się cechować inteligencją) 
ii. Automatyczne wnioskowanie (imitacja zdolności samodzielnego rozumowania) 
iii. Robotyka (projektowanie i konstrukcja robotów i algorytmów nimi sterujących) 
iv. Uczenie się maszynowe (tworzenie zestawów reguł w oparciu o informacje wej- 
ściowe) 
j) Grafika komputerowa 
i. Podstawy grafiki komputerowej (algorytmy generowania i filtrowania obrazów) 
ii. Przetwarzanie obrazów (wydobywanie informacji z obrazów) 


iii. Interakcja człowiek-komputer (zagadnienia tzw. interfejsów służących do komuni- 
kacji z komputerem) 


Do tej listy można też dopisać wiele dziedzin z tzw. pogranicza informatyki, np. zbiory 
rozmyte, bioinformatykę, statystykę obliczeniową itd. 
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Na koniec tego paragrafu warto zwrócić uwagę (choć zdania na ten temat są podzielone), że 
to, co określamy mianem informatyki czy też nauk informacyjnych, jest przez wielu uzna- 
wane za pojęcie szersze od angielskiego computer science, czyli nauk o komputerach. Informa- 
tyka jest więc ogólnym przedmiotem dociekań dotyczących przetwarzania informacji, w tym 
również przy użyciu urządzeń elektronicznych. „Informatyka jest nauką o komputerach w takim 
stopniu jak astronomia nauką o teleskopach” (Dijkstra). 


1.2 Organizacja komputerów 


Jak zdążyliśmy nadmienić, podstawy teoretyczne działania współczesnych elektronicznych kom- 
puterów zawarte zostały w pracy C. Shannona (1937 r.). Mogą być one pojmowane jako zestaw 
odpowiednio sterowanych przełączników, tzw. bitów (ang. binary digits). Stany tych przełącz- 
ników mogą być dwa: 


— prąd nie płynie (co oznaczamy jako 0), 

— prąd płynie (co oznaczamy przez 1). 

Na współczesny komputer osobisty (PC) składają się: 

a) jednostka obliczeniowa (CPU, od ang. central processing unit): 


— jednostki arytmetyczno-logiczne (ALU, ang. arithmetic and logical units), 
— jednostka do obliczeń zmiennopozycyjnych (FPU, ang. floating point unit), 
— układy sterujące (CU, ang. control units), 
— rejestry (akumulatory), pełniące funkcję pamięci podręcznej; 

b) pamięć RAM (od ang. random access memory) — zawiera dane i program, 

c) urządzenia wejściowe (ang. input devices), 

d) urządzenia wyjściowe (ang. output devices). 


Jest to tzw. architektura von Neumanna (por. rys. B). Najważniejszą jej cechą jest to, że 
w pamięci operacyjnej znajdują się zarówno instrukcje programów, jak i dane; to od kontekstu 
zależy, jak powinny być one interpretowane. 


Rysunek 2: Architektura współczesnych komputerów osobistych. 
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2 Reprezentacje liczb całkowitych 


Rozważmy najpierw reprezentacje nieujemnych liczb całkowitych. Systemy liczbowe można 
podzielić na pozycyjne i addytywne. Oto przykłady każdego z nich. 


Pozycyjne systemy liczbowe: 


a) system dziesiętny (decymalny, indyjsko-arabski) — o podstawie 10, 
b) system dwójkowy (binarny) — o podstawie 2, 

c) system szesnastkowy (heksadecymalny) — o podstawie 16, 

d) ws 


Addytywne systemy liczbowe: 


a) system rzymski, 
b) system sześćdziesiątkowy (Mezopotamia). 


Nas szczególnie będą interesować systemy pozycyjne, w tym używany na co dzień system 
dziesiętny. 


2.1 System dziesiętny 


System dziesiętny (decymalny, ang. decimal), przyjęty w Europie od XVI w., to pozycyjny 
system liczbowy o podstawie 10. Do zapisu liczb używa się w nim następujących symboli 
(cyfr): 0, 1, 2, 3, 4,5,6, 7,8,9. 

Jako przykład rozważmy liczbę 194510 (przyrostek ¡y wskazuje, że chodzi nam o liczbę 
w systemie dziesiętnym). Jak w każdym systemie pozycyjnym istotne jest tu, na której pozycji 
stoi każda z cyfr. 


1945, 
= 1 0 0 0 
+ 9 0 0 
+ 4 0 
+ 5 
1945., 
= 1 x 1000 
+ 9 x100 
+ 4 x10 
+ 5 xl 
1945, 
= 1 xi 
+ 9 x10? 
+ 4 x10! 
+ 5 x10 
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Numerację pozycji cyfr zaczynamy od 0. Dalej oznaczamy je w kierunku od prawej do 
lewej, czyli od najmniej do najbardziej znaczącej. 


1945, 


= x10? 
9 x10? 
4 x10! 
| x107 
3 2 1 


Ogólnie rzecz biorąc, mając daną liczbę n cyfrową, zapisaną jako ciąg cyfr dziesiętnych 
bn—1bn—2 . . . b1bo, gdzie b; € {0,1,...,9}dla¿ =0,1,...,n— 1, jej wartość można określić za 
pomocą wzoru 


n—1 
Y b, x 10. 
1=0 


Liczbę 10 występującą we wzorze nazywamy podstawą systemu. 

Jak łatwo się domyślić, można wprowadzić systemy pozycyjne o innych podstawach. My 
szczególnie skupimy się na systemach o podstawie 2 i 16, które są przydatne z punktu widzenia 
informatyki. 


2.2 System dwójkowy 


Dwójkowy pozycyjny system liczbowy (binarny, ang. binary) to system pozycyjny o podsta- 
wie 2. Używanymi weń symbolami (cyframi) są tylko 0, 1. Zauważmy, że cyfry te mają ,,prze- 
łącznikową” interpretację: prąd nie płynie/płynie, przełącznik wyłączony/włączony. Zatem n 
cyfrowa liczba w tym systemie może opisywać stan n przełączników naraz. 

Jako przykład rozpatrzmy liczbę 1010015. 


101001, 


x25 


+ 0 x24 
+ 1 x23 
+ 0 x2? 
+ 0 x21 
+ 1 x2 


101001, 


x32 
0 x16 
1 x8 

x4 

0 x2 
1 xi 


+++ ++ 
© 
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Zatem 1010012 = 3210 + 810 + 110 = 4140. Jest to przykład konwersji (zamiany) podstawy 


liczby dwójkowej na dziesiętną. 


Przy dokonywaniu operacji na liczbach binarnych przydatna może się okazać tablica 


przedstawiająca wybrane potęgi liczby 2. 


Tablica 1: Wybrane potęgi liczby 2. 


k ze 
0 1 
1 2 
2 4 
3 8 
4 16 
5 32 
6 64 
7 128 
8 256 
9 512 
10 1024 
11 2 048 
12 4096 
13 8 192 
14 16 384 
15 32 768 
16 65 536 
31 2 147 483 648 
32 4 294 967 296 
63 | 9223372036 854 775 808 
64 | 18446 744073 709551 616 


Konwersja liczb do systemu dziesiętnego jest stosunkowo prosta. Przyjrzyjmy się jak prze- 
kształcić liczbę dziesiętną na dwójkową. Listing [I]podaje algorytm, który moze byé w tym celu 


wykorzystany. 


Dla przykładu, rozważmy liczbę 194549. Tablica [2| opisuje kolejno wszystkie kroki po- 
trzebne do uzyskania wyniku w postaci dwójkowej. Możemy z niej odczytać, iż 194519 = 


11110011001. 
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Listing 1: Konwersja liczby dziesiętnej na dwójkową. 


// Wejscie: liczba kEeN. 


s| { 


1// Wyjście: postać tej liczby w systemie dwójkowym (kolejne 
cyfry będą ,,wypisywane °? w kolejności od lewej do prawej). 
3 niech i — największa liczba naturalna taka, że k>2'; 
dopóki (i>0) 
jeśli (k>2) 
{ 
wypisz (1); 
k=k-2% 
) 


w_przeciwnym_przypadku 
wypisz(0); 


i=i-1; 
s|} 


// (wynik został ,,wypisany””) 


Tablica 2: Przykład — konwersja liczby 194510 na dwójkowa. 


k % i 
1945 | 2048 

1945 | 1024 10| 1 
1945-1024= 921 | 512 9|1 
921-512= 409 | 256 8| 1 
409-256= 153| 128 7 ]1 
153-128 = 25 64 6|0 
25 32 50 
25 16 4/1 
25-16= 9 8 311 
9-8= 1 4 2|0 
1 2 10 
1 1 OJ1 
1-1= 0 1 


2.3 System szesnastkowy 


System szesnastkowy (heksadecymalny, ang. hexadecimal) jest systemem pozycyjnym o pod- 
stawie 16. Uzywanymi symbolami (cyframi) są: 0,1,...,9,A,B,C,D,E,F. Jak widzimy, wykorzy- 
stywać będziemy także kolejne litery alfabetu. 

Jak można zauważyć, liczba w postaci dwójkowej prezentuje się na kartce lub na ekranie 
dosyć. .. okazale. Jako że jest to czasem problematyczne, w zastosowaniach informatycznych 
zwykło się zapisywać liczby właśnie jako szesnastkowe, gdyż 16 = 2%. Dzięki temu można 
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Tablica 3: Cyfry systemu szesnastkowego i ich wartości w systemie dwójkowym oraz dziesięt- 
nym. 


BIN DEC HEX BIN DEC HEX 
0000 0 0 1000 8 8 
0001 1 1 1001 9 9 
0010 2 2 1010 10 A 
0011 3 3 1011 11 B 
0100 4 4 1100 12 C 
0101 5 5 1101 13 D 
0110 6 6 1110 14 E 
0111 7 7 1111 15 F 


użyć jednego symbolu zamiast 4 w systemie dwójkowym. Fakt ten sprawia również, że kon- 
wersja z systemu binarnego na heksadecymalny i odwrotnie jest bardzo prosta. 
Tablica[B]przedstawia cyfry systemu szesnastkowego i ich wartości w systemie dwójkowym 
oraz dziesiętnym. 
Zatem, np. 10111102 = 5E16, gdyż grupując kolejne cyfry dwójkowe czwórkami, od prawej 
do lewej, otrzymujemy 
RCA 
01011110. 
Podobnie postępujemy dokonując konwersji w przeciwną stronę, np. 2D;6 = 1011013, bowiem 


0010 1101 
RZ 
2 D y. 


Jak widzimy, dla uproszczenia zapisu początkowe zera możemy pomijać bez straty znaczenia 
(wyjątek w $[2.4]). 

Na koniec, konwersji na i z postaci dziesiętnej możemy dokonywać za pośrednictwem opa- 
nowanej już przez nas postaci binarnej. 


2.4 System U2 reprezentacji liczb ze znakiem 


Interesującym jest pytanie, w jaki sposób można reprezentować w komputerze liczby ze zna- 
kiem, np. —10012? Jak widzimy, znak jest nowym (trzecim) symbolem, a elektroniczny prze- 
łącznik może znajdować się tylko w dwóch stanach. 

Spośród kilku możliwości rozwiązania tego problemu, obecnie w większości komputerów 
osobistych używany jest tzw. kod uzupełnień do dwóch, czyli kod U2. Niewątpliwą jego do- 
datkową zaletą jest to, iż pozwala na bardzo łatwą realizację operacji dodawania i odejmowania 
(zob. ćwiczenia). 

W tej notacji najstarszy bit determinuje znak liczby. I tak: 


— najstarszy bit równy O oznacza liczbę nieujemną, a z kolei 
— najstarszy bit równy 1 liczbę ujemną. 
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Jeśli dana jest n cyfrowa liczba w systemie U2 postaci bqn_1bn-2. ..b1bg, gdzie b; € {0,1} 
dlai=0,1,...,n— 1, to jej wartość wyznaczamy następująco: 
n—2 
—bm-1 X gril + y b; X de 
i=0 
Zatem dla liczby n bitowej najstarszy bit mnożymy nie przez 27", lecz przez —2"”!, 
Dla przykładu, 01102 = 1103 oraz 


10lly = —1 x 23 +0 x 2? +1 x2!+1 x2 = —5i = —101>. 


Co ważne, tutaj najstarszego bitu równego 0 nie można dowolnie pomijać w zapisie. Mamy 
bowiem 0b„—2. . . bou2 Æ bn—2 . - . bou2. Jednakże, 0b„—2 . . boya = 00...Obn_2...bov2. Z dru- 
giej strony, można udowodnić, że również 1b„-2...bou2 = 11... lbp_2...bove. 

Łatwo pokazać, że użycie k bitów pozwala na zapis liczb —2*7!,...,2*=1 — 1, np. liczby 8 
bitowe w systemie U2 mogą mieć wartości od —128 do 127. 


3 Reprezentacje liczb rzeczywistych 


Ostatnim ważnym dla nas zagadnieniem jest problem reprezentacji liczb rzeczywistych, a wła- 
Ściwie odpowiedniego jego podzbioru, przydatnego podczas dokonywania obliczeń. Rozwa- 
żymy tutaj tzw. reprezentację stałoprzecinkową i zmiennoprzecinkową. Ta ostatnio jest po- 
wszechnie używana we współczesnych komputerach. 


3.1 System stałoprzecinkowy 


W reprezentacji stałoprzecinkowej liczb rzeczywistych w systemie dwójkowym liczba bitów 
używanych do zapisu części „dziesiętnej” i części ułamkowej jest z góry ustalona. 

Rozpatrzmy liczbę postaci by_1bn_2...by,bę_1...bg. Jest to liczba n bitowa, w której na 
część ułamkową przypada k bitów. Jej wartość można wyznaczyć ze wzoru 


n—-1 
5 bi X 2%. 
1=0 


Tym samym np. 1011,1012 =8+2+14 L ł i = 11 10. Niestety, ustalenie liczby k z góry 
uniemożliwia przybliżanie liczb o (relatywnie) dużej i małej wartości z zadowalającą dokład- 


nością. 


3.2 System zmiennoprzecinkowy 


Współczesne komputery osobiste są w stanie przechowywać 1 przetwarzać liczby o długości 

dokładnie 8 (tzw. bajt, ang. byte), 16,32,64 lub nawet 128 bitów. Własność ta umożliwia zdefi- 

niowanie postaci zmiennopozycyjnej liczb rzeczywistych (ang. floating-point numbers). 
Reprezentowana liczba jest przechowywana w pamięci wg następującego schematu: 


KZEXMAZZ, 


gdzie: 
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— s € {0,1} — znak (1 bit, interpretowany jako —1>), 
— m — mantysa (odpowiednio znormalizowana), 
— e — wykładnik. 


Liczba cyfr mantysy i wykładnika jest ustalona z góry, np. w komputerach osobistych jest 
to 23+8 (32 bitowa) bądź 52+11 (64 bitowa liczba zmiennoprzecinkowa). 

Dokładne omówienie zagadnienia reprezentacji i arytmetyki liczb zmiennopozycyjnych wy- 
kracza poza ramy naszego wykładu. Więcej szczegółów, w tym rachunek błędów arytmetyki 
tych liczb, przedstawiony będzie na wykładzie z metod numerycznych] 

Warto jednak zapamiętać kilka problemów związanych z użyciem liczb zmiennopozycyj- 
nych. 


a) Nie da się reprezentować całego zbioru liczb rzeczywistych, a tylko jego podzbiór (a 
nawet właściwie jest to podzbiór liczb wymiernych!). Wszystkie liczby są zaokrąglane 
do najbliższej reprezentowalnej. 

b) Arytmetyka zmiennoprzecinkowa nie spełnia w pewnych szczególnych przypadkach wła- 
sności łączności i rozdzielności. 

c) Dodatkowo, w tym systemie zdefiniowane są trzy specjalne elementy: 


— oo (Inf, ang. infinity), 
— —oo(-Inf), 
— nie-liczba (NaN, ang. not a number). 


Np. 1/0 = Inf, 0/0 = NaN. 


Zainteresowanym osobom można polecić angielskojęzyczny artykuł: D. Goldberg, What Every Computer 
Scientist Should Know About Floating-Point Arithmetic, ACM Computing Surveys 21(1), 1991, 5—48. Do pobrania 
ze strony http://dlec.sun.com/pdf/800-7895/800-7895.pdf. 
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4 Cwiczenia 


Zadanie 2.1. Przedstaw następujące liczby całkowite nieujemne w postaci binarnej, dziesiętnej 
i szesnastkowej: 1010. 1810. 18:16, 101», 10150, 10156, ABCDEF'6, 6413531210, 1101011110012, 
FFFFFOC2;yg. 


Zadanie 2.2. Dana jest n cyfrowa liczba w systemie U2 zapisana jako ciąg cyfr bn—1bn—2 . . . b1bo, 
gdzie b, € 40,1) dlaż = 0,1,...,n — 1. Pokaż, że powielenie pierwszej cyfry dowolną liczbę 
razy nie zmienia wartości danej liczby, tzn. bq_1bn_2... bybg = bn-1bdn-1 - - - On—10n—2 . - . b1bo- 


Zadanie 2.3. Przedstaw następujące liczby dane w notacji U2 jako liczby dziesiętne: 10111100, 
00111001, 1000000111001011. 


Zadanie 2.4. Przedstaw następujące liczby dziesiętne w notacji U2 (do zapisu użyj 8, 16 lub 
32 bitów): —12, 54, —128, —129, 53263, —32000, —56321, —3263411. 


x Zadanie 2.5. Rozważmy operacje dodawania i odejmowania liczb nieujemnych w reprezenta- 
cji binarnej. Oblicz wartość następujących wyrażeń korzystając z metody analogicznej do spo- 
sobu „szkolnego” (dodawania i odejmowanie słupkami). Pamiętaj jednak, że np. 12 + 12 = 102. 


a) 1000, + 111», e) 1111, — 00015, 
b) 1000, + 1111», f) 1000 — 0001, 

c) 1110110» + 11100111», g) 10101010, — 11101, 

d) 111115 + 11111111, h) 11001101» — 100101110. 


x Zadanie 2.6. Okazuje się, że liczby w systemie U2 można dodawać i odejmować tą samą 
metodą, co liczby nieujemne w reprezentacji binarnej. Oblicz wartość następujących wyrażeń 
i sprawdź otrzymane wyniki, przekształcając je do postaci dziesiętnej. Uwaga: operacji dokonuj 
na dwóch liczbach n bitowych, a wynik podaj również jako n bitowy. 


a) 1000y2 + O111yv>, d) 11001101y2 — 10010111vo. 
b) 10110110y2 + 11100111v2, e) 10001101y2 — 01010111yvo. 
c) 10101010y2 — 00011101yz, 


x Zadanie 2.7. Niech dana będzie liczba w postaci U2. Pokaż, że aby uzyskać liczbę do niej 
przeciwną, należy odwrócić wartości jej bitów (dokonać ich negacji, tzn. zamienić zera na je- 
dynki i odwrotnie) i dodać do wyniku wartość 1. Ile wynoszą wartości O101yva, 1001y2 i 0l1llu2 
po dokonaniu tych operacji? Sprawdź uzyskane rezultaty za pomocą konwersji tych liczb do 
systemu dziesiętnego. 
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5 Wskazówki do ćwiczeń 


Odpowiedź do zadania |2.1 


a) 1010 = 1010; = Aig. 

b) 1819 = 10010 = 1216. 

c) 1816 = 11000, = 2410. 

d) 1012 = 510 = 546. 

e) 10110 = 1100101 = 6516. 

f) 10116 = 1000000012 = 25710. 

g) ABCDEF;6 = 101010111100110111101111 = 1125937510. 
h) 6413531210 = 11110100101010000010010000> = 3D2A09055. 
i) 110101111001, = D79,6 = 344910. 


jJ) FFFFFOC2;6 = 11111111111111111111000011000010> = 42949633949. 


Odpowiedź do zadania 2.3] 


a) 10111100y2 = —68. 
b) 00111001y2 = 5710. 
c) 1000000111001011y2 = —32309. 


Odpowiedź do zadania 2.4] 


a) —1210 = 11110100v2. 

b) 5410 = 001101104. 

c) —12810 = 10000000». 

100 = Ta 

e) 5326310 = 00000000000000001101000000001111v>. 

f) —3200010 = 11111111111111111000001100000000v>. 
g) —563219 = 11111111111111110010001111111111y». 
h) —3263411; = 11111111110011100011010001001101y>. 


Odpowiedź do zadania 2.5] 


a) 1000 + 111, = 11117. 
b) 10002 + 1111, = 101112. 
c) 11101103 + 111001112 = 1010111012, gdyż 


+ 1 
= 1 0 


O|Ę = 
elo = 
—|OC- 
=|= = = 
O|- = 
—|- © 


d) 111112 + 11111111 = 100011110». 
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e) 
f) 
g) 


Odpowiedź do zadania |2.6 
a) 1000y2 + 0111y2 = 1111y2 = —110 = —810 + 010. 
b) 
c) 10101010v2 — 00011101y2 = 10001101y2 = —11510 = —8610 — 2910. 
d) 


e) 


1111, — 00012 = 1110. 


1000, — 0001, = 111». 
101010103 — 11101, = 10001101», gdyż 
1 -1 =l 
1 0 1 1 0 1 
: 110 
= 1000110 
11001101, — 10010111, = 110110», gdyż 
-1 1 -1 -1 
1100110 
-10010 11 
=0011011l 


10110110y2+11100111y2 = 110011101y2 = 10011101y2 = —9910 = —7410+(—25)i0. 


11001101v2 — 10010111y2 = 00110110y2 = 5410 = —5110 — (—10510). 


10001101y2—01010111y2 = 00110110y2 = 5410 £ —11510—8710 = —20240! Wystąpiło 
tzw. przepełnienie, czyli przekroczenie zakresu wartości dla liczby 8-bitowej (ale mamy 
28 — 20210 = 5410). Można zapisać operandy na większej liczbie bitów, np. szesnastu: 


1111111110001101y2 — 0000000001010111y2 = 1111111100110110y2 = —20210. 


Aby lepiej zrozumieć to zjawisko, rozważ „kolejne” 4 bitowe liczby w systemie U2. 


U2 DEC 
0000 0 
0001 1 
0010 2 
0011 3 
0100 4 
0101 5 
0110 6 
0111 7 
1000 -8 
1001 -7 
1010 -6 
1011 -5 
1100 -4 
1101 -3 
1110 -2 
1111 -1 
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1 Deklaracja zmiennych w C++ 


W poprzednich rozdziałach, czytając pseudokody algorytmów, często napotykaliśmy na in- 
strukcję podobną do następującej. 


niech xeR; 


Tym samym mieliśmy na myśli „niech od tej pory x będzie zmienną ze zbioru liczb rzeczywi- 
stych”. 

Instrukcja tego typu to tzw. deklaracja zmiennej. Zmienna służy do przechowywania do- 
wolnych wartości ze zbioru, z którego pochodzi. Taki zbiór nazywamy typem zmiennej (w przy- 
padku powyższego x było to R). 

Mamy dwie możliwości użycia zmiennej: możemy jej przypisać wartość bądź przecho- 


wywaną weń wartość odczytać. W tym sensie zmienną można porównać do szuflady, która 
ma swoją etykietkę (nazwę, identyfikator), jednoznacznie ja identyfikujaca. Do takiej szuflady 
można coś włożyć (jednocześnie pozbywając się wcześniejszej zawartości) bądź coś z niej wy- 
jaé. 

Dzięki zmiennym pisane przez nas programy (ale także np. twierdzenia matematyczne) nie 
opisują jednego, konkretnego przypadku szczególnego pewnego interesującego nas problemu, 
lecz wszystkie możliwe w danym kontekście. Dzięki takiemu mechanizmowi mamy więc moż- 
liwość uogólniania naszych wyników. Przykładowo, ciąg twierdzeń 


Twierdzenie 1. Pole koła o promieniu 1 wynosi r. 
Twierdzenie 2. Pole koła o promieniu 2 wynosi 4r. 
Twierdzenie 3. ... 


przez każdego matematyka (nawet przyszłego matematyka) byłby uznany, łagodnie rzecz uj- 
mując, za nieudany żart. Jednakże uogólnienie powyższych wyników przez odwołanie się do 
de facto zadeklarowanej zmiennej sprawia, że wynik staje się interesujący. 


Twierdzenie 4. Niech r € R,. Pole koła o promieniu r wynosi mr”. 


Zauważmy, że w matematyce zdanie deklarujące zmienną „Niech r € R...” nie musi się 
pojawiać w takiej formie. Większość autorów podchodzi do tego, rzecz jasna, w sposób ela- 
styczny. Mimo to, pisząc np. „Pole koła o promieniu r > 0 wynosi mr*.” bądź nawet „Pole koła 
o promieniu r wynosi mr*.”, domyślamy się, że chodzi właśnie o powyższą konstrukcję. 

Wiemy już, że komputery jednak nie są tak domyślne jak my. W języku C++ wszystkie 
zmienne muszą zostać zadeklarowane przed ich pierwszym użyciem. Składnia najprostszej 
instrukcji służącej do osiągnięcia tego celu jest ściśle określona: 


typ ident; 


gdzie typ jest typem zmiennej ($[1.1), a ident jej identyfikatorem ($[1.2). 

Deklaracja zmiennej powoduje przydzielenie jej pewnego miejsca w pamięci RAM kom- 
putera (więcej szczegółów na ten temat na wykładzie o wskaźnikach). Warto zapamiętać, że 
zmienna nie jest inicjowana automatycznie. Dopóki nie przypiszemy jej jawnie jakiejś warto- 
ści, będzie przechowywać ,,$mieci”. 


AiPP-III, s. 2. 


Uwaga 


Każda instrukcja w C++ musi być zakończona średnikiem! 


1.1 Typy liczbowe proste 

W niniejszym paragrafie omówimy dostępne w języku C++ typy zmiennych liczbowych (są to 
tzw. typy proste). 

1.1.1 Typy całkowite 


Wśród typów całkowitych wyróżniamy: 


Nazwa typu | Liczba bitów! 


char 8 bitów 
short 16 bitów 
int 32 bity 


long long | 64 bity 


My najczęściej będziemy korzystać z typu int. Jest to liczba 32-bitowa w systemie U2. 
Zatem może być ona traktowana jako zbiór int = [-2%, —23t + 1,...,231 — 1} CZ. 

Co więcej, do każdego z ww. typów można zastosować modyfikator unsigned — wtedy 
otrzymujemy liczbę nieujemną (w zwykłym systemie binarnym). np. unsigned int = (0,1,..., 
232 — 1). 

Stałe całkowite, np. 0, —5, 15210, są reprezentantami typu int danymi w postaci dzie- 
siętnej. Co więcej, można używać stałych w postaci szesnastkowej, poprzedzając liczby przed- 
rostkiem Ox, np. 0x53a353 czy też Oxabcdef, oraz ósemkowej, korzystając z przedrostka 0, np. 
0777. Nie ma, niestety, możliwości definiowania bezpośrednio stałych w postaci dwójkowej. 


1.1.2 Typy zmiennoprzecinkowe 


Z kolei wśród typów zmiennoprzecinkowych wyróżniamy: 


Nazwa typu | Liczba bitów! 


float 32 bity 
double 64 bity 


Mamy] float C R, double c R. W przypadku liczby 32 bitowej najmniejsza reprezen- 
towalna wartość dodatnia to ok. 1,18 x 107%, a największa — ok. 3,4 x 10%, Dla liczby 64 
bitowej mamy odpowiednio ok. 2,23 x 10 3% i 1,80 x 10308, 

Stałe zmiennoprzecinkowe wprowadzamy podając zawsze część dziesiętną i część ułam- 
kową rozdzieloną kropką (nawet gdy część ułamkowa jest równa 0), np. 3.14159, 1.0 albo 
—0.000001. Domyślnie należą one do typu double. Dodatkowo, liczby takie można wprowa- 
dzać w formie notacji naukowej, używając do tego celu separatora e, który oznacza „razy 
dziesięć do potęgi”. I tak liczba —2.32e—4 jest stałą o wartości —2,32 x 1074. 


!Podana liczba bitów dotyczy Microsoft Visual C++. 
2Przypomnijmy, że typy zmiennopozycyjne przechowują też wartości specjalne: Inf, -Inf oraz NaN. 
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1.1.3 Typ logiczny 


Oprócz powyższych, wśród typów prostych, dostępny jest jeszcze typ bool, służący do repre- 
zentowania zbioru wartości logicznych. 

Zmienna tego typu może znajdować się w dwóch stanach, określonych przez stałe logiczne 
true (prawda) oraz false (fałsz). 


1.2 Identyfikatory 


Każda zmienna musi mieć jakąś nazwę, aby można się było do niej odwołać. Identyfikatory 
mogą się składać z następujących znaków: 


abcde fghijklmnopqrstuvwxy z 
ABCD EF GH I J KLMNOPOQORSTUVWKY Z 


oraz z poniższych, pod warunkiem, że nie są one pierwszym znakiem identyfikatora: 
0123456 7 8 9 


Zatem przykładami poprawnych identyfikatorów są: i, suma, wyrazenieZeWzoru32, _2mla 
oraz zmienna_pomocniczal. Z kolei zgodnie z podaną regułą np. 3523aaa, ala ma kota oraz 
:—) nie mogą być identyfikatorami. 

Identyfikatorami nie mogą być też zarezerwowane słowa kluczowe (np. int, double). Dla- 
tego specjalnie wyróżniamy je w tym tekście pogrubioną czcionką. 


Uwaga 
W języku C++ wielkość liter ma znaczenie! Np. 7 oraz i to identyfikatory dwóch różnych 
zmiennych. Z tego też powodu identyfikator Int nie jest tożsamy z typem int. 


Na koniec naszych rozważań o identyfikatorach pewna rada praktyczna. Dobrą praktyką jest 
nadawanie takich nazw zmiennym, by w pewnym sensie były samoobjaśniające się. Zmienne 
poleKwadratu, delta, wspX w pewnych kontekstach posiadają tę cechę. Często używamy też 
identyfikatorów i, j, k jako liczników w pętlach. Z, drugiej strony, np. zmienna albo ratunku 
niezbyt jasno mówią do czego same mogą służyć. 


2 Operatory 


2.1 Operator przypisania 


Operator przypisania, =, służy do nadania zmiennej wartości. Składnia: 


ident = wyr; 


gdzie ident to identyfikator zmiennej, której wartość chcemy ustalić, a wyr to wyrażenie (ang. 
expression), którego wykonanie (obliczenie, ewaluacja) zwraca pewną wartość odpowiedniego 
typu. Wyrażeniami mogą być m.in. 


— stała całkowita, zmiennoprzecinkowa, logiczna (np. 7, —4.le7, false ), 
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— identyfikator zmiennej bądź stałej symbolicznej (wartość zostanie z niej pobrana), 

— wartość zwracana przez funkcję (więcej o funkcjach w kolejnych wykładach), 

— wynik zastosowania operatorów arytmetycznych, logicznych, relacyjnych lub bitowych 
(zob. $P.3]i dalej) na innych wyrażeniact | 


Prześledźmy poniższy kod. 


int i; // deklaracja zmiennej 

int j; 

i = 4; // przypisanie wartości 4, czytaj: i staje się 4 
J = i; /* przypisanie wartości takiej , 


jaką aktualnie przechowuje zmienna i */ 


cout << j; // wypisze na ekran "4" 


Uwaga 


Tekst drukowany kursywą i szarą czcionką to komentarze. Komentarze są albo zawarte 
między znakami /* i */ (i wtedy mogą zajmować wiele wierszy kodu), albo następują po znakach 
// (i wtedy sięgają do końca aktualnego wiersza). 

Treść komentarzy jest ignorowana przez kompilator i nie ma wpływu na wykonanie pro- 
gramu. Jednakże opisywanie tworzonego kodu jest bardzo dobrym nawykiem, bo sprawia, że 
jest on bardziej zrozumiały dla jego czytelnika. 


Uwaga 


W języku C++ istnieje specjalny obiekt o nazwie cout, umożliwiający wypisywanie war- 
tości zmiennych różnych typów na ekran monitora (a dokładnie, na tzw. standardowe wyjście). 
Drukowaną zmienną „wysyła” się do niego za pomocą operatora <<. 


Istnieje jeszcze inny sposób przypisywania wartości zmiennym, dzięki któremu dane nie 
muszą być określone podczas pisania kodu. Informacje te można określić w trakcie działania 
programu. Można o nie bowiem poprosić użytkownika, aby je podał za pomocą klawiatury (za 
pośrednictwem obiektu cin i operatora >>). 


int x; 
cout << "Podaj wartość x: "; 
cin >> x; /* wczytanie wartości x z tzw. standardowego wejścia 


(domyślnie jest to klawiatura komputera) */ 
cout << "Teraz x=" << x << endl; 


U 


Przy okazji omawiania operatora przypisania, warto wspomnieć o możliwości definiowania 
stałych symbolicznych, według składni: 


const typ ident = wartość; 


3Zauważmy, że definicja wyrażenia jest rekurencyjna; działania na wyrażeniach też są wyrażeniami. 
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Przy tworzeniu stałej symbolicznej należy od razu przypisać jej wartość. Cechą zasadniczą 
stałych jest to, że raz zainicjowane, nie mogą być później zmienione (warto je stosować, by 
wyeliminować możliwość pomyłki programisty). 

Przykłady: 


N 


w 


ÉS 


const int dlugoscCiagu = 10; 

const double pi = 3.14159; 

const double G = 6.67428e—11; /* to samo, co 6,67428 x 10711. x/ 
pi = 4.3623; // błąd — pi jest "tylko do odczytu" 


N 


2.2 Operatory rzutowania. Hierarchia typów 


Zastanówmy się, co się stanie, jeśli zechcemy przypisać zmiennej jednego typu wartość innego 
typu. Rozważmy inny przykład: 


double x; // deklaracja zmiennej 
int i = 4, j; /* deklaracja dwóch zmiennych , 
zmienna i zostanie zainicjowana wartością 4, 


zmienna j jest niezainicjowana */ 


== E, /x inny typ! niejawna konwersja na 4.0, 


tzw. promocja */ 


x = 3.14159; 
RZE: // błąd! (typ docelowy jest słabszy) 
j = (int)x; // jawna konwersja (tzw. rzutowanie) — teraz OK 


" " 


cout << J <<", << x; // wypisze na ekran "3, 3.14159" 


Przypisanie wartości typu double do zmiennej typu int zakończy się błędem kompilacji. 
Jest to spowodowane tym, że w tym miejscu mogłaby wystąpić utrata informacji. W języku 
C++ obowiązuje następująca hierarchia typów: 


typ logiczny typy całkowite typy zmiennoprzecinkowe 


rwn AAA A 2 2 AAA, AAA AAA 
bool C char C short C int C long long € float C double . 


Promocja typu, czyli konwersja na typ silniejszy, bardziej „pojemny”, odbywa się automa- 
tycznie. Dlatego w powyższym kodzie instrukcja x=i; wykona się poprawnie. 

Z, drugiej strony, rzutowanie do słabszego typu musi być zawsze jawne! Należy expli- 
citć powiedzieć kompilatorowi, że postępujemy Świadomie. Służy do tego tzw. operator rzu- 
towania. Taki operator umieszczamy przed rzutowanym wyrażeniem. Oprócz zastosowanej 
w powyższym przykładzie składni (typ)wyrażenie, dostępne są dodatkowo dwie równoważne: 
typ(wyrażenie ) oraz static_cast <typ>(wyrazenie). 
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Warto zapamiętać, że w przypadku rzutowania typu zmiennopozycyjnego do całkowitolicz- 
bowego następuje obcięcie części ułamkowej, np. (int )3.14 da w wyniku 3, a (int )—3.14 zaś 
wartość —3. 

Z, kolei dla konwersji do typu (bool), wartość równa O da zawsze wynik false, a różna od 
0 będzie przekształcona na true. 


Uwaga 


Przy konwersji typu bool do całkowitoliczbowego wartość true zostaje zamieniona zawsze 
na 1, a false na O. 


2.3 Operatory arytmetyczne 


Wśród operatorów arytmetycznych można wyróżnić operatory binarne (wymagające dwóch 
operandów) oraz unarne (przyjmujące jeden operand). Można je (poza wyszczególnionymi wy- 
jątkami) stosować dla wyrażeń typu całkowitego i zmiennoprzecinkowego. 


Dostępne operatory binarne przedstawia poniższa tabela. 


Operator Znaczenie 
+ dodawanie 
= odejmowanie 
* mnożenie 
dzielenie 
% reszta z dzielenia (tylko całkowite) 


Uwaga 


W C++ nie ma operatora potęgowania! 


Przykłady: 


cout << 2+2 << endl; 

cout << 3.0/2.0 << endl; 

cout << 3/2 << endl; 

cout << T%4 << endl; 

cout << 1/0.0 << endl; /x uzgodnienie typów — to samo, co 
1.0/0.0 x/ 


Jak widzimy, zastosowanie operatorów do wyrażeń dwóch różnych typów powoduje konwersję 
operandu o typie słabszym do typu silniejszego operandu. Wynikiem wykonania rozpatrywanej 
części kodu będzie: 


(dlaczego?) 
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Zajmijmy się teraz operatorami unarnymi. 


Operator Znaczenie 


N 


+ unarny plus (nic nie robi) 
— unarny minus (zmiana znaku na przeciwny) 
Przykład: 
int i = —4; // zmiana znaku stałej 4 
i =<i; // zmiana znaku zmiennej i 
cout << i; // wynik: 4 


Ponadto zdefiniowane są jeszcze dwa operatory unarne, które można zastosować tylko do 
zmiennych typu całkowitego. Operator inkrementacji (++) służy do zwiększania wartości zmien- 
nej o 1, a operator dekrementacji (——) zmniejszania o 1. 

Zdefiniowane są aż dwa warianty powyższych operatorów: przedrostkowy (ang. prefix) oraz 
przyrostkowy (ang. suffix), co — nie wiedzieć czemu — powoduje rokrocznie przerażenie stu- 


denckiej braci. 


A przecież zasada jest bardzo prosta! W przypadku operatora inkrementacji/dekrementacji 
przedrostkowego, wartością całego wyrażenia jest wartość zmienionej zmiennej (operator ten 
działa tak, jakby najpierw aktualizował stan zmiennej, a potem zwracał swoją wartość w wy- 
niku). Z kolei dla operatora przyrostkowego wartością wyrażenia jest stan zmiennej sprzed 
zmiany (najpierw zwraca wartość zmiennej, potem aktualizuje jej stan). 


Poświęćmy kilka chwil na przeanalizowanie następującego przykładu. 


int j = 5, j = 9; 


=j; // wynik: 
j——; // 
j = ++i; // wynik: 
j = i++; // 
j = ++100; // błąd, 


100 


j=l; 


to samo to j = j—l; 


to samo to j 


j==6 (przedrostkowy) 


i==/ (przyrostkowy) 


jest stałą 


Zatem j = ++i; można zapisać jako dwie instrukcje: i = i+l; j =i; 
Z, drugiej strony, j = i++; jest równoważne operacjom: j = i; i =i+l; 


Niewątpliwą zaletą tych dwóch operatorów jest możliwość napisania zwięzłego fragmentu 
kodu. Niestety, jak widać, odbywa się to kosztem czytelności. 


2.4 Operatory relacyjne 


Operatory relacyjne służą do porównywania wartości innych wyrażeń. Wynikiem takiego po- 
równania jest wynik typu bool. Są to operatory binarne. 
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Operator Znaczenie 


== czy równe? 
l= czy różne? 


< czy mniejsze? 

<= czy nie większe? 
> czy większe? 

>= czy nie mniejsze? 


Nie należy mylić operatora porównania (==) i operatora przypisania (=) — jest to częsty 
błąd. 


Przykłady: 


bool w/ = l==1; // true 
int w2 = 2>=3; // 0, czyli (int)false 
bool w3 = (0.0 == (((1e34 + le-34) — le34) — le—34)); // false! 


// Uwaga na porównywanie liczb zmiennoprzecinkowych — błędy! 


2.5 Operatory logiczne 


Operatory logiczne służą do działań na wyrażeniach typu bool (lub takich, które są sprowa- 
dzalne do bool). 


Operator Znaczenie 


! negacja (unarny) 
Il alternatywa (lub) 
&& koniunkcja (i) 


Wyniki powyższych działań na wszystkich możliwych wartościach operandów zestawione 
są w poniższych tzw. tablicach prawdy. 


! Il false | true SSL false | true 
false | true false | false | true false || false | false 
true false true true | true true false | true 

Kilka przyktadów: 
bool w/ = (!true); // false 
bool w2 = (1<2 II 0<1); // true —> (true || false) 
bool w3 = (1.0 ££ 0.0); // false —> (true && false) 


Uwaga 


Okazuje się, że komputery są czasem leniwe. W przypadku wyrażeń składających się 
z wielu podwyrażeń logicznych, obliczane jest tylko to, co jest potrzebne do ustalenia wyniku. 
I tak w przypadku koniunkcji, jeśli jeden operand ma wartość false, to wyznaczanie wartości 
drugiego jest niepotrzebne. Podobnie jest dla alternatywy 1 jednego operandu prawdziwego. Na 
przykład: 
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int a = 0; 


bool p = (true || (++a==0)); // leniwy 

cout << a; // a==0 

bool q = (true && (++a==0)); // tutaj już musi policzyć 
cout << a; IL == 


2.6 Operatory bitowe 


Operatory logiczne działają na poszczególnych bitach operandów typu całkowitego (jak już 
zdążyliśmy się przekonać, każda informacja, w tym i liczby, są przechowywane w pamięci 
komputera jako ciąg bitów). Zwracany wynik jest też liczbą całkowitą. 


Operator Znaczenie 


= bitowa negacja (unarny) 
l bitowa alternatywa 
& bitowa koniunkcja 
di bitowa alternatywa wyłączająca (albo, ang. exclusive-or) 
<<k przesunięcie w lewo o k bitów (ang. shift-left) 
>> k przesunięcie w prawo o k bitów (ang. shift-right) 


Operator bitowej negacji zamienia każdy bit liczby na przeciwny. Operatory bitowej alter- 
natywy, koniunkcji i alternatywy wyłączającej zestawiają bity na odpowiadających sobie pozy- 
cjach dwóch operandów, według następujących reguł. 


~ IJO/1) [&O0O 1) [2 [0 1 


O0J1) [0 0|1| |0 
1/0) |1)/1)1) |1]0/1]| |1010 


© 
© 
© 
— 


Np. 0xb6 ^ 0x5f == Oxe9, gdyż 


Operatory przesunięcia zmieniają pozycje bitów w liczbie. Bity, które nie „mieszczą się” 
są tracone. W przypadku operatora <<, na zwalnianych pozycjach wstawiane są zera. Operator 
>> wstawia zaś na zwalnianych pozycjach bit znaku. Dla przykładu, dokonajmy przesunięcia 
bitów liczby 8 bitowej o 3 pozycje w lewo. 


wynik 
10101011 101 01011000 
Oxab <<3 == 0x58 


Oto kolejne przykłady: 
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N 


int pl = 1 | 2; // 3, bo 00000001 | 00000010 —> 00000011 


int p2 = -0xfO; // 15, bo ~11110000 => 00001111 
int p3 = Oxb ^ 5; // 14, bo 00001011 ^ 000000101-—> 00001110 
int p4=5 << 2; // 20, bo 00000101 << 2 => 00010100 
int pó= 75 1; // 3, bo 00000111 >> 1 => 00000011 
int pó = —128 >> 4;// —8, bo 10000000 >> 4 => III11000 


Uwaga: n<<k równoważne jest (o ile nie nastąpi przepełnienie) n x 2*, a n>>k równoważne 
jest całkowitoliczbowej operacji n x 2* (k krotne dzielenie całkowite przez 2). 


Uwaga 


Zauważmy, że operatory << i >> mają inne znaczenie w przypadku użycia ich na obiek- 
tach, odpowiednio, cout i cin. Jest to praktyczny przykład zastosowania tzw. przeciążania ope- 
ratorów, czyli zróżnicowania ich znaczenia w zależności od kontekstu. Więcej na ten temat 
dowiemy się na II semestrze. 


2.7 Operatory łączone 


Dla skrócenia zapisu, zostały też wprowadzone tzw. operatory łączone, łączące operacje aryt- 
metyczne lub bitowe i przypisanie: 


+=, —=, *=, f=: %=, > =, ne, <<=, >: 


Np. x += 10; znaczy to samo, cox = x + 10; 


2.8 Priorytety operatorów 


W niektórych z powyższych przykładów zaobserwowaliśmy, że w jednym wyrażeniu można 
używać wielu operatorów. Kolejność wykonywania działań jest jednak Ściśle określona. Można 
ją wymusić za pomocą nawiasów okrągłych (...). 

Domyślnie jednak wyrażenia obliczane są według priorytetów, zestawionych w kolejności 
od największego do najmniejszego w tab.|1| W przypadku operatorów o tym samym priorytecie, 
operacje wykonywane są w kolejności od lewej do prawej. 


Przykłady: 
int pl = 2+4x3/2; // pl = (2+((4x3)/2]); 
pl += 2>3 == |true; /* pl += ((2>3) == (/true)); 
czyli to samo, co pl++; */ 
cout << pl; // 9 
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Tablica 1: Priorytety operatorów. 


13 | ++, —— (przyrostkowe) 

12 | ++, —— (przedrostkowe), +, — (unarne), !, - 
11 | x, /, % 

10 | +, — 

9 | <<, >> 

8 | <, <=, >, >= 

7 ==, l= 

6 | © 

5 A 

4 1 

3 | ZE 

2/1 

| |=, +=, —=,*=, |=, %=, |=, &=, ^=, <<=, >>= 
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3 Cwiczenia 


Zadanie 3.1. Wyraź następujące liczby zmiennoprzecinkowe dane w notacji naukowej w zwy- 
kłej postaci dziesiętnej. 


a) 42166, c) 2.314e—4, 
b) 1.95323e2, d) 4.235532e—2. 


Zadanie 3.2. Wyznacz wartość następujących wyrażeń. Ponadto określ typ zwracany przez 
każde z nich. 


a) 10.0+15.0/2+4.3, P 10+17/3.0+4, 
b) 10.0+15/2+4.3, g) 3*4%06+6, 

c) 3.0*4/6+6, h) 3.0x4%6+6, 
d) 20.0—2/6+3, 1) 10+17%3+4. 


e) 10+17x3+4, 


Zadanie 3.3. Dane są 3 zmienne zadeklarowane w sposób następujący. 


double a = 10.6, b = 13.9, c = —3.42; 


Oblicz wartość poniższych wyrażeń. 


a) int(a), e) int(a+b)xc, 
b) int(c), f) double(int(a+b)/int(c)), 
c) int(a+b), g) double(int(a) )/c. 


d) int(a)+b+c, 


Zadanie 3.4. Korzystając z przekształceń logicznych (np. praw De Morgana, praw rozdzielno- 
$ci) uprość następujące wyrażenia (zakładamy, że a,b,c,d są typu double, a p,q,r typu bool). 


a) Up), e) !l(a>=b && b>=c && a>=C), 
b) !p śe !ą, ) (a>b && a<c) Il (a<c && a>d), 
c) (p Il tą Il tr), 9 p II !p. 


d) !(b>a && b<c), 


Zadanie 3.5. Jaki wynik dadzą następujące operacje bitowe wykonane na danych typu short? 


a) OxOFCD |OxFFFF, e) 14 <<4, 
b) 364 & 01323, f) 0xf5a3 >>8, 
c) -163, g) 0x3a9f >> 7. 


d) 0xFC93 ^ 0x201D, 


Zadanie 3.6. Niech dana będzie zmienna typu int. W jaki sposób dokonać zmiany znaku war- 
tości tej zmiennej na przeciwny nie używając operatora —? 
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Zadanie 3.7. Policyjny fotoradar emituje fale elektromagnetyczne o częstotliwości f, Hz. Wią- 
zka tych fal jest odbijana od nadjeżdżającego z prędkością v samochodu i, jako że auto jest w ru- 
chu, powraca do urządzenia ze zmienioną częstotliwością f, Hz. Polscy funkcjonariusze testują 
właśnie najnowszy model brytyjskiej „suszarki”. Związek pomiędzy omawianymi zmiennymi 
można wyrazić równaniem 
Í o f e 
fat Je 
Prędkość jednak jest podawana w milach na godzinę. Wiedząc, że 1 mila to 1609,344 m, napisz 
fragment kodu w języku C++, który dla danego fei f, poda prędkość nadjeżdżającego samo- 
chodu w km/h. Ponadto wypisz na ekran wartość logiczną mówiącą, czy została przekroczona 
dopuszczalna prędkość, wynosząca w tym miejscu 50 km/h. 

Jaki będzie wynik działania tej procedury dla f, = 2 x 107” Hz i f, = 2.0000004 x 1079 
Hz? 


v = 6,685 x 108 


x Zadanie 3.8. Danych jest 6 zmiennych ax,ay,bx,by,cx,cy typu double, reprezentujących 
współrzędne 3 punktów w R?: a = (ax,ay), b = (bx,by), e = (cx,cy). Napisz fragment kodu, 
który wyznaczy kwadrat promienia okręgu przechodzącego przez a, b i c. Dany jest on wzorem 


2  la— c] |b — cf? ja — b|? 
- 4lla=c)x(b=c)2 ” 


r 


gdzie np. |a — b| = y (ax — bx)? + (ay — by)? oraz a x b = axby — aybx. 


Zadanie 3.9. Dane są a,b,c,d,e, f € R takie, że układ dwu równań liniowych względem 
niewiadomych x, y € R: 


ar+by = c, 
dr +ey = f. 


jest oznaczony. Zaproponuj fragment kodu w języku C++, który wyznaczy jego rozwiązanie, 
tzn. obliczy wartość zmiennych x,y na podstawie pewnych wartości zmiennych a,b,c,d,e, f. 
Do reprezentacji liczb rzeczywistych użyj typu double. 


AiPP-III, s. 14. 


4 Wskazówki do ćwiczeń 


Odpowiedź do zadania 3.2] 


a) 10.0+15.0/2+4.3 == 21.8 (double). 
b) 10.0+15/2+4.3 == 21.3 (double). 
c) 3.0x4/6+6 == 8.0 (double). 

d) 20.0—2/6+3 == 23.0 (double). 

e) 10+17*3+4 == 65 (int). 

f) 10+17/3.0+4 — 19.6667 (double). 
g) 3*4%6+6 == 6 (int). 

h) 3.0x4%6+6 — błąd kompilacji. 

1) 10+17%3+4 == 16 (int). 


Odpowiedź do zadania|3.5| 


a) OxOFCD | 0xFFF. 
b) 364 & 0x323 == 288 (0x0720). 

c) -163 == — 164 (OxFF5C) (dlaczego -164?). 
d) 0xFC93 ^ 0x201D == —9074 (0xDC8E). 

e) 14 << 4 == 224 (0x00E0). 

f) 0xf3a3 >> 8 == 0xFFFS (bit znaku == 1). 
g) 0x3a9f >> 7 == 117 (00075). 


== —1 (OxFFFP). 


Wskazówka do zadania B.6| Zmienna typu int jest reprezentowana w kodzie U2. 


Odpowiedź do zadania |3.8 


// Dane 
double 
double 
double 
// Dane 
double 


double 
double 
double 
double 
double 
double 
double 
double 


wejściowe: 


ax = 
bx 


cx = 


., ay = ...; // współrzędne punktu a. 
., by = ...; // współrzędne punktu b. 
«> CY = ...; // współrzędne punktu œ. 


wyjściowe: 


F2: 


amcx 
amcy 
bmcx 
bmcy 
a_c 
b_c 
a_b 


// 


kwadrat promienia okręgu przech. 


ax—cx; //a=C — wsp. x 
ay—cy; //a=C — wsp. y 
bx—cx; //b—c —— wsp. x 
by-cy; //b—c —— wsp. y 


amcx*xamcxtamcy*xamcy; //|a— e]? 


= bmcxxbmcx+bmcy*bmcy; //|b— el? 


ilwekt 


= (ax-bx)*(ax—bx)+(ay—by)*(ay—by); 


= amcx*bmcy — amcyx*bmcx; 


r2 = a _ cxb_cxa b/4.O/ilwekt/ilwekt:; 


przez a,b,c. 


//|a— bl? 
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1 Bezpośrednie następstwo instrukcji 
Poznaliśmy już następujące zagadnienia: 


— deklaracja zmiennych i stałych nazwanych, 

— przypisywanie wartości, 

— obliczanie wartości wyrażeń z użyciem operatorów, 

— wypisywanie wartości zmiennych na ekran i wczytywanie ich wartości z klawiatury. 


Zauważmy, że do tej pory nasz kod w języku C++ nie miał żadnych rozgałęzień. Wszyst- 
kie instrukcje były wykonywane jedna po drugiej. Owo tzw. bezpośrednie następstwo można 
zobrazować jak na rys. |l| Jest to tak zwany schemat blokowy algorytmu (ang. control flow 
diagram, czyli schemat przepływu sterowania). 


Rysunek 1: Schemat blokowy bezpośredniego następstwa instrukcji. 


W kolejnych paragrafach poznamy konstrukcje, które pozwolą nam zmieniać przebieg pro- 
gramu. Dzięki temu będziemy mogli dostosowywać jego działanie do różnych szczególnych 
przypadków w zależności od wartości przetwarzanych danych. 


2 Instrukcja warunkowa 


Podczas zmagania się z różnymi problemami prawie zawsze spotykamy sytuację, gdy musimy 
dokonać jakiegoś wyboru. Na przykład, rozwiązując równanie kwadratowe inaczej postępu- 
jemy, gdy ma ono dwa pierwiastki rzeczywiste, a inaczej, gdy nie ma ich wcale. Policjant 
mierzący fotoradarem prędkość nadjeżdżającego samochodu inaczej postąpi, gdy stwierdzi, że 
dopuszczalna prędkość została przekroczona o 50 km/h, niż gdyby się okazało, że kierowca je- 
dzie prawidłowo. Student przed sesją głowi się, czy przyłożyć się bardziej do programowania, 
algebry, do obu na raz (jedynie słuszna koncepcja) czy też dać sobie spokój i iść na imprezę itp. 

W języku C++ takim mechanizmem używanym do dokonywania wyborów jest instrukcja 
warunkowa if (od ang. jeśli). Wykonuje ona pewien fragment kodu wtedy i tylko wtedy, gdy 
pewien dany warunek logiczny jest spełniony. Jej składnia jest następująca. 


if (warunek) instrukcjaT; 
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gdzie warunek (koniecznie ujęty w nawiasy!) jest pewnym wyrażeniem sprowadzalnym do typu 
bool. Stosowny schemat blokowy przedstawia rys. 


Gz) 


Rysunek 2: Schemat blokowy instrukcji warunkowej if. 


Jeśli warunek nie jest spełniony, możliwe jest wykonanie innej instrukcji poprzez dodanie 
słowa kluczowego else (od ang. w przeciwnym przypadku) według składni (por. rys. B): 


if (warunek) instrukcjaT; 
else instrukcjaF; 


T 


ED 


Rysunek 3: Schemat blokowy instrukcji warunkowej if... else. 


Jeżeli chcemy wykonać warunkowo kilka instrukcji, tworzymy w tym celu tzw. blok in- 
strukcji, ograniczony nawiasami klamrowymi /... |. Aby zwiększyć czytelność kodu, instruk- 
cje w bloku powinniśmy wyróżnić wcięciem. Jest to jedna z zasad dobrego pisania programów. 


Rozważmy przykład, w którym wyznaczane jest minimum z dwóch liczb całkowitych. 
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int x, y; 
cin >> x >> y; // wprowadź x i y z klawiatury 


if (x < y) // tutaj nie ma średnika! 
cout << X; 

else // tutaj nie ma średnika! 
cout << p; 


Instrukcje warunkowe można, rzecz jasna, zagnieżdżać. Oto przykład służący do znajdywa- 
nia minimum z trzech liczb. 


int x, y, z; 
cin >> x >> y >> Z; 


if (x < y) I // klamerka może się znaleźć też w osobnym wierszu 
if (z < x) 
cout << Z; 
else 
cout << x; 


) 
else 
{ 
if (z < y) 
cout << z; 
else 
cout << y; 
) 


Z zagniezdzaniem instrukcji warunkowych trzeba jednak uważać. Słowo kluczowe else do- 
tyczy najbliższej instrukcji if. Zatem poniższy kod zostanie wykonany przez komputer inaczej 
niż sugerują to wcięcia. 


int x, Y, Zi 
čin >> X >> y >> Z; 


if (x != 0) 
if (y > 0 & z > 0) 
cout << yxz/x; 
else // else dotyczy if (y > 0 śe z > 0) 
cout << ":—("; 


Aby dostosować ten kod do widocznej intencji programisty, należałoby objąć fragment 
if (y >0%6é2>0) cout << yx*z/x; 


nawiasami klamrowymi. 
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3 Pętle 


Oprócz instrukcji warunkowej if, rozgałęziającej przebieg sterowania „w dół”, możemy sko- 
rzystać z tzw. pętli. Umożliwiają one wykonywanie tej samej instrukcji bądź bloku instrukcji 
wielokrotnie, być może na innych danych, dopóki pewien warunek logiczny jest spełniony. 

Pomysł ten jest oczywiście bliski naszemu życiu. Pętle znajdują zastosowanie, jeśli zacho- 
dzi konieczność np. wykonania pewnych podobnych operacji na każdym z elementów danego 
zbioru. Np. życząc sobie zsumować wydatki poczynione na przyjemności podczas każdego dnia 
wakacji, rozpatrujemy dzień pierwszy, potem dzień drugi, a potem dzień trzeci, a potem itd., 
czyli tak naprawdę dzień i-ty, gdzie ¿ = 1,2,...,n. Dalej, a ileż to razy słyszeliśmy słowa 
mądrej mamy „dopóki (Inauczysz_się) zostajesz_w_domu;”. To jest również przykład (nieko- 
niecznie świadomie użytej) pętli. 


W języku C++ zdefiniowane zostały trzy instrukcje realizujące tego typu ideę. 


— while, 
— for oraz 
— do ... while. 


Przyjrzymy się im dokładnie w poniższych podrozdziałach. 


3.1 Pętla while 


Najprostszą konstrukcją realizującą pętlę jest instrukcja while. Składnia: 


while (warunek) instrukcja; 


gdzie warunek jest pewnym wyrażeniem sprowadzalnym do typu bool. Najlepiej będzie, gdy 
istnieje możliwość jego zmiany na skutek wykonywania podanej instrukcji . W przeciwnym 
wypadku stworzymy program, który nigdy się nie zatrzyma. 

Schemat blokowy przedstawionej pętli zobrazowany jest na rys. 


T 


false 


(errn) 


Rysunek 4: Schemat blokowy pętli while. 
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Tak samo jak w przypadku instrukcji if, zamiast pojedynczej instrukcji we wszystkich oma- 
wianych pętlach może pojawić się cały ich blok, ujęty w nawiasy klamrowe. 


Uwaga 


Sam średnik (;) jest instrukcją pustą, dlatego pętla 


¡while (true) ; 


będzie w nieskończoność nie robić nic. Nie bierzmy z niej przykładu. 


Oto sposób na wypisanie na ekran kolejno liczb 1,2,...,100 i zsumowanie ich wartości 
(przypomnijmy sobie w tym miejscu problem młodego Gaussa z pierwszego wykładu). 


¡int i = 1; 
int suma = 0; 


N 


« While (i<=100) // ten warunek nie będzie kiedyś spełniony... 


5| { 

6 cout << i << endl, 

7 suma += i; 

8 i++; /* ... gdyż jest zależny od wartości zmiennej i, 
9 która się w tym miejscu zmienia */ 

10) 


p|cout << "Suma=" << suma << endl; 


3.2 Pętla for 


W pewnych zastosowaniach mamy czasem do czynienia z następującym schematem. 


instrukcjaP; // instrukcja inicjująca 
while (warunek) 
( 
instrukcjaX; 
instrukcjaA; // instrukcja aktualizująca 


Zauważmy, że przebiega on według następującego schematu. Najpierw przygotuj się przed 
wykonaniem pewnej pracy, wykonując inicjującą instrukcjęP . Potem, dopóki warunek pozo- 
staje spełniony, wykonuj instrukcjęX , jednakże aktualizując za pomocą instrukcjiA pewne 
dane przed każdą kolejną iteracją. 


Taką konstrukcję można zapisać równoważnie, korzystając z pętli for, której schemat blo- 
kowy przedstawia rys.|5| A oto jej składnia. 


for (instrukcjaP; warunek; instrukcjaA) 
instrukcjaX; 
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| instrukcjaP 


true false 


warunek 


instrukcjaX 


instrukcjaA 


(ici) 


Rysunek 5: Schemat blokowy petli for. 


Wróćmy na chwilę do przykładu omówionego w poprzednim paragrafie. Można go także 
rozwiązać, korzystając z wprowadzonej właśnie pętli. 


int suma = 0; 


for (int i=l; i<=100; ++i) 

{ 
cout << i << endl, 
suma += i; 


1 


cout << "Suma=" << suma << endl; 


Jeśli chcemy umieścić więcej niż jedną instrukcję inicjującą bądź aktualizującą, należy 
każdą z nich oddzielić przecinkiem (nie średnikiem, ani nie tworzyć dlań bloku). Zatem ten 
kod jest równoważny poniższemu, chociaż traci on znacznie na czytelności. 


int suma = 0; /* chcemy, by ta zmienna była dostępna również 


po zakończeniu działania pętli for... */ 


for (int i=l; i<=100; suma += i, ++i) // albo nawet suma += i++ 
cout << i << endl; 


1 


cout << "Suma=" << suma << endl; 


3.3 Pętla do... while 


Inną odmianą pętli jest konstrukcja do ... while, określona według składni: 
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do 
instrukcja ; 


while (warunek); 


Jak widać na schemacie blokowym (rys.[6), uzywamy jej, gdy chcemy zapewnié, ze instrukcja 
zostanie wykonana co najmniej jeden raz. 


instrukcja 


Rysunek 6: Schemat blokowy pętli do... while. 


A oto kolejna wersja rozpatrywanego przez nas przykładu. 


int i = 1, suma = 0; 
2| do 
al 
suma += 1; 
cout << i << endl; 
++i; 
} 
while (i<=100); 


3.4 break i continue 


W pewnych (dość rzadkich) szczególnych przypadkach zachodzi konieczność natychmiastowej 
zmiany przebiegu wykonywania pętli. 

Instrukcja break służy do natychmiastowego wyjścia z pętli (niezależnie od wartości wa- 
runku testowego). 

Z kolei instrukcja continue służy do przejścia do kolejnej iteracji pętli (ignorowane są in- 
strukcje następujące po continue). W przypadku pętli for instrukcja aktualizująca jest jednak 
wykonywana. 

Schemat blokowy z rys.|/|przedstawia zmianę przebiegu programu za pomocą omawianych 
instrukcji na przykładzie pętli for. 
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break _--:-)-. continue 


Rysunek 7: Instrukcje break i continue a pętla for. 


W przypadku zastosowania kilku pętli zagnieżdżonych, instrukcje te dotyczą tylko jednej 
(wewnętrznej) pętli. 
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4 Cwiczenia 
Zadanie 4.1. Wyraź następujące pętle dane w sposób opisowy za pomocą instrukcji for. 


a) dlaż =0,1,...,n — 1 wypisz 2 (dla pewnego n € N), 

b) dla ¿ = n,n — 1,...,O wypisz 2 (dla pewnego n € N), 

c) dla j = 1,3,...,2k — 1 wypisz j (dla pewnego k € N), 

d) dla i =1,2,4,7,...,n wypisz 1 (dla pewnego n € N), 

e) dla j =1,2,4,8,16,...n wypisz j (dla pewnego n € N), 

f) dla j = 1, 2,4,8,16,..., 2% wypisz j (dla pewnego k € N), 

g) dla z =a,a+0,a + 20,...,b wypisz x (dla pewnych a,b,ó € R,a <b,ó > 0). 


Zadanie 4.2. Zaprogramuj w języku C++ algorytm Euklidesa do wyznaczania największego 
wspólnego dzielnika dwóch liczb (zob. zestaw zadań nr 1). 


Zadanie 4.3. Spośród liczb 1,2,..., 100 wypisz na ekran wszystkie podzielne przez 7, tzn. 
7,14,21,.... 


Zadanie 4.4. Spośród liczb 1,2,..., 100 wypisz na ekran wszystkie podzielne przez 2 lecz 
niepodzielne przez 5, tzn. 2,4,6,8,12,.... 


Zadanie 4.5. Spośród liczb 1,2,..., 100 wypisz na ekran co drugą podzielną przez 5 lub po- 
dzielną przez 7, tzn. 5, 10, 15, 21, 28,.... 


x Zadanie 4.6. Napisz fragment kodu, który sprawdzi, czy następujące liczby są pierwsze: 
7,93,97,6687,6689, 6691. 


Zadanie 4.7. Napisz fragment kodu, który znajduje minimum z danych liczb catkowitych do- 
datnich. Liczby odczytuj z klawiatury, póki użytkownik nie wprowadzi 0. 


Zadanie 4.8. Napisz fragment kodu, który znajduje różnicę między maksimum a minimum 
z danych liczb rzeczywistych nieujemnych. Liczby odczytuj z klawiatury, póki użytkownik nie 
wprowadzi liczby ujemnej. 


Zadanie 4.9. Napisz fragment kodu, który aproksymuje wartość liczby r za pomocą wzoru 


i- 4 
sód oo AE E 
T ( 3'3 cs) 


Wypisz kolejne przybliżenia korzystając z 1,2,3,...,25 początkowych wyrazów tego szeregu. 


Zadanie 4.10. Napisz fragment kodu, który aproksymuje wartość liczby e za pomocą wzoru 


ei Ad a 
1! | 


gdzie n! = 1 x 2x --- x n. Wypisz wynik dopiero wtedy, gdy różnica pomiędzy kolejnymi 


wyrazami szeregu będzie mniejsza niż 107”. 


Zadanie 4.11. Napisz fragmenty kodu, które posłużą do wyznaczenia wartości następujących 
wyrażeń. 
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a) 2” dla pewnego n € N, 
OBIE 
5 i 
d) IE- 141? a 
e) e” = 5% $ dla pewnego x € R, 
Ð n(1+2) 2 010, EL gn dla pewnego z € [-1, 1), 


n=l n 
g) sinz ~ y an dla pewnego z € R, 
h) cosx œ 55109 A x?" dla pewnego z € R, 


100 (2n)! 


sea zana? dla pewnego z € (—1, 1). 


1) aresina Y 
Przypomnijmy, że zgodnie z umową np. 0,4 = 1+2+ -:: + 10 oraz np. I$; zal = 


Loa 
a XA 


| 


Zadanie 4.12. Korzystając ze wzoru na przybliżoną wartość funkcji sin podanego w zad. 
utwórz kod, który wydrukuje tablice przybliżonych wartości sin x dla x = ET, K=U, Lom 
i pewnego n, np. n = 10. Wynik niech będzie postaci podobnej do poniższej. 


x sin(x) 
0.0000000 0.0000000 
0.3141593  0.3090170 


3.1415927 0.0000000 


Zadanie 4.13. Pobierz wartości zmiennych a,b typu double z klawiatury. Będą one definio- 
wać równanie względem niewiadomej x € R postaci a x + b = 0. Zaproponuj fragment kodu 
w języku C++, który wyznaczy jego rozwiązanie. Poprawnie identyfikuj przypadek, gdy dane 
równanie nie jest równaniem liniowym. 


Zadanie 4.14. Dla danych a,b, c,d, e, f € R zaproponuj kod w języku C++ do rozwiązywania 
układu dwóch równań liniowych względem niewiadomych x,y € R: 


ar+by = c, 
dr +ey = f. 


Algorytm ten powinien poprawnie identyfikować przypadki (np. wypisując stosowny komuni- 
kat na ekranie), w których dany układ jest sprzeczny bądź nieoznaczony. Współczynniki układu 
pobierz z klawiatury. Do reprezentacji zbioru R użyj typu double. 


x Zadanie 4.15. Napisz fragmenty kodu, które sprawdzą, czy następujące zdania logiczne są 
tautologiami. 


a) pA —p, d) =(pA q) € =p V q, 
b) =(=p) = p, e) pv(qAr) $(pVqg)A(pVr), 
©) pvVq *pVG, ) (pq) e (p Vq) A^ (pV q). 


x Zadanie 4.16. Dla danej zmiennej typu int napisz program, który wypisze na ekran jej war- 
tość w postaci binarnej. 
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5 Wskazówki do ćwiczeń 
Wskazówka do zadania|4.6] Dla liczby k co najwyżej należy sprawdzić, czy dzieli sig ona bez 
reszty przez każdą z liczb 2,3,..., vk. 


Wskazówka do zadania Użyj zagnieżdżonych pętli for do sprawdzenia wszystkich 
możliwych podstawień wartości logicznych pod zmienne p, q, r: 


for (int p=0; p<=l; ++p) // bool(0) to false , bool(1) to true 
for (int q=0; q<=1; ++q) 
for (int r=0; r<=l; ++r) 


// sprawdzenie warunku... 


Tutaj może się przydać (ale nie musi) instrukcja break. 


Wskazówka do zadania 4.16] Liczba typu int jest 32 bitowa. Należy użyć pętli, która będzie 
w każdej iteracji wypisywać wartość kolejnego bitu (0 lub 1). Wyrażenie x & (1<<k) przyj- 
muje wartość 0, gdy k-ty bit jest równy 0 oraz wartość różną od 0, gdy k-ty bit jest równy 1 
(dlaczego ?). 
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1 Tablice jednowymiarowe 


Do tej pory przechowywaliśmy dane używając pojedynczych zmiennych. Były to tzw. zmienne 
skalarne (atomowe). Jedna jednostka danych odpowiadała pojedynczej zmiennej. 

Często jednak mamy dany cały ciąg zmiennych tego samego typu. Dla przykładu, rozważmy 
fragment programu dokonujący podsumowania rocznych zarobków pewnej doktorantki. 


double zarobkil , zarobki? , /*... */ , zarobkil2; 
// deklaracja 12 zmiennych 
zarobkil = 1399.0; // styczeń 
zarobki2 = 1493.0; // luty 
ELO Gna 
zarobkil2 = 999.99; // grudzień 


double suma = 0.0; 


suma += zarobkil; 
suma += zarobki2; 
//... 


suma += zarobkil2; 


" 


cout << "Zarobiłam w 2010 r. " << suma << Zła "2 


N 


A 


Dochód z każdego miesiąca przechowywany jest w oddzielnej zmiennej. Wszystkie z nich 
są tego samego typu. Łatwo zauważyć, że nie jest to zbyt wygodne. 


W tym wykładzie zajmiemy się tablicami jednowymiarowymi o ustalonym rozmiarze. 
Tablice jednowymiarowe (ang. one-dimensional arrays) są reprezentacją znanych nam obiek- 
tów matematycznych: ciągów skończonych bądź wektorów. 

Deklaracja tablicy jednowymiarowej następuje według poniższej składni: 


typ identyfikator |rozmiar |; 


gdzie rozmiar € N, rozmiar>0, jest stałą (!). 


Dostęp do elementów tablicy odbywa się za pomocą operatora indeksowania [-]. Elementy 
są numerowane od 0 do n — 1, gdzie n to rozmiar tablicy. Na przykład: 


int t[5]; // deklaracja tablicy (5 elementów typu int) 


// t[0] to pierwszy element (!) 
// t[4] to ostatni element (!) 


Operator indeksowania przyjmuje jako parametr dowolną wartość całkowitą (np. stałą bądź 
wyrażenie arytmetyczne). Każdy element tablicy traktujemy tak, jakby był zwykłą zmienną, 
z którą do tej pory mieliśmy do czynienia. 


Uwaga 


W języku C++ nie ma mechanizmów sprawdzania poprawności indeksów! Następujący 
kod najczęściej nie spowoduje natychmiastowego wystąpienia błędu. 
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int t[5]; 
t[—100] LTS: 4% s= 
3 1[10000] = 25326; // :-( 


N 


Jednakże, zmieni on wartość komórek pamięci wykorzystywanych przez inne funkcje. Skutki 
tego działania mogą się objawić w innym miejscu programu, powodując nieprzewidywalne 
i trudne do wykrycia błędy. 

Więcej na temat tego zagadnienia w wykładzie dotyczącym dynamicznego przydziału pa- 
mięci. 


Wróćmy do przykładu dotyczącego naszej zaradnej koleżanki. Oto w jaki sposób można 
wykorzystać tablice do podsumowania jej zarobków. Zauważmy analogię pomiędzy tym a po- 
przednim fragmentem kodu. 


¡double zarobki[12];  // deklaracja 12 elementowej tablicy 


s |zarobki[0] 1399.0; // styczeń 
a|zarobki[1] = 1493.0; // luty 

5/7 ... 

s zarobki[11] = 999.99; // grudzień 


s|double suma = 0.0; 
o for (int i=0; i<=1l; ++i) // rozważ wszystkie elementy... 
10 suma += zarobki[i|; 


" " 


<< suma << 


e cout << "Zarobitam w 2010 r. zł."; 


Przetworzenie całej tablicy za pomocą pętli for i dokonanie na każdym elemencie tej samej 
operacji arytmetycznej sprawia, że kod staje się bardziej zwięzły i czytelny. Z drugiej strony, 
nietrudno go teraz byłoby poprawić tak, żeby np. uwzględniał zarobki z całego okresu jej stu- 
diów. 


Ponadto, dla wygody, przy tworzeniu tablicy można przypisać wartości jej elementom w na- 
stępujący sposób: 


typ identyfikator[n] = { 


wartość0 , wartośćl , /* ... */, wartośćn 
J; 
Pozwala to jeszcze bardziej uprościć program dla naszej doktorantki: 
ıı const int n = 12; /x stała! »/ 
2| double zarobki|n| = | 1399.0, 1493.0, /*...*/, 999.99 |; 


double suma = 0.0; 


a 
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00 
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for (int i=0; i<n; ++i) 
suma += zarobkili]; 


" " 


cout << "Zarobiłam w 2010 r. << SUMA << ZŁ 


O tablicach dowolnego rozmiaru (nie definiowanego z góry na etapie pisania kodu) oraz 
sposobach przekazywania tablic innym funkcjom dowiemy się więcej z wykładu na temat dy- 
namicznej alokacji pamięci. 


Na marginesie, aby wypisać zawartość całej tablicy, można użyć następującego fragmentu 
kodu: 


const int n = ...; 
int sablal = | «-+ Y 


for (int i=0; i<n; ++i) 
cout << t[i] << endl; 


Z drugiej strony, aby wczytać elementy tablicy z klawiatury, można napisać: 


const int n = ...; 
int tabln]; 


" 1 


cout << "Podaj wartości " << n << elementów." << endl; 
for (int i=0; i<n; ++i) 


cin >> t[i]; 


2 Sortowanie tablic 


Rozpatrzmy problem sortowania tablic jednowymiarowych. Jest to jedno z ciekawszych za- 
gadnień, które pojawiło się na początku rozwoju informatyki teoretycznej, a w którym zastoso- 
wanie znajdują właśnie tablice. Można je sformalizować w sposób następujący. 

Dana jest tablica £ rozmiaru n zawierająca elementy, które można porównywać za pomocą 
operatora relacyjnego <=. Należy zmienić kolejność (tj. dokonać permutacji) elementów 1 tak, 
by zachodziło 


t[0] <=+[1], +[1] =t[2], ..., t[n—2] <=t[n—l]. 


Oczywiście rozwiązanie takiego problemu nie musi być jednoznaczne. Dla tablic zawiera- 
jących elementy t[i] i t[j] takie, że dla ¿Aj zachodzi t[i] <=t[j] oraz t[j] <= t[i], może 
istnieć wiele rozwiązań spełniających powyższy warunek. 


W niniejszym paragrafie omówimy trzy elementarne algorytmy: 


a) sortowanie przez wybór, 
b) sortowanie przez wstawianie, 
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c) sortowanie bąbelkowe. 


Cechują się one tym, że w pesymistycznym (,,najgorszym”) przypadku liczba operacji porów- 
nań elementów tablicy jest proporcjonalna do n°. Inną ważną cechą, którą poddamy analizie, 
będzie liczba potrzebnych przestawień tych elementów. 

Bardziej wydajne (i bardziej złożone, np. sortowanie szybkie, przez łączenie, przez kop- 
cowanie) algorytmy będą omówione w semestrze III. Niektóre z nich wymagają co najwyżej 
kn log, porównań dla pewnego k. 


2.1 Sortowanie przez wybór 


W algorytmie sortowania przez wybór (ang. selection sort) dokonujemy wyboru elementu naj- 
mniejszego spośród jeszcze nieposortowanych, póki cała tablica nie zostanie uporządkowana. 
Ideę tę formalizuje następujący pseudokod: 


dla i=0,1,...,n—2 
( 
// t[O], ..., t[i—-1] są uporządkowane względem <= 
// (nadto, są juź na swoich ostatecznych miejscach) 
j = indeks najmniejszego elementu 
spośród t[il, ..., t[n—l]; 
zamień elementy t[i] i t[j]; 
} 


Jako przykład rozważmy, krok po kroku, przebieg sortowania tablicy int t[5] = {4,1,3,5,2}; 
Kolejne iteracje działania tego algorytmu ilustrują rys. 

W kroku I (rys.|1) mamy i==0. Dokonując wyboru elementu najmniejszego spośród £[0], 
..., L[4] otrzymujemy ¡==1. Zamieniamy więc elementy z[0] i +[1] miejscami. 

W kroku II (rys. |2) mamy i==1. Wybór najmniejszego elementu wśród t[1],..., t[4] daje 
j==4 Zamieniamy miejscami zatem £[1] it[4]. 

Dalej (rys. B), i==2. Elementem najmniejszym spośród t[2],..., t[4] jest +[7] dla j==2 
Zamieniamy miejscami zatem niezbyt sensownie £ [2] i t[2]. Komputer, na szczęście, zrobi to 
bez grymasu. 

W ostatnim kroku (rys. |4) i==3 i j==4, dzięki czemu możemy uzyskać ostateczne rozwią- 
zanie (rys. [5). 


Implementację sortowania przez wybór w języku C++ i analizę liczby porównań oraz prze- 
stawień elementów pozostawiamy jako ćwiczenie. 
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Rysunek 1: Sortowanie przez wybór — przykład — iteracja I. 
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Rysunek 2: Sortowanie przez wybór — przykład — iteracja II. 
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Rysunek 3: Sortowanie przez wybór — przykład — iteracja III. 
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Rysunek 4: Sortowanie przez wybór — przyktad — iteracja IV. 


Rysunek 5: Sortowanie przez wybór — przykład — rozwiązanie. 
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2.2 Sortowanie przez wstawianie 


Algorytm sortowania przez wstawianie (ang. insertion sort) jest metodą często stosowaną 
w praktyce do porządkowania małej liczby elementów (do ok. 20-30) ze względu na swą pro- 
stotę i szybkość działania. 

W niniejszej metodzie w i-tym kroku elementy £ [0], ..., t[i—1] są już wstępnie upo- 
rządkowane względem relacji <=. Pomiędzy nie wstawiamy £[i] tak, by nie zaburzyć porządku. 


Formalnie rzecz ujmując, idea ta może być wyrażona za pomocą pseudokodu: 


dla ¡=1,2,...,n-1 


( 
// t[O], ..., t[i—-1] są uporządkowane względem <= 
// (ale niekoniecznie jest to ich ostateczne miejsce) 
j = indeks takiego elementu spośród t[0],...,t[i], że 
t[k] <= t[i] dla każdego k < j oraz 
t[l] > t[i] dla każdego I > j; 
wstaw t[i] przed t[j|; 
) 


gdzie przez operację „wstaw t[i] przed t[j]”,dla0 < j < i rozumiemy ciąg działań, mający 
na celu przestawienie kolejności elementów tablicy: 


t [0] = (6-1) | tl] e Ei | tli] | t[i+1] = t[n—1] 
tak, by uzyskać: 
t [0] = ei | tli] t[j] e t[i—1] | t[i+1] en t[n—1] 


Powyższy pseudokod może być wyrażony w następującej równoważnej formie: 


dla ¡=1,2,...,n-1 


( 
// t[O], ..., t[i—-1] są uporządkowane względem <= 
// (ale niekoniecznie jest to ich ostateczne miejsce) 
znajdź największe j ze zbioru (0,...,i) takie, że 
j == 0 lub t[j—l] <= t[i]; 
wstaw t[i] przed t[j|; 
} 


Jako przykład rozpatrzmy znów tablicę int 1[5] = {4,1,3,5,2}; 
Przebieg kolejnych wykonywanych kroków przedstawiają rys. 


Implementację sortowania przez wstawianie i analizę liczby porównań oraz przestawień 
elementów pozostawiamy jako ćwiczenie. 
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Rysunek 6: Sortowanie przez wstawianie — przykład — iteracja I. 
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Rysunek 7: Sortowanie przez wstawianie — przykład — iteracja II. 
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Rysunek 8: Sortowanie przez wstawianie — przykład — iteracja III. 
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Rysunek 9: Sortowanie przez wstawianie — przykład — iteracja IV. 
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2.3 Sortowanie bąbelkowe 


Sortowanie bąbelkowe (ang. bubble sort) jest interesującym przykładem algorytmu pojawia- 
jącego się w większości podręczników akademickich dotyczących podstawowych sposobów 
sortowania tablic, który jednakże prawie wcale nie jest stosowany w praktyce. Jego wydajność 
jest bowiem bardzo słaba w porównaniu do dwóch metod opisanych powyżej. Z drugiej strony, 
posiada on sympatyczną „hydrologiczną” (nautyczną?) interpretację, która urzeka wielu wykła- 
dowców, w tym i skromnego autora niniejszego opracowania. Tak umotywowani, przystąpmy 
więc do zapoznania się z nim. 

W tym algorytmie porównywane są elementy tylko ze sobą bezpośrednio sąsiadujące. Jeśli 
okaże się, że nie zachowują one odpowiedniej kolejności względem relacji <=, element ,,cigz- 
szy” wypychany jest „w górę”, niczym pęcherzyk powietrza (tytułowy bąbelek) pod powierzch- 
nią wody. 


A oto pseudokod: 
dla i=n—1 ,...,l 
( 
dla ¡=0,...,¿-1 
( 
// porównuj elementy parami 
jeśli (t[j] > t[j+1]) 
zamień t[j] i t[j+l]; 
// tzn. "wypchnij" cięższego "'bąbelka" w górę 
} 
// tutaj elementy t[i],...,t[n—1] są już na swoich miejscach 
} 


Dla przykładu rozpatrzmy ponownie tablicę int 1[5] = (4,1,3,5,2); 
Przebieg kolejnych wykonywanych kroków przedstawiają rys. 


Implementację sortowania bąbelkowego w języku C++ i analizę liczby porównań oraz prze- 
stawień elementów pozostawiamy jako ćwiczenie. 
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Rysunek 10: Sortowanie bąbelkowe — przykład — krok I. 
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Rysunek 11: Sortowanie bąbelkowe — przykład — krok II. 
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Rysunek 12: Sortowanie bąbelkowe — przykład — krok III. 


© e 


AiPP-V, s. 18. 


DE 


bo 


Rysunek 13: Sortowanie bąbelkowe — przykład — krok IV. 
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Rysunek 14: Sortowanie bąbelkowe — przykład — rozwiązanie. 
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3 Cwiczenia 


Zadanie 5.1. Niech dana będzie tablica zadeklarowana jako int tab[n], dla pewnego n. Napisz 
kod, który przesunie każdy element o indeksie > 0 o jedną komórkę w lewo, a element pierwszy 
wstawi na miejsce ostatniego. 


Zadanie 5.2. Za pomocą tylko jednej pętli for znajdź w tablicy double tab[n] element naj- 
mniejszy i największy. 


Zadanie 5.3. Napisz fragment kodu, który w tablicy bool tab[n] zliczy, ile razy występuje 
wartość true oraz dokona negacji wszystkich elementów. 


Zadanie 5.4. Napisz kod, który znajduje najmniejszy element w tablicy int tab[n]. Następnie 
wypełnij tym elementem wszystkie komórki o parzystych indeksach oraz elementem przeciw- 
nym do niego komórki o indeksach nieparzystych. 


Zadanie 5.5. Dla danej tablicy liczb rzeczywistych t rozmiaru n napisz kod, który wyznaczy 


wartość średniej arytmetycznej wszystkich elementów, danej wzorem = > 
Zadanie 5.6. Niech dany będzie n-elementowy ciąg liczb rzeczywistych a = (a,,a2,..., an). 
Średnią geometryczną nazywamy wartość 


Napisz program, który wczyta do tablicy n = 8 wartości z klawiatury oraz następnie policzy 
wartość ich średniej geometrycznej. 


Zadanie 5.7. Niech dany będzie n-elementowy ciąg liczb rzeczywistych a = (a,,02,...,dn). 
Średnią ważoną względem wektora wag w = (wy, W2, . . . , Wn) o elementach nieujemnych oraz 
takiego, że );_, w, = 1, nazywamy wartość 


WMwy(a) = 5 QiWi. 
i=1 


Napisz program, który dla n = 5 wczyta z klawiatury wektor wag w i sprawdzi, czy spełnia on 
postawione wyżej założenia oraz wyznaczy wartość średniej ważonej ciągu (—2, —1,0, 1, 2). 


Zadanie 5.8. Niech dany będzie n-elementowy ciąg liczb rzeczywistych a = (a,,a2,..., an). 
Operatorem maks-min względem wektora wag w = (wy, W9, ..., Wn) składającego się z war- 
tości rzeczywistych, nazywamy wartość 


MaxMinw(a) = max (min{a;, wi}). 


Jn. 


Napisz program, który dla n = 5 wczyta z klawiatury ciąg a, następnie wyznaczy wartość 
operatora maks-min względem wektora wag (1,2,...,n). 
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Zadanie 5.9. Niech dany będzie wektor o elementach rzeczywistych x = (£1, ..., £n). Napisz 
program, który wczyta wartości jego elementów z klawiatury (dla n = 9) oraz policzy wartość 
jego normy euklidesowej, wg wzoru 


Zadanie 5.10. Niech dane będą dwa n-elementowe wektory o elementach rzeczywistych x = 
(X1,...,Ty)iy = (Y1, -.-, Yn). Napisz program, który wczyta wartości ich elementów z kla- 
wiatury (dla n = 5) oraz policzy wartość ich odległości w metryce supremum, wg wzoru 


Zadanie 5.11. Niech dany będzie wielomian rzeczywisty stopnia n, w(x) = w[0]z”+w[l]x"+ 
. ..W[n]zx”, w[n] Æ 0, którego współczynniki przechowywane są w n + 1 wymiarowej tablicy 
o elementach typu double. Napisz program, który wczytuje wartość współczynników dla n = 3 
oraz wartość x i wyznacza wartość w(x) wg powyższego wzoru. 


x Zadanie 5.12. Zmodyfikuj program z zad. tak, by korzystał z bardziej efektywnego 
obliczeniowo wzoru na wartość w(x), zwanego schematem Hornera: 


w(x) = (---(((w[n]z + w[n—1]) z + w[n—2]) z + w[n—3]) ---) x + w[0]. 


x Zadanie 5.13. Niech dane będą wielomiany w(x) i v(x) stopnia, odpowiednio, n i m. Napisz 
program, który wyznaczy wartości współczynników wielomianu u(x) stopnia n + m, będącego 
iloczynem wielomianów w i v. Dokonaj obliczeń dla w(x) = 1* + 41? — x + 2 oraz v(x) = 
zf + z” + 10. 


Zadanie 5.14. Zaimplementuj algorytm sortowania przez wybór dla danej tablicy o n elemen- 
tach typu int. Oblicz, ile jest potrzebnych operacji porównań oraz przestawień elementów w za- 
leżności od n. 


Zadanie 5.15. Zaimplementuj algorytm sortowania przez wstawianie dla danej tablicy o n ele- 
mentach typu int. Oblicz, ile jest potrzebnych operacji porównań oraz przestawień elementów 
w zależności od n dla tablicy już posortowanej oraz dla tablicy posortowanej w kolejności od- 
wrotnej. 


Zadanie 5.16. Zaimplementuj algorytm sortowania bąbelkowego dla danej tablicy o n ele- 
mentach typu int. Oblicz, ile jest potrzebnych operacji porównań oraz przestawień elementów 
w zależności od n dla tablicy już posortowanej oraz dla tablicy posortowanej w kolejności od- 
wrotnej. 


x Zadanie 5.17. Dana jest tablica o n>1 elementach typu double. Napisz funkcję, która ob- 
liczy wariancję jej elementów używając tylko jednej pętli for. Wariancja elementów ciągu 
X = (11,...,T,) dana jest wzorem 


AAi l py 2 
L => 2,(—X) , (1) 


i=1 


gdzie X jest średnią arytmetyczną ciągu x. 
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4 Wskazówki do ćwiczeń 


Odpowiedź do zadania 


u(x) = £ + £7 +4e + 30% + 11% + 2a” + 40a? — 10x + 20. 


Odpowiedź do zadania |5.14 


¡¡void zamien(intś x, intá y) 
| 

3 int f = x; 

4 x = Yy; 

5 y = t; 

6| } 

R 

s int main() 

o| { 

0 // za 

1 const int n = ...; 

2 int jla] "| sss J 

2 

4 for (int i=0; i<n—l; ++i) 
5 ( 

6 // znajdź indeks najmn. el. 
7 int j= i; 

8 for (int k=i+1; k<n; ++k) 
9 if (t[k] < t[j]) 

20 j S 

21 

2 zamiení(t[il, t[j]); 

2 } 

24 // 

|) 


t[i], t[n=l]; 


spośród 
elementów 


// porównanie 


// zamiana elementów 


Analiza liczby porównań elementów. Liczba ta jest ta sama niezależnie od kolejności da- 
nych wejściowych i wynosi (n — 1) + (n — 2) +--+ 1 = n(n — 1)/2 = n? /2 — n/2. 

Liczba zamian elementów w powyższej implementacji jest również niezależna od permu- 
tacji danych wejściowych i wynosi n — 1. W prosty sposób można zmniejszyć tę liczbę tak, 
by w przypadku ciągu już posortowanego nie wykonywać żadnych przestawień. 


Odpowiedź do zadania |5.15 
Przykładowa implementacja: 


¡int main() 


2 
3 // . 
4 const int n = ...; 
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5 int t[n] = £ ... ); 

6 

7 for (int i=l; i<n; ++i) 

3 ( 

9 // znajdź największe j ze zbioru [0,...,i)] takie, że 
0 // j == 0 lub t[i] > t[j-1]; 

1 int j; 

2 for (j=i; j>0; ==) 

3 if (t[j—1] <= t[i]) // porównanie 
4 break; 

5 

6 // wstaw t[i] przed t[j]; 

7 for (int k = i—l; k>=j; k—-) 

8 zamien(t[k|, t[k+1]); // zamiana 
9 ) 

20 // 

a|} 


Dla danych już posortowanych liczba porównań wykonywanych przez powyższy kod wy- 
nosi n — 1, a liczba przestawień 0. 

Dla danych odwrotnie posortowanych liczba przestawień i porównań równa jest 1 + 2 + 
-+ (n — 1) = n(n — 1)/2. 


Odpowiedź do zadania|5.16 


Jint main() 

| 

3 // zi 

4 const int n = ...; 

5 int slal = era $ 

6 

7 for (int i=n—-l; i>0; —-1) 

3 ( 

9 for (int j=0; j<i; ++) 

o ( 

1 // porównuj elementy parami 

2 if (t[j] > t[j+1]) // porownanie 
3 zamien(t[j]|, t[j+1]); // zamiana 
4 } 

5 ) 

6 

7 // 

s|} 


Liczba porównań jest niezależna od permutacji ciągu wejściowego i wynosi (n — 1) + (n — 
2) +- +1 = n(n — 1)/2. 
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Liczba przestawień dla danych posortowanych wynosi 0. Dla danych odwrotnie posortowa- 
nych równa jest n(n — 1)/2. 


Wskazówka do zadania Należy znaleźć wzory (tzw. rekurencyjne) na średnią arytme- 
tyczną i wariancję k początkowych elementów ciągu, zależne od elementu x, oraz średniej 
arytmetycznej i wariancji elementów (11,...,Tx-1). 
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1 Funkcje 


1.1 main() 


Najprostszy pełny kod źródłowy działającego programu w języku C++ można zapisać w sposób 
następujący. 


¡int main() 


( 


3 return 0; // zakończenie programu (,,sukces `’) 


N 


Każdy program w języku C++ musi zawierać definicję obiektu o nazwie main(), inaczej jego 
kompilacja nie zakończy się powodzeniem. main() (od ang. główny) stanowi punkt startowy 
każdego programu. Instrukcje w nim zawarte zostają wykonane jako pierwsze. 

De facto, obiekt ten jest funkcją. 

W niniejszej części materiałów przyjrzymy się sposobom deklaracji i użycia funkcji, które 
stosuje się, aby uprościć podprogram główny i zwiększyć czytelność kodu. 


Uwaga 


Znakomita większość pisanych przez nas programów będzie jednak wpisywać się w poniż- 
szy schemat. 


|#include <iostream> // Wczytanie biblioteki 
2 // systemowej iostream 
3 using namespace std; // Udostepnienie zawartych 


4 // w niej narzędzi 


s int main() 
1| { 
8 // kod 
9 // 


11 return 0; // zakończenie programu (,, sukces ””) 


2|} 


Zwróćmy uwagę na dwie pierwsze instrukcje. Służą one do wczytania i udostępnienia biblio- 
teki iostream, w której zawarte są m.in. definicje obiektów cout i cin. Dzięki nim możemy 
wypisywać komunikaty na ekranie oraz pobierać dane z klawiatury. 


1.2 Definiowanie funkcji 


Na etapie projektowania programu często można wyróżnić wiele modułów (części, podpro- 
gramów), z których każda jest odpowiedzialna za pewną logicznie wyodrębnioną, niezależną 
czynność. Fragmenty mogą cechować się dowolną złożonością. Mogą same np. składać się 
z kolejnych podmodułów. 
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Zadaniem programisty jest powiązanie tych fragmentów w spójną całość, m.in. poprzez 
zorganizowanie odpowiedniego przepływu sterowania i wymiany informacji (proces taki może 
przypominać budowanie domku z klocków różnych kształtów). 


Rozważmy dwa przykładowe zalążki programów, obrazujące pewne sytuacje z tzw. życia. 
Najpierw przyjrzyjmy się czynnościom potrzebnym do uruchomienia samochodu. 


¡(int main() 
| { 
3 ustawFotel(); 

4 ustawLusterka(); 

; zapnijPasy(); 

6 przekręćKluczyk(); 


1 sprawdźKontrolki(): 
8 włączŚwiatła(): 
9 return 0; 


10| ) 


Programiście, który wyróżnił ciąg czynności potrzebnych do wykonania pewnego zadania po- 
zostaje tylko szczegółowe określenie (implementacja), na czym one polegają. Zauważmy, że 
dzięki takiemu sformułowaniu rozwiązania możliwe jest łatwiejsze zapanowanie nad złożono- 
ścią kodu. 

Drugi przykład dotyczy organizacji nauki w semestrze pewnego pilnego studenta. 


¡int main () 

| 

3 EF aws 

4 do { 

5 if (dzieńTygodnia == niedziela) 
6 continue; 

; 

8 pouczSięTrochę(),; 

9 

0 if (zmęczony) { 

1 if (godzinNaukiDziś > 5) 

2 mozeszWreszcielséNaRandke () ; 
3 else 

4 odpocznijChwilę(); 

5 ) 

6 } while (!nauczony); 

7 // 

s|) 


Najprostszym praktycznym sposobem podziału programu na moduły jest użycie funkcji] 


IW semestrze II poznamy jeszcze inny sposób, pozwalający na pisanie naprawdę dużych programów: tzw. 
programowanie obiektowe. 
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Funkcja to odpowiednio wydzielony fragment kodu, wykonujący pewne instrukcje, do któ- 
rego można się odwoływać z innych miejsc programu. 


Projektując funkcję należy wziąć pod uwagę, jak ma ona wchodzić w interakcję z innymi 
funkcjami (np. main()), to znaczy jakiego typu dane wejściowe powinna ona przyjąć i jakiego 
typu dane wyjściowe będą wynikiem jej działania. 

Bezpośrednim skutkiem działania funkcji jest uzyskanie jakiejś konkretnej wartości na 
wyjściu (może to też być ,,nic”, które jest reprezentowane przez typ void). Skutkiem pośred- 
nim działania funkcji może być zmiana kontekstu, tj. stanu komputera (np. wypisanie czegoś 
na ekran). Idea ta została zobrazowana na rys. 


Rysunek 1: Schemat przepływu danych w funkcjach. 


FUNKCJA 


KONTEKST 


Składnia definicji funkcji: 


typZwracany nazwaFunkcji(listaArgumentów) // (nagłówek) 


( 
// kod funkcji (ciało) 


) 


gdzie: 


— typZwracany to przeciwdziedzina, np. int, bool, void itp., 

— nazwaFunkcji jest identyfikatorem (zob. wykład o deklaracji zmiennych), 

— a listaArgumentów deklaruje dziedzinę, tj. parametry wejściowe (tak jak deklaruje się 
zmienne), oddzielając je przecinkami. 


Do zwracania wartości przez funkcję (należących do określonej przeciwdziedziny) służy 
instrukcja return. Działa ona podobnie jak break w przypadku pętli, tj. natychmiast przerywa 
wykonywanie funkcji, w której aktualnie znajduje się sterowanie. Dzięki temu następuje powrót 
do miejsca, w którym nastąpiło wywołanie danej funkcji. 

Jak już wspomnieliśmy, za typZwracany można przyjąć też void, czyli „nic” (0). Jedynym 
skutkiem działania takiej funkcji jest zmiana kontekstu. W tym przypadku instrukcja return 
może zostać pominięta. Po wykonaniu ostatniej instrukcji następuje automatyczny powrót do 
funkcji wywołującej. 


1.3 Definicja a deklaracja 


Rozpatrzmy definicję funkcji służącej do wyznaczania kwadratu liczby rzeczywistej. Wyraźmy 
ją za pomocą notacji matematycznej. 
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N 


Niech 


kwadrat : R > R 
nazwa funkcji dziedzina przeciwdziedzina 
taka, że 
kwadrat(x) := z x z. 
Pierwsza część — „Niech... — powyższej definicji nazywa się deklaracją (prototypem). 


Za jego pomocą orzekamy, że mamy zamiar określić funkcję o danej nazwie, posiadającą pewną 
dziedzinę i pewną przeciwdziedzinę. W tym momencie nie wiadomo jednak jeszcze, co i jak 
dokładnie taka funkcja ma robić. 

Dlatego dalej — „taka, że ...” — następuje właściwa część definicji, która podaje do- 
kładny sposób (algorytm) przetworzenia danych wejściowych celem uzyskania spodziewanego 
wyniku. 

Powyższa definicja określa to postępowanie w następujący sposób. 

a) Weź argument wejściowy (będący liczbą rzeczywistą). Niech w tym kontekście będzie on 

identyfikowany jako x (równie dobrze mógłby to być każdy inny symbol, np. y, a, Q1). 
b) Wyznacz wartość iloczynu z X z. 
c) Rezultat obliczeń z p. b) (liczbę rzeczywistą) zwróć jako wynik. 


Przyjrzyjmy się implementacji rozpatrywanej funkcji i jej przykładowemu zastosowaniu. 


#include <iostream> 
using namespace std; 


// ,,Niech... 
double kwadrat(double x) 
{ // ,,taka, że..." 
return xxx; 
) 
int main() // funkcja bezargumentowa 


( 
double y = 0.5; 
cout << kwadrat( y) << endl; // wywotanie funkcji 
cout << kwadrat(2.0) << endl; 
return 0; 


Jako ze program czytany jest (przez nas jak i przez kompilator) od góry do dołu, definicja 
każdej funkcji musi się pojawić przed jej pierwszym użyciem. 

Jednakże czasem dobrze jest (np. dla czytelności) umieścić ją w innym miejscu (poniżej, 
w innymi pliku itp.). Można to zrobić, pod warunkiem, że obiektom z niej korzystającym udo- 
stępnimy jej deklarację. Intuicyjnie, innym obiektom wystarczy przekazać informację, że dana 
funkcja istnieje. Dokładna specyfikacja jej działania nie jest im de facto potrzebna. 

Deklaracji funkcji dokonujemy przez podanie jej nagłówka zakończonego $rednikiem w miej- 
scu, w którym nie została ona jeszcze użyta (ale poza inną funkcją), według składni: 
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typZwracany nazwaFunkcji(listaArgumentów); 


Co ważne, w deklaracji musimy użyć takiej samej nazwy funkcji oraz typów argumentów wej- 
ściowych i wyjściowych jak w definicji (która ostatecznie musi się gdzieś pojawić). Nazwy tych 
argumentów mogą być jednak inne, a nawet mogą zostać pominięte. 


Dlatego nasz przykład również mógłby być zapisany tak: 


#include <iostream> 
2| using namespace std ,; 


double kwadrat(double); // deklaracja 


int main() 
{ 
double x = kwadrat (2.0); 
cout << x << endl, 
cout << kwadrat(x)+kwadrat(8.0) << endl; 
cout << kwadrat(kwadrat(0.5)) << endl; 
return 0; 


s double kwadrat(double x) // definicja 
( 


return xxx; 


Uwaga 


Warto mieć na uwadze, że instrukcja 


cout << kwadrat(kwadrat(0.5)); 


zostanie wykonana w sposób podobny do następującego: 


double zmiennaPomocniczal = kwadrat(0.5); 
double zmiennaPomocnicza2 = kwadrat(zmiennaPomocniczal); 
cout << zmiennaPomocnicza2; 


Wartości pośrednie są obliczanie i odkładane „na boku”. Tym samym, podanie funkcji kwadrat 
(czy dowolnej innej przyjmującej jako parametr typ double) jako argumentu „kwadrat(0.5)” 
jest tożsame z przekazaniem jej wartości 0.25. 

Jest to o tyle istotne, iż w przeciwnym przypadku „kwadrat(0.5)” musiałby być wyliczony 
dwukrotnie (mamy przecież return xxx; w definicji). Nic takiego się jednak nie dzieje. 
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1.4 Biblioteki funkcji 


Zbiór funkcji podręcznych można umieścić w osobnych plikach, tworząc tzw. bibliotekę funk- 
cji. Dzięki temu program staje się bardziej czytelny, można łatwiej zapanować nad jego zło- 
żonością, a ponadto otrzymujemy sposobność wykorzystania pewnego kodu ponownie w przy- 
szłości. Jest to ważne, gdyż raz napisany i przetestowany zestaw funkcji oszczędza nam sporo 
pracy. 


Aby zatem stworzyć bibliotekę, tworzymy najczęściej dwa dodatkowe pliki. 


— plik nagłówkowy (ang. header file), nazwabiblioteki .h — zawierający tylko de- 
klaracje funkcji, 

— plik źródłowy (ang. source file), nazwabiblioteki .cpp — zawierający tylko defi- 
nicje funkcji. 


Plik nagłówkowy należy dołączyć do wszystkich plików, które wykorzystują funkcje z danej 
biblioteki (również do pliku źródłowego biblioteki) za pomocą dyrektywy: 


Hinclude "nazwabiblioteki.h" 


Zwróćmy uwagę, że nazwa naszej biblioteki, znajdującej się wraz z innymi plikami two- 
rzonego projektu, ujęta jest w cudzysłowy, a nie w nawiasy trójkątne (które ładują biblioteki 
systemowe). 

Zauważmy, że podczas laboratoriów już korzystaliśmy z jednej biblioteki. Był to zestaw 
funkcji służący do rysowania. W funkcji main() wywoływaliśmy je celem stworzenia cieszą- 
cych oczy obrazków. Z pomocą jednej biblioteki można by było napisać bardzo dużo progra- 
mów, jeden rysujący np. domek, drugi kotka itd. 


Dla ilustracji przyjrzymy się, jak wyglądałby program składający się z biblioteki zawierają- 
cej dwie następujące funkcje dane w notacji matematycznej. 


Niech 
min: Z x Z > Z, 
max: Z x Z > Z, 
takie, że 
naaa n dlan<m, 
"7 | m dam<n, 
oraz 


je n dlan > m, 
| m dlam>n. 


Oto zawartość pliku podreczna.h. 


#pragma once // na początku każdego pliku nagłówkowego! 


int min(int i, int j); 
int max(int i, int j); 
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Uwaga 


Dyrektywa 


fpragma once 


powinna się znaleźć na początku każdego pliku nagłówkowego. Zapobiega to ponownemu do- 


4: * * 


Ta dyrektywa jest dostępna tylko w kompilatorze Visual C++. W innych narzędziach uzy- 
skanie tego efektu jest nieco bardziej skomplikowane: 


#ifndef _ PODRECZNA H 
Hdefine _ PODRECZNA _H 


// deklaracje .... 


Hendif 


Definicje funkcji z naszej biblioteki znajdą się w pliku podreczna.cpp. 


Hinclude "podreczna.h" /x dołącz plik z deklaracjami x/ 


int min(int x, int y) /* nazwy parametrów nie muszą być takie 


same jak w deklaracji x/ 


( 
if (x <= y) return x; 
else return y; 
) 
int max(int w, int z) 
( 
if (w >= z) 
return w; 
return z; 
) 


Gdy biblioteka jest gotowa, można przystąpić do napisania funkcji głównej, która będzie 
posiłkować się nimi do własnych celów. Oto zawartość pliku program. cpp. 


#include <iostream> /x dołącz bibliotekę systemową */ 
using namespace std; 


Hinclude "podreczna.h" /* dołącz własną bibliotekę */ 


int main(void) // tak też można zadeklarować 


// funkcję bezargumentową 
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int x, y; 


" 


cout << "Podaj dwie liczby. "; 
cin >> x >> y; // wczytaj z klawiatury 


cout << "Max=" << max(x,y) << endl; 
cout << "Min=" << minlx,y) << endl, 


return 0; 


2 Przegląd funkcji w bibliotekach systemowych 


Elementem standardu języka C++ jest też wiele przydatnych bibliotek systemowych zawierają- 
cych często wykorzystywane funkcje. Poniżej omówimy te, które interesować nas będą najbar- 
dziej. 

2.1 Funkcje matematyczne 

Biblioteka <cmath> udostępnia wybrane funkcje matematyczne] Ich przegląd zamieszczamy 


w tab. [I] 2]i B] 


Tablica 1: Funkcje trygonometryczne dostępne w bibliotece <cmath>. 


Deklaracja Znaczenie 

double cos (double); cosinus (argument w rad.) 
double sin (double); sinus (argument w rad.) 

double tan (double); tangens (argument w rad.) 
double acos(double); arcus cosinus (wynik w rad.) 
double asin (double); arcus sinus (wynik w rad.) 
double atan (double); arcus tangens (wynik w rad.) 
double atan2(double y, double x); arcus tangens y/x (wynik w rad.) 


Dla przypomnienia, funkcja „sufit” określona jest jako 
[x| =[minieZ:i>2), 
a funkcja „podłoga” zaś jako 


|z| =([maxi€ Z:i< 2). 


?Pełna dokumentacja biblioteki <cmath> w języku angielskim dostępna jest do pobrania ze strony 
http://www.cplusplus.com/reference/clibrary/cmath/. 
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Tablica 2: Funkcja wykładnicza, logarytm, potęgowanie w bibliotece <cmath>. 


Deklaracja Znaczenie 

double exp(double); funkcja wykładnicza 
double log (double); logarytm naturalny 
double /og/0(double); logarytm dziesiętny 
double sqrt (double); pierwiastek kwadratowy 
double pow(double x, double y); zy 


Tablica 3: Funkcje dodatkowe w bibliotece <cmath>. 


Deklaracja Znaczenie 

double fabs (double); wartość bezwzględna 
double ceil (double); „sufit” 

double floor (double); „podłoga” 


2.2 Liczby pseudolosowe 


W bibliotece < csźdlib > znajdują się funkcje służące do generowania liczb pseudolosowych. 
Zawarta jest tam funkcja, dająca za pomocą czysto algebraicznych metod w wyniku liczby, 
które można traktować jako (które wyglądają jak) liczby losowe. 

Generator należy zainicjować przed użyciem funkcją void srand(int z), gdzie z>1 to tzw. 
ziarno. Jedno ziarno generuje zawsze ten sam ciąg liczb. Można także użyć aktualnego czasu 
systemowego do zainicjowania generatora. Dzięki temu podczas każdego kolejnego urucho- 
mienia programu otrzymamy inny ciąg. 


Hinclude <cstdlib > 


¿¡Hinclude <ctime> // tu znajduje sig funkcja time() 
Zd zę 
s srand (time (0)); // za każdym razem inne liczby 
// 
Funkcja 
int rand(); 


generuje całkowite liczby losowe z rozkładu dyskretnego jednostajnego określonego na zbiorze 
{0,1,..., RAND_MAX — 1}. 
Jednostajność oznacza, że prawdopodobieństwo ,,wylosowania” każdej z liczb jest takie samo. 


Uwaga 
RAND_MAX jest stałą zdefiniowaną w bibliotece < cstdlib >. 
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N 


N 


ÉS 


u 


Jeśli chcemy uzyskać liczbę np. ze zbioru (1, 2,3), możemy napisad'| 


cout << (rand() % 3) + 1; 


Z kolei, by uzyskać liczbę rzeczywistą z przedziału [0, 1), piszemy 


cout << ((double)rand() / (double )RAND MAX); 


Uwaga 


Korzystając z własnoręcznie napisanej funkcji generującej liczbę rzeczywistą € [0, 1] 


double /os01() 


( 
return ((double)rand() / (double)RAND_MAX) ; 


) 


łatwo napisać funkcję ,.losujaca” liczbę ze zbioru La, a + 1,...,b), gdzie a,b € Z: 


int JosAB(int a, int b) 


( 
double ab = los0l()x(b-a+l)+a; // liczba rzeczywista 
// z przedziału [a,b+l) 
return (int)(floor(ab)); // ,,podtoga”” z ab 
) 


2.3 Asercje 


Biblioteka cassert udostępnia funkcję o następującej deklaracji: 


void assert(bool); 


Umożliwia ona sprawdzenie dowolnego warunku logicznego. Jeśli nie jest on spełniony, 
nastąpi zakończenie programu. W przeciwnym wypadku nic się nie stanie. 

Taka funkcja może być szczególnie przydatna przy testowaniu programu. Zabezpiecza ona 
m.in. przed danymi, które teoretycznie nie powinny się w danym miejscu pojawić. 


Dla przykładu rozpatrzmy „bezpieczną” funkcję wyznaczającą pierwiastek z nieujemnej 
liczby rzeczywistej. 


Hinclude <cassert > 
Hinclude <cmath> 


double pierwiastek(double x) 


( 


3Podana metoda nie cechuje się zbyt dobrymi własnościami statystycznymi. Szczegóły poznamy jednak do- 
piero na laboratoriach ze statystyki matematycznej w semestrze VI. 
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6 M Dane: 2>0 
7 \\ Wynik: yz 


9 assert(x>=0); // jeśli nie, to błąd — zakończenie programu 
10 return sqgrt(x); 

ul} 

Uwaga 


Sprawdzanie wszelkich warunków za pomocą assert () można wyłączyć globalnie za po- 


mocą dyrektywy 


#define NDEBUG 


umieszczonej na początku pliku źródłowego bądź nagłówkowego. 
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3 Cwiczenia 


Zadanie 6.1. Napisz funkcję parzysta, która sprawdza czy dany argument typu int jest liczbą 
parzystą czy nieparzystą. Zwróć wynik typu bool. 


Zadanie 6.2. Napisz funkcję silnia, która dla danego n € N zwraca wartość 1 x 2x :::xn. 
Zadanie 6.3. Napisz funkcję max, która dla danych a,b, c € Z zwraca ich maksimum. 


Zadanie 6.4. Napisz funkcję med, która znajduje medianę (wartość środkową) trzech liczb 
rzeczywistych, np. med(4, 2, 7) = 4 i med(1, 2,3) = 2. 


Zadanie 6.5. Napisz funkcję nwd zwracającą największy wspólny dzielnik dwóch liczb natu- 
ralnych. 


Zadanie 6.6. Napisz funkcję nww zwracającą najmniejszą wspólną wielokrotność dwóch liczb 
naturalnych. 


Zadanie 6.7. Napisz funkcję o nazwie bmi, która jako argument przyjmuje wzrost (w m) 
i wagę (w kg) pacjenta, a jako wynik zwraca jego wskaźnik masy ciała (BMI), określony jako 
BMI = waga/wzrost”. (Ciekawostka: wg WHO BMI od 18,5 do 25,0 jest uznawana za wartość 
prawidłową.) 


Zadanie 6.8. Napisz funkcję odl, która przyjmuje współrzędne rzeczywiste dwóch punktów x/, 
y1,x2, y2 i zwraca ich odległość euklidesową daną wzorem |x—y| = y (xl — y1)? + (x2 — y2)?. 


Zadanie 6.9. Napisz funkcję odlsup, która przyjmuje współrzędne rzeczywiste dwóch punk- 
tów x1,y1,x2,y2 i zwraca ich odległość w metryce supremum, tj. 


dm(x,y) = max{ |x] — yl|, |x2 — y2|}. 


Zadanie 6.10. Napisz funkcję odlLp, która przyjmuje współrzędne rzeczywiste dwóch punk- 
tów x1,y1,x2,y2 i zwraca ich odległość w metryce L”, gdzie p € |1, oo) jest także parametrem 
funkcji, wg wzoru 


[ix — ylle = 4x1 — y1 + |x2 — y2}. 


Zadanie 6.11. >>2, E” gi jest rozwinięciem funkcji In (1+1) dla (—1, 1] w szereg Taylora. 
Napisz funkcję lognat02, która dla danego x € (0, 2] zwraca przybliżenie wartości ln z, a dla 
x £ (0,2] wartość NaN. 


Zadanie 6.12. Wiemy, że szereg > ;*, = o dia |x| < 1 jest zbieżny i jego suma równa 


1-2) (i)? (47) 
jest v1 + w. Napisz funkcję pierw02, która dla danej liczby rzeczywistej x € [0,2] znajduje 


przybliżenie jej pierwiastka na podstawie podanego wzoru, a dla liczb z g [0,2] zwraca NaN. 


Zadanie 6.13. Dane są rozwinięcia następujących funkcji w szereg Taylora. 


a) e” =), L dla pewnego x € R, 


i=0 n! 


b) sin z = NW = ZAGŃ dla pewnego z € R, 
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c) cosg = X Sr x?” dla pewnego z € R, 


100 (2n)! 


se cera dla pewnego z € (—1, 1). 


d) arcsin z = )) 


Napisz funkcje w C++, które przybliżają wartości powyższych sum bądź zwracają NaN dla 
argumentów poza obszarem zbieżności. 


Zadanie 6.14. Napisz funkcję zaokr, która dla liczby x € R wyznacza jej zaokrąglenie dzie- 
siętne — z dokładnością do podanej liczby cyfr dziesiętnych k — jako |x 10* +0,5|/10*, gdzie 
|u] jest funkcją „podłoga”. 


Zadanie 6.15. Niech dane będą a,b,c typu double. Zmienne te definiują równanie względem 
niewiadomej x € R postaci 


ax? +Hbi+c=0. 


Zaproponuj funkcję w języku C++, która wyznaczy jego rozwiązanie i wypisze wynik na ekran. 


Poprawnie identyfikuj przypadki, gdy dane równanie nie ma rozwiązań w R, a także, gdy nie 
jest ono równaniem kwadratowym. 


Zadanie 6.16. Napisz funkcję implementującą grę w „Zgadulę”. Losuje ona liczbę całkowitą 
z zakresu od 1 do 100. Użytkownik próbuje odgadnąć liczbę wprowadzając swe typy z klawia- 
tury, póki jej nie zgadnie. Za każdym razem otrzymuje komunikat zwrotny, np. ,,za mało” bądź 
„za dużo”. 


Zadanie 6.17. Napisz funkcję implementującą grę w „Zgadulę” zawierającą elementy ,,sztucz- 
nej inteligencji”. Losuje ona liczbę całkowitą z zakresu od 1 do 100. Następnie komputer sam 
próbuje odgadnąć tę liczbę, przy okazji wypisując swe typy na ekranie. Zaproponuj prosty algo- 
rytm, który (nie oszukując!) znajdzie poprawne rozwiązanie w średnio jak najmniejszej liczbie 
kroków. 
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4 Wskazówki do ćwiczeń 


Wskazówka do zadania By znaleźć przybliżenie wartości funkcji albo rozpatrz tylko 
n początkowych wyrazów szeregu, np. n = 30, albo przerwij obliczenia dopiero, gdy moduł 
kolejnego dodawanego wyrazu jest mniejszy niż założone ô, np. 6 = 1076. 
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1 Rozszerzenie wiadomości o funkcjach 


W niniejszym paragrafie rozszerzymy naszą wiedzę związaną z funkcjami o bardziej zaawan- 
sowane zagadnienia. 


1.1 Zasięg zmiennych 


Przyjrzyjmy się zasięgowi definiowanych przez nas zmiennych. Zasięg określa, jak długo dany 
obiekt istnieje i w jaki sposób dane zmienne są widoczne z innych miejsc programu. 

Wszystkie zmienne, które definiujemy wewnątrz każdej funkcji to tzw. zmienne lokalne. 
Tworzone są one, gdy następuje wywołanie funkcji, a usuwane, gdy następuje jej opuszczenie. 

Zmienne lokalne nie są współdzielone między funkcjami. I tak, zmienna x w funkcji f() to 
zmienna inna niż x w funkcji g(). Funkcja nie ma możliwości bezpośredniego odwołania się 
do zmiennej lokalnej innej funkcji, nawet przez się wywoływanej. Zatem jedynym sposobem 
na wymianę danych między funkcjami jest zastosowanie parametrów wejściowych i wartości 
zwracanych. 

Parametry funkcji spełniają rolę zmiennych lokalnych, których wartości są przypisywane 
automatycznie przy wywołaniu funkcji. 


Przyjrzyjmy się rys. 


= int f(int n, int m) 
= int x = 8; { n 
y L---|int y = 4; int x = nxm, ,  |---- m 
Z int z = f(x, y); return x; X 
32 i 


Rysunek 1: Zasięg zmiennych. 


Jak powiedzieliśmy, x po lewej i x po prawej to dwie różne zmienne. Zmienne lokalne n, 
m,x są tworzone na użytek wewnętrzny aktualnego wywołania funkcji f(). Mają one „pomóc” 
funkcji w spełnieniu swego zadania, jakim jest uzyskanie na wyjściu pewnej wartości, która 
jest konieczna do działania bloku kodu po lewej stronie. Zmienne te zostaną skasowane zaraz 
po tym, gdy funkcja zakończy swoje działanie (nie są one już do niczego potrzebne). W tym 
sensie można traktować taką funkcję jako czarną skrzynkę, gdyż to, jakie procesy wewnątrz 
niej zachodzą, nie ma bezpośredniego wpływu na obiekt, który posiłkuje się f© do osiągnięcia 
swoich celów. 

xi y są przekazane do f () przez wartość. To znaczy, że parametry są obliczane przed wywo- 
łaniem (w tym przypadku pobierane są po prostu wartości przechowywane w tych zmiennych), 
a wyniki tych operacji są kopiowane do argumentów wejściowych. Widoczne są one w f() jako 
zmienne lokalne, odpowiednio, n i m. 

Z, wartością zwracaną za pomocą instrukcji return kompilator postępuje w sposób analo- 
giczny, tzn. nie przekazuje obiektu x jako takiego, lecz wartość, którą on przechowuje. Jeśli 
w tym miejscu stałoby np. złożone wyrażenie arytmetyczne albo stała, reguła ta byłaby oczy- 
wiście zachowana. 
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Uwaga 


Warto zapamiętać, że zmienne przekazywane przez wartość są kopiowane, co w przy- 
padku wielu wywołań funkcji na dużych obiektach (należących do złożonych struktur danych, 
zob. wykład X) może być czasochłonne. 


Ponadto, zmienne lokalne nie są przypisane na stałe samym funkcjom, tylko ich wywo- 
łaniom. Obiekty utworzone dla jednego wywołania funkcji są niezależne od zmiennych dla 
innych wywołań. Jest to bardzo istotne w przypadku techniki zwanej rekurencją, która polega 
na tym, że funkcja wywołuje siebie samą. Więcej szczegółów podamy w kolejnej części niniej- 
szego wykładu. 


Uwaga 


W języku C++ dostępne są także zmienne globalne. Jednakże stosowanie ich nie jest zale- 
cane. Nie będziemy ich zatem omawiać. 


1.2 Przekazywanie parametrów przez referencję 


Standardowo parametry wejściowe funkcji zachowują się jak zmienne lokalne — ich zmiana 
nie jest odzwierciedlona ,,na zewnątrz”. 


Rozważmy następujący przykład. 


void zamien(int x, int y) 


( 


int í = x; 
= y; 
y= 1; 


int main() 


( 


int a= 1, m= 2: 

zamien(n, m); // przekazanie parametrów przez wartość 
// (skopiowanie wartości) 

cout << n. <<", << m << endl: 

return 0; 


Wynikiem tej operacji jest, rzecz jasna, napis 1, 2 na ekranie. Funkcja zamien() z punktu 
widzenia main() nie robi zupełnie nic. 


W pewnych szczególnych przypadkach uzasadnione jest zatem przekazywanie parametrów 
w inny sposób — przez referencję (odniesienie). Dokonuje się tego poprzez dodanie znaku & 
po nazwie typu zmiennej na liście parametrów funkcji, wg składni: 
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typź identyfikator 


Przekazanie parametrów przez referencję umożliwia nadanie bezpośredniego dostępu (rów- 
nież do zapisu) do zmiennych lokalnych funkcji wywołującej (być może pod inną nazwą). 
Obiekty te nie są kopiowane; udostępniane są takie, jakie są (lecz być może pod inną nazwą). 

Co ważne, w taki sposób można tylko przekazać zmienną! Przekazanie wartości zwracanej 
przez jakąś inną funkcję, stałej albo złożonego wyrażenia arytmetycznego zakończyłoby się 
niepowodzeniem. 


Tym samym dopiero teraz możemy przedstawić prawidłową wersję funkcji zamien(). 


| void zamien(intś x, intá y) 


ali 


3 int £ = x; 
4 A = Yy; 

5 y = t; 
6| } 


s int main () 


| 


10 int n= 1l,m= 2; 

1 zamien(n, m); // przekazanie parametrów przez referencję 
12 cout << n << ", " << m << endl; 

13 return 0; 


Zmienna n widoczna jest tutaj pod nazwą x. Jest to jednak de facto ta sama zmienna — są 
to dwa różne odniesienia do tej samej, współdzielonej komórki pamięci. 
Uzyskujemy prawidłowy wynik: 2, 1. 


Istnieje jeszcze jedno ważne zastosowanie zmiennych przekazywanych przez referencję. 
Może ono służyć do obejścia ograniczenia związanego z tym, że funkcja za pomocą instrukcji 
return może zwrócić co najwyżej jedną wartość. 

Popatrzmy na poniższy przykład. Jest to funkcja zwracająca część całkowitą ilorazu oraz resztę 
z dzielenia dwóch liczb. 


| void iloreszta(int x, int y, int& iloraz , int& reszta) 


| 


3 iloraz =x l y; 


4 reszta K yi 


Jint main() 


s| { 


9 int n= 7, m= 2; // wejście 


10 int i, r; /x zmienne, które wykorzystamy 
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do przekazania wyniku x/ 
iloreszta(n, m, i, r); 
" " " " 


cout << m <<. "=" ES i << "e" << m << "PT << Fi 
return 0; 


Wynikiem tego programu będzie więc 7=3x2+1. 


1.3 Parametry domyślne funkcji 


Parametry domyślne to argumenty, których jawne pominięcie przy wywołaniu funkcji powo- 
duje, że zostaje im przypisana pewna z góry ustalona wartość. 
Składnia deklaracji parametru domyślnego: 


typ identyfikator = wartość 


Parametry z wartościami domyślnymi mogą być tylko przekazywane przez wartość. Mogą 
się one pojawić tylko na końcu listy parametrów (może być ich wiele). 


Jeśli rozdziela się definicję i deklarację funkcji, parametry domyślne powinny się pojawić 
tylko w jednym miejscu w kodzie programu. Jednakże dla czytelności lepiej jest, gdy pojawiają 
się w deklaracji. 

Oto kilka przykładów prawidłowych deklaracji funkcji: 


— void f(int x=3); 

void f(int x=3, int y=2, int z=5); 
— void f(int x, int y=3, int z=2); 
— void f(int x, int y, int z=2); 
void f(int& x, int y=2); 


Nieprawidłowe deklaracje funkcji: 


— void f(int x, int y=3, int z); 
— void f(int x=3, int y); 
— void fint x, int& y=2); 


Dla ilustracji rozważmy funkcję wyznaczającą pierwiastek liczby rzeczywistej, domyślnie 
o podstawie 2. 


#include <cmath> 


double pierwiastek (double x, double p=2) 
{ // pierwiastek , domyślnie kwadratowy 
assert(p >= 1); 
return pow(x, 1.0/p); 
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int main(void) 

( 
cout << pierwiastek(10) << endl; // 3.1622 
cout << pierwiastek (10, 2.0) << endl; // 3.162278 
cout << pierwiastek(10, 3.0) << endl; // 2.1544 
return 0; 


1.4 Przeciazanie funkcji 


W języku C++ można nadawać te same nazwy (identyfikatory) różnym funkcjom, pod warun- 
kiem, że różnią się one co najmniej typem lub liczbą parametrów. Jest to tzw. przeciążanie 
funkcji (ang. function overloading). Funkcje takie mogą się różnić zwracanym typem, nie jest 
to jednak warunek dostateczny rozróżniania funkcji przeciążonych od siebie. 

Ma to sens, gdy funkcje wykonują podobne (w sensie interpretacyjnym) czynności, jed- 
nakże na danych różnego typu. 


Oto kilka przykładów: 


int modul(int x); 
double modul(double x); 


void f(); 
int f(int x, int y=2); 
int f(int x); // Błąd! — szczególny przypadek powyższego 


bool g(int x, int y); 
char g(int x, int y); // Błąd! — nie różnią się argumentami 


2 Rekurencja 


Z, rekurencją (bądź rekursją, ang. recursion) mamy do czynienia wtedy, gdy w definicji pew- 
nego obiektu pojawia się odniesienie do niego samego. 

Rozważmy najpierw rys. |2| Przedstawia on pewną znaną, zacną damę, trzymającą obraz, 
który przedstawia ją samą trzymającą obraz, na którym znajduje się ona sama trzymająca ob- 
raz... Podobny efekt moglibyśmy uzyskać filmując kamerą telewizor pokazujący to, co właśnie 
nagrywa kamera. 

A teraz popatrzmy na rys.|3|przedstawiający prosty fraktal. Przypatrując się dłużej tej struk- 
turze widzimy, że ma ona bardzo prostą konstrukcję. Wydaje się, że narysowanie kształtu ta- 
kiego drzewa odbywa się następująco. Rysujemy odcinek o pewnej długości. Następnie ob- 
róciwszy się nieco w lewo/prawo, rysujemy nieco krótszy odcinek. W punktach, w których 
zakończyliśmy rysowanie, czynimy podobnie itd. 
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Rysunek 2: Rekurencyjna Mona Lisa. Źródło:  http://www.dominiek.eu/blog/wp- 
content/uploads/2007/11/megamonalisa_recursion.jpg. 


Rysunek 3: Fraktalne drzewo. Źródło: http://diffusedproductions.net/images/fractal-tree2.jpg. 


W przypadku języków programowania mówimy o rekurencji, gdy funkcja wywołuje samą 
siebie. Jednakże takie postępowanie spowodowałoby być może zawieszenie się komputera. 
Istotną cechą rekurencji jest więc warunek stopu, czyli zdefiniowanie przypadku, w którym 
zagłębianie się rekurencyjne zostaje przerwane. 


Uwaga 


Każde wywołanie rekurencyjne powoduje utworzenie nowego zestawu zmiennych lokal- 
nych! 


Rekurencja, jak widzimy, jest bardzo prostą techniką, której opanowanie pozwoli nam two- 
rzyć bardzo ciekawe i często łatwe do zaprogramowania algorytmy. Wiele bowiem obiektów 
(np. matematycznych) jest właśnie ze swej natury zdefiniowana rekurencyjnie. 

W poniższych paragrafach rozważymy kilka ciekawych przykładów takich zagadnień. 
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2.1 Przykład: silnia 
W poniższej definicji silni pojawia się odniesienie do... silni. Zauważmy jednak, że obiekt ten 
jest dobrze określony, gdyż podany został warunek stopu. 

Mí. =- T; 

n! = n(n-1)!dlan>0. 


Funkcję w języku C++ służącą do obliczenia silni można napisać, bazując wprost na defini- 
cji. 


¡int silnia(int n) 


| 


3 assert(n >= 0); 
4 if (n == 0) return 1; 
5 else return nxsilnia(n—1); 


Dla porównania przyjrzyjmy się, jak by mogła wyglądać analogiczna funkcja nierekuren- 
cyjna. 


jint silnia2(int n) 


| 


3 assert(n >= 0); 

4 int w = l; 

5 for (int i=l; i<=n; ++i) 
6 w x= 1; 

7 return w; 


To, którą wersję uważamy za czytelniejszą, zależy oczywiście od nas samych. 


2.2 Przykład: NWD 


Poznaliśmy już algorytm wyznaczania największego wspólnego dzielnika dwóch liczb natural- 

nych. Okazuje się, że NWD można także określić poniższym równaniem rekurencyjnym. Niech 
l<n<m. 

m dla n = 0 

NWD = j 

(n, m) NWD(m modn,n) dlan > 0. 


Przekłada się ona bezpośrednio na następujący kod w C++. 


¡int nvd(int n, int m) 

2l { 

3 assert(n >= m && n >= 0); 
4 if (n == 0) return m; 
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5 else return nwd(m % n, n); 


s| } 


2.3 Przykład: wieże z Hanoi 


Danych jest n krążków o różnych średnicach oraz trzy słupki. Początkowo wszystkie krążki 
znajdują się na słupku nr 1, w kolejności od największego (dół) do najmniejszego (góra). Sytu- 
ację wyjściową dla n = 4 przedstawia rys. 


Rysunek 4: Sytuacja wyjściowa dla 4 krążków w problemie wież z Hanoi. 


Cel: przeniesienie wszystkich krążków na słupek nr 3. 
Zasada #1: krążki przenosimy pojedynczo. 
Zasada #2: krążek o większej średnicy nie może się znaleźć na mniejszym. 


Zastanówmy się, jak by mógł wyglądać algorytm służący do rozwiązywania tego problemu. 
Niech będzie to funkcja Hanoi(k, A, B, C), która przekłada k krążków ze słupka A na słupek C 
z wykorzystaniem słupka pomocniczego B, gdzie A,B,C € {1,2,3}. Wykonywany ruch będzie 
wypisywany na ekranie. 

Aby przenieść n klocków ze słupka A na C, należy przestawić n — 1 mniejszych elementów 
na słupek pomocniczy B, przenieść n-ty krążek na C i potem pozostałe n— 1 krążków przestawić 
z B na C. Jest to, rzecz jasna, sformułowanie zawierające elementy rekurencji. 


Wywołanie początkowe: Hanoi(n, 1, 2, 3). 


Kod w języku C++: 


¡void Hanoi(int k, int A, int B, int C) 
> 
3 if (k>0) // warunek stopu 

4 ( 

5 Hanoi(k-1, A, C, B); 

6 cout << A << "_>" << "C" << endl; 
7 Hanoi(k-1, B, A, ©); 
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Ciekawym zagadnieniem jest wyznaczenie liczby przestawień potrzebnych do rozwiązania 
łamigłówki. W zaproponowanym algorytmie jest to: 


EE = 1, 
L(n) = L(n-1)+1+L(n—1)dlan>0. 


Można pokazać, że jest to optymalna liczba przestawień. 
Rozwiążmy powyższe równanie rekurencyjne, aby uzyskać postać jawną rozwiązania: 


Ln) = 2L(n-1)+1, 
L(n)+1 = 2(L(n-1)+1). 


Zauważmy, że L(n) + 1 tworzy ciąg geometryczny o ilorazie 2. Zatem 


LD+1=2 
L(2)+1 = 4, 


L(3)+1 = 8, 


L(n)+1 Z 2”, 


Więc 
L(n) =2" —1. 

Zagadka Wież z Hanoi stała się znana w XIX wieku dzięki matematykowi E. Lucasowi. Jak 
głosi tybetańska legenda, mnisi w świątyni Brahmy rozwiązują tę łamigłówkę przesuwając 64 
złote krążki. Podobno, gdy skończą oni swe zmagania, nastąpi koniec świata. Zakładając jed- 
nak, że nawet gdyby wykonanie jednego ruchu zajmowało 1 sekundę, to na pełne rozwiązanie 
potrzeba by wtedy aż 2%* — 1 = 18446744073709551615 sekund, czyli około 584542 miliardów 
lat. To ponad 400 razy dłużej niż szacowany wiek Wszechświata! 


2.4 Przykład: liczby Fibonacciego 
Jako ostatni problem rozpatrzmy równanie, do którego doszedł Leonard z Pizy (1202), rozwa- 
żając rozmnażanie się królików. 

Sformułował on następujący uproszczony model szacowania liczby par w populacji. Na 
początku mamy daną jedną parę królików. Para osiąga płodność po upłynięciu jednej jednostki 
czasu od narodzenia. Z każdej płodnej pary rodzi się w kolejnej jednostce jedna para potomstwa. 

Liczbę par królików w populacji w chwili n można opisać za pomocą tzw. liczb Fibonac- 
ciego. 

F(0) = 1, 
R) = 1, 
F(n) = F(n-1)+F(n-2)dlan> 1. 


W wyniku otrzymujemy zatem ciąg 1, 1, 2,3,5,8, 13,21, 34,.... 


Bezpośrednie przełożenie powyższego równania na kod w C++ może wyglądać jak niżej. 
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int fibrek(int n) 
| { 
if (n > 1) return fibrek(n—l)+fibrek(n—2); 
else return l; 


Bezpośredni algorytm rekurencyjny jest jednak nieefektywny! Przypatrzmy się, w jaki spo- 
sób przebiega wywołanie dla fibrek (5): 


fibrek (5); 
fibrek (4); 
fibrek (3); 
fibrek (2); 
fibrek ; 


fibrek 
fibrek(1 


(5) 
(4) 
(3) 
(2) 
(1) 
(0) 
(1) 
fibrek(2); 
(1) 
(0) 
(3) 
(2) 
(1) 
(0) 
(1) 


fibrek 
fibrek 
fibrek 
fibrek 
fibrek 
fibrek 


fibrek(1 


Większość wartości jest niepotrzebnie liczona wielokrotnie. 


Rozważmy, ile jest potrzebnych operacji arytmetycznych (dodawanie) potrzebnych do zna- 
lezienia F, za pomocą powyższej funkcji. Liczbę tę można opisać równaniem: 


L(0) = 0, 
E(1) = 0, 
L(n) = Ln-1)+L(n-2)+1dlan > 1. 


Rozpatrując kilka kolejnych wyrazów, zauważamy, że L(n) = F(n — 1) — 1 dlan > 0. 

Okazuje się, że dla dużych n, L(n) jest proporcjonalne do c” dla pewnej stałej c. Zatem 
liczba kroków potrzebnych do uzyskania n-tej liczby Fibonacciego algorytmem rekurencyj- 
nym rośnie wykładniczo, tj. bardzo szybko. Dla przykładu czas potrzebny do policzenia F4g na 
komputerze autora tych refleksji to 1,4 s, dla F45 to już 15,1 s, a dla F5o to aż 2 minuty i 47 s. 

Tym razem okazuje się jednak, że rozwiązanie rekurencyjne jest nieefektywne. Zatem po- 
zostaje jedynie użycie wersji iteracyjnej, które wykorzystuje tylko co najwyżej n operacji aryt- 
metycznych! 
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3 Cwiczenia 


Zadanie 7.1. Można pokazać, że >, a jest rozwinięciem funkcji e” w szereg Taylora. Napisz 
dwie wersje funkcji Exp, które dla danego x € R znajdują przybliżenie jego eksponensu na 
podstawie podanego wzoru. 


a) W jednej rozpatrz tylko n początkowych wyrazów szeregu, np. n = 30, 
b) W drugiej przerwij obliczenia dopiero, gdy moduł kolejnego dodawanego wyrazu jest 
mniejszy niż założone 6, np. ô = 1076. 


Rada: Niech n i 6 będą parametrami funkcji z wartościami domyślnymi. 


Zadanie 7.2. Napisz kilka przeciążonych wersji funkcji swap, które przestawiają wartości 
dwóch argumentów wejściowych o typach int, double i bool. 


Zadanie 7.3. Napisz nierekurencyjną funkcję służącą do znalezienie n-tej liczby Fibonacciego. 


Zadanie 7.4. Napisz funkcję swap, która za pomocą ciągu przypisań przestawia wartości czte- 
rech zmiennych rzeczywistych (a, b, c, d), tak by na wyjściu otrzymać (c, b, d, a). 


Zadanie 7.5. Napisz funkcję sort, która porządkuje niemalejąco wartości trzech argumentów 
wejściowych (liczby całkowite). 


Zadanie 7.6. Zaproponuj funkcję wyznaczającej wartość tzw. funkcji 91 McCarthy ego. 


TE n — 10 dla n > 100, 
-~ [M(M(n+11)) dla n < 100. 


Ciekawostka: okazuje się, że M (n) = 91 dla każdego n < 101 oraz M(n) = M(n) — 10 dla 
n > 101. 
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1 Dynamiczna alokacja pamięci 


1.1 Organizacja pamięci komputera 


Jak już wiemy, w pamięci operacyjnej komputera przechowywane są zarówno programy jak 
i dane. Jej jednostką podstawową jest komórka o rozmiarze jednego bajta. Każda komórka 
pamięci posiada swój adres. Reprezentowany jest on we współczesnych komputerach jako 32- 
lub 64-bitowa liczba całkowita (bez znaku). 

Z, punktu widzenia każdego programu można wyróżnić następujący podział puli adresowej 
pamięci (w architekturze von Neumanna): 


— kod programu — dane interpretowane są jako instrukcje procesora, 

— stos — przechowywane są wartości zmiennych lokalnych funkcji, 

— sterta — znajdują się dane dynamicznie przydzielane na prośbę programu (zob. dalej), 

— część niedostępna — zarządzana przez system operacyjny (m.in. dane innych progra- 
mów). 


Zatem każdy program przechowuje informacje potrzebne do wykonywania swych czynno- 
Ści na stosie (ang. stack) i stercie (ang. heap). 

Stos jest częścią pamięci operacyjnej, na której dane są umieszczane i kasowane w porządku 
„ostatni na wejściu, pierwszy na wyjściu” (LIFO, ang. last-in-first-out, zob. rys.[1). Umieszcza- 
nie i kasowanie danych na stosie odbywa się automatycznie. Każda wywoływana funkcja two- 
rzy na stosie miejsce dla swoich zmiennych lokalnych. Gdy funkcja kończy działanie, usuwa 
z niego dane (to dlatego zmienne lokalne przestają wtedy istnieć). 


Początek stosu - - - - 


Rysunek 1: Umieszczanie i kasowanie danych na stosie. 
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1.2 Wskaźniki 


Każda zmienna ma przyporządkowaną komórkę (bądź komórki) pamięci, w której przechowuje 
swoje dane, np. zmienna typu int zajmuje 4 takie komórki (4 bajty). 


Fizyczny adres zmiennej (czyli numer komórki) można sprawdzić za pomocą operatora &. 


int x; 
cout << "x znajduje się pod adresem 
// np. Oxe3d30dbc 


" 


<< Mi; 


Istnieje specjalny typ danych do przechowywania informacji o adresach innych zmiennych, 
zwany wskaźnikami. Składnia deklaracji wskaźnika na zmienną typu typ (czyli deklaracji 
zmiennej przechowującej adres w pamięci jakiejś jednostki danych typu typ) jest następująca: 


typ* zmienna; // * oznacza: wskaźnik 


Uwaga 


Istnieje specjalne miejsce w pamięci o adresie O (NULL), do którego odwołanie się powo- 
duje wystąpienie błędu. Często używa się tego adresu np. dla niezainicjowanych wskaźników 
w celu oznaczenia, że nie wskazują one „nigdzie”. 


Na wskaźnikach został określony tzw. operator wyłuskania, «, dzięki któremu możemy 
odczytać, co się znajduje pod danym adresem pamięci. 


Przyjrzyjmy się poniższemu przykładowi. Tworzone są dwie zmienne, jedna typu całkowi- 
tego, a druga wskaźnikowa. Ich rozmieszczenie w pamięci (na stosie, są to bowiem zmienne 
lokalne jakiejś funkcji) przedstawia rys.[2| Każda z tych zmiennych umieszczona jest pod pew- 
nym adresem w pamięci RAM. Początkowy numer komórki można oddczytać za pomocą ope- 
ratora &. 


Listing 1: „Wyłuskanie” danych spod danego adresu. 


int x = 100; 
int» wskx = Sx; 


cout << wskx << endl; // np. 0xdf590544 
cout << «*wskx << endl; // 100 


Wypisanie wartości wskaźnika oznacza wypisanie adresu, na który wskazuje. Wypisanie „wyłu- 
skanego” wskaźnika zaś powoduje wydrukowanie wartości komórki pamięci, na którą pokazuje 
wskaźnik. Jako że zmienna typu int ma rozmiar 4 bajtów, adres następnej zmiennej (wskx) jest 
o 4 jednostki większy od adresu x. 


AiPP-VIII, s. 3. 


wskx 


(int*) 0xd£590544 
Oxdf590548 - - - - 


p (int) 100 


s Oxdf590544 ---- 


Rysunek 2: Zawartość pamięci w programie [I] 


Aby jeszcze lepiej zrozumieć omawiane zagadnienie, rozważmy fragment kolejnego pro- 
gramu. 


Listing 2: Proste operacje z użyciem wskaźników. 


intx w; 
w = Sx; 
xw = l; 
w =&y; 
*w = 2; 


Zawartość pamięci po wykonaniu kolejnych linii kodu przedstawia rys.B] Tym razem za po- 
mocą operatora wyłuskania zapisujemy dane do komórek pamięci, na które pokazuje wskaźnik 
w. 
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pe 
w (int*) ? 
Oxdf590508- - - 
y (int) ? 
Oxdf590504- - - 
x (int) ? 
0xdf590500 - - - 
wW (int*) Oxdf590500 
0xdf590508 - - - 
y (int) ? 
Oxdf590504- - - 
> (int) ? 
--> 0xdf590500- - - 
s (int*) Oxdf590500 
Oxdf590508- - - : 
y (int) ? 
Oxdf590504- - - 
A (int) 1 
«> 0xdf590500- - - 
w (int*) 0xdf590504 
Oxdf590508- - - : 
y (int) ? 
--> 0xdf390504- - - 
> (int) 1 
Oxdf590500- - - 
w (int*) 0xdf590504 
Oxdf590508- - - ; 
y (int) 2 
--> 0xdf590504- - — 
F (int) I 
0xdf590500 - - - 
YN 


Rysunek 3: Zawartość pamięci po wykonaniu kolejnych linii kodu z przykładu [2] 
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1.3 Tablice a wskaźniki 


Zapewne ciekawi nas, w jaki sposób zorganizowana jest pamięć na stosie, gdy deklarowane są 
tablice jednowymiarowe. Możemy to sprawdzić w następujący sposób. 


Listing 3: Organizacja pamięci dla tablic. 
¡int [3] = f1,2,3); 
2) cout << kt; // np. 0x9b64e300 


Y 


Kolejne elementy tablicy w pamięci zawsze następują po sobie (tablica tworzy spójny ciąg 
bajtów), zob. rys.[4] 


AAS 
t [2] (int) 3 
0x9b64e308 - - -- 
t] (int) 2 
0x9b64e304 ---- 
[0] (int) I 
bo?»  0x9b64e300 ---- 


Rysunek 4: Organizacja pamięci w przypadku tablic w przykładzie 3] 


Ponadto, zwróćmy uwagę na to, co się dzieje, gdy wykonamy następujący kod. 


¡cout << St; // np. 0x9bó4e300—znajduje sie tu 
2| cout << t; // 0x9b64e300—wskazuje tu 
3| cout << Er [0]; 7/7 0x9b64e300—tu pierwszy element 


Okazuje się, że zmienną tablicową (dotyczy to tablic o stałym rozmiarze) można traktować jako 
wskaźnik (ale nie odwrotnie). 

Typ int [3] (szczegółowy) jest sprowadzalny do typu int: (uniwersalny). Synonimem int + 
jest int [] , który oznacza „jakaś tablica”, „wskaźnik na jakiś ciąg (być może jednoelementowy) 
danych typu int”. Dlatego też zapis +1 i t[0] jest równoważny. Co więcej, *(1+k) znaczy to 
samo, co £[k]. 

Z, powyższych uwag wynika, że tablicę (dowolnego) rozmiaru można przekazać funkcji 
właśnie za pomocą wskaźnika. Koniecznie jednak trzeba pamiętać, aby także dostarczyć funkcji 
rozmiar tablicy, bowiem wskaźnik to tylko adres pierwszego elementu. 


Przykład: funkcja wyznaczająca sumę elementów tablicy. 


ıı double suma(doublex t, int n) // albo ,,double[] t’? 
| { 

3 double s = 0.0; 

4 for (int i=0; i<n; ++i) 

5 s += £[1]; 

6 return s; 
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int main(void) 
( 
double punkty[4] = { 10.0, 11.0, 12.0, 9.5 ) ; 
cout << suma(punkty, 4); // przekazanie tablicy do funkcji 
return 0; 


Na marginesie, powyższa funkcja może być zapisana również w dwóch następujących po- 
staciach. Nie jest to jednak zalecany sposób pisania kodu, gdyż trudno zrozumieć intencję jego 
autora. 


double suma(doublex t, int n) 


{ 
double s = 0.0; 
for (int i=0; i<n; ++i) 
S += *(t+i); // ROBI SIĘ GORĄCO!!! :-) 
return s; 
) 
double suma(doublex t, int n) 
( 
double s = 0.0; 
for (int i=0; i<n; ++i) 
s += *(t++); // OLABOGA!!! 
return s; 
) 


1.4 Przydział i zwalnianie pamięci ze sterty 


Oprócz Ściśle określonej na etapie pisania programu ilości danych na stosie, można również 
dysponować pamięcią na stercie. Część tej pamięci jest przydzielana (alokowana) dynamicznie 
podczas działania programu za pomocą operatora new. Po użyciu należy ją zwolnić za pomocą 
operatora delete. 


Uwaga 


Zaalokowany obiekt będzie istniał w pamięci nawet po wyjściu z funkcji, w której go 
stworzyliśmy! Dlatego należy pamiętać, aby go usunąć w pewnym miejscu kodu. 


Oto składnia instrukcji służących do alokacji i dealokacji pamięci na jeden obiekt. 


typ* obiekt = new obiekt; // przydział (zwraca wskaźnik) 
// ... 
delete obiekt; // zwolnienie 
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m 
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Tworzonemu pojedynczemu obiektowi można od razu przypisać wartość, np. 


int» n = new int(7); 


Można również przydzielić pamięć na wiele obiektów następujących po sobie, czyli na ta- 


blicę. 

int n = 4; // tutaj już nie musi być stała 

typ* obiekt = new obiekt[n]; // przydział (zwraca wskaźnik) 
// 


delete |] obiekt; // zwolnienie 


Elementom tak tworzonej tablicy nie można niestety przypisać od razu wartości. Trzeba w tym 
celu skorzystać np. z operatora przypisania. 


Oto fragment kodu, który napisał tata Jasia celem zmotywowania go do nauki. 


int n; 
double: godzinyNauki; 


cout << "Ile dni się uczyłeś do kolokwium?"; 
cin >> n; 


godzinyNauki = new double[n]; // utwórz tablicę o n elementach 
cout << "Ile godzin się uczyłeś każdego dnia?"; 
for (int i=0; i<n; i++) 


cin >> godzinyNauki[i|]; 


e rens 
cout << "I tak za mało :—)"; 


delete [] godzinyNauki; // tablica juz nie jest potrzebna dalej 


Zauważmy, że w przeciwieństwie do tablic z poprzedniego wykładu tym razem n nie jest stałą. 
Rozmiar tablicy zostaje ustalony na etapie działania programu. Pobierany jest on z klawiatury, 
uzależniając go od życzenia jego użytkownika. 


2 Łańcuchy znaków 


Do tej pory nie zastanawialiśmy się wspólnie, w jaki sposób można reprezentować w naszych 
programach napisy. Często są one nam potrzebne, np. gdy chcemy zakomunikować coś waż- 
nego użytkownikowi czy też przechować bądź przetworzyć informacje o nienumerycznym cha- 
rakterze. 


AiPP-VIII, s. 8. 


2.1 Kod ASCII 


Pojedyncze znaki drukowane przechowywane są najczęściej jako typ char (ang. character). 
Oczywiście pamiętamy, że tego typu używaliśmy do przechowywania bajtów, czyli 8 bitowych 
liczb całkowitych. 


Istnieje ogólnie przyjęta umowa (standard), że liczbom z zakresu 0—127 odpowiadają Ściśle 
określone znaki, tzw. kod ASCII (ang. American Standard Code for Information Interchange). 
Zestawiaja je tablice [1H4] 

Szczęśliwie w języku C++ nie musimy pamiętać, która liczba odpowiada jakiemu sym- 
bolowi. Aby uzyskać wartość liczbową symbolu drukowanego, należy ująć go w pojedyncze 
cudzysłowy. 


char cl = A; 

char c2 = "nv; // znak nowej linii 

cout << cl << c2; // "A" i przejście do nowej linii 
cout << (int)cl << (char)59 << (int)c2; //  "65;13" 


Jak widzimy, domyślnie wypisanie na ekran zmiennej typu char jest równoważne z wydruko- 
waniem symbolu. Można to zachowanie zmienić, rzutując ją na typ int. 

Ponadto, przyglądając się uważniej tablicy ASCII, warto zanotować następujące prawidło- 
wości. 


Mamy następujący porządek leksykograficzny: A” <'B' <...<'Z' <'a'<'b' 
A ssa" ZY A 

— Symbol cyfry c € (0,1,..., 9) można uzyskać za pomocą wyrażenia ' 0' +c. 

Kod ASCII n-tej wielkiej litery alfabetu łacińskiego to "A' +n-1. 

— Kod ASCII n-tej małej litery alfabetu łacińskiego to 'a' +n-1. 
Zamiana litery 1 na małą literę następuje za pomocą operacji 1+32 == 1+0x20. 
— Zamiana litery 1 na wielką literę następuje za pomocą operacji 1-32 == 1-0x20. 


Pozostałe symbole odpowiadające wartościom liczbowym (0x80—0xFF) nie są określone 
przez standard ASCII. Zdefiniowane są one przez inne kodowania, np. CP-1250 (Windows) 
bądź ISO-8859-2 (Internet, Linux) zawierają polskie znaki diakrytyczne. Jak widzimy, sprawa 
polskich „ogonków ” jest nieco skomplikowana. Zatem na poczatkowym etapie pro 


gramowania uzywajmy tylko liter alfabetu lacinskiego w programach, 
ktore przetwarzaja napisy. 


Uwaga 


Istnieją jeszcze inne standardy kodowania, zwane UNICODE (UTF-8, UTF-16, ... ), w któ- 
rych jednemu znakowi niekoniecznie musi odpowiadać jeden bajt. Obsługa ich jednak jest nieco 


skomplikowana, zatem nie będziemy się nimi zajmować w tym wykładzie. 
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Tablica 1: Kod ASCII cz. I — znaki kontrolne 


DEC HEX Znaczenie DEC HEX Znak 
0 00  NullQA0) 16 10 Data Link Escape 
1 01 Start Of Heading 17 11 Device Control 1 
2 02 Start of Text 18 12 Device Control 2 
3 03 End of Text 19 13 Device Control 3 
4 04 End of Transmission 20 14 Device Control 4 
5 05 Enquiry 21 15 Negative Acknowledge 
6 06 Acknowledge 22 16 Synchronous Idle 
7 07 — Bella) 23 17 End of Transmission Block 
8 08 Backspace (ub) 24 18 Cancel 
9 09 Horizontal Tab 25 19 End of Medium 
10 OA Line Feed (Nr) 26 1A Substitute 
11 OB Vertical Tab (At) 27 IB Escape 
12 OC Form Feed 28 1C File Separator 
13 OD Carriage Return (\n) 29 1D Group Separator 
14 OE Shift Out 30 1E Record Separator 
15 OF Shift In 31 1F Unit Separator 


Tablica 2: Kod ASCII cz. II 
DEC HEX Znaczenie DEC HEX Znak 


32 20 Spacja 48 30 0 
39. (26.71 49 31 1 
4 20 w" 50 32 2 
35. DE 51 33 3 
36 24 $ 52 34 4 
37 25 % 53 35 5 
38 26 8 54 36 6 
300 57. % 55 37 7 
40 28 ( 56 38 8 
41 29 ) 57 39 9 
42 2A * 58 3A : 

A3 OR = 59 3B ; 

4 20, 60 3C < 
45 2D - 61 3D = 
46  2E . 62 3E > 
41 2F / 63  3F ? 
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Tablica 3: Kod ASCII cz. III 


DEC HEX Znaczenie 


DEC HEX Znak 


64 40 0 80 50 P 
65 41 A 81 51 Q 
66 42 B 82 52 R 
67 43 C 83 33 S 
68 44 D 84 54 T 
69 45 E 85 55 - U 
70 46 F 86 56 V 
71 47 G 87 57 W 
72 48 H 88 58 X 
73 49 I 89 59 Y 
74 4A J 90 SA Z 
75 4B K 91 5B I 
76 4C L 92 5C i 
77 4D M 93 5D | 
78 4E N 94 SE ^ 
79 4F O 95 SB ..- 
Tablica 4: Kod ASCII cz. IV 
DEC HEX Znaczenie DEC HEX Znak 
96 60 á ` 112 70 p 
97 6l a 113 71 q 
98 62 b 114 72 r 
99 63 c 115 73 s 
100 64 d 116 74 t 
101 65 e 117 75 u 
102 66 f 118 76 v 
103 67 g 119 77 w 
104 68 h 120 78 x 
105 69 i 121 79 y 
106 6A j 122 7A z 
107 6B k 123 7B 4 
108 6C 1l 124 7C | 
109 6D m 125 TD | 
110 6E n 126 JE ~x 
111 6F o 127 TF Delete 
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2.2 Reprezentacja napisów 


Wiemy już, w jaki sposób obsługiwać pojedyncze znaki. Najprostszym sposobem reprezento- 
wania ciągów symboli drukowanych, czyli napisów, są tablice elementów typu char zakoń- 
czone umownie bajtem o wartości zero (znak *10”). 

Napisy można utworzyć używając cudzystowów (" ... "). Są to jednak tablice tylko do od- 
czytu. Nie wiadomo bowiem, w jakim miejscu w pamięci zostaną one umieszczone. 


char» napisl = "Pewien napis."; // 13 znaków + bajt 0 
// «prawiex równowaznie: 
char napis2[14] = 
"Pp" e” By” Pq? "e? nm * > E 
% > 3 > > 
> 


EJ » % > $ 4.3 ES El $ bd 


> 
Na ay Pe lg Sy Ag M7 „JĄ 


// Znaki w zmiennej napisl są tylko do odczytu! 
napisl [1] *k*; // nie wiadomo co się stanie 
napis2[1] = °k’; /7 ok 


2.3 Operacje na łańcuchach znaków 


Jako że napisy są zwykłymi tablicami, implementacja podstawowych operacji na nich jest dość 
prosta. W niniejszym paragrafie rozważymy kilka z nich, resztę pozostawiając jako ćwiczenie. 
Najpierw przyjrzymy się wypisywaniu. 


char* napis = "Jakiś napis."; // tablica znaków zakończona 


zerem 


// Zatem: 
cout << napis; 


// Jest równoważne: 

int i=0; 

while (napis[i] != *0”*) v/ dopóki nie koniec napisu 
cout << napis[i|]; 
++i; 


Długość napisu można sprawdzić w następujący sposób. 


int dlug(char* napis) 


// Zwraca długość napisu (bez bajtu zerowego). Zgodnie 
// z umową, mimo że jest to tablica , potrafimy jednak 
// sprawdzić , gdzie się ona kończy 
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int i=0; 
while (napis[i] != 0) 
++i; 


return i; 


} 
Ostatnim przykładem będzie tworzenie dynamicznie alokowanej kopii napisu. 

char» kopia(charx napis) 
{ 

int n = dlug (napis), 

char» nowy = new char[n+1]; // o jeden bajt więcej! 

// nie zapominamy o skopiowaniu bajtu zerowego! 

for (int i=0; i<n+l; ++i) 

nowy[ 1] = napis[i|]; 

return nowy; // dalej nie zapomnijmy o dealokacji pamięci 

} 


2.4 Biblioteka cstring 


Biblioteka <cst ring> definiuje wiele funkcji przetwarzających łańcuchy znaków |] Wybrane 
funkcje zestawia tab. 


l Zobacz angielskojęzyczną dokumentację dostępną pod adresem http://www.cplusplus.com/reference/ 


clibrary/cstring/ 
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Tablica 5: Wybrane funkcje biblioteki <cstring>. 


Deklaracja Opis 
int strlen (char: nap); Zwraca długość napisu. 
char: strcpy (char: cel, char: zrodlo ); Kopiuje napisy. Uwaga: długość ta- 


blicy cel nie może być mniejsza niż 
długość napisu zrodlo +1. 


char: strncpy (char cel, char: zrodlo, int n); Kopiuje co najwyżej n znaków 
z jednego napisu do drugiego. 

char: strcat (char: cel, char: zrodlo ); Łączy napisy cel i zrodlo. 

char: strncat (char* cel, char: zrodlo, int n); Dołącza do cel co najwyżej n zna- 
ków z napisu zrodlo. 

int stremp(char* napl, char: nap2); Porównuje napisy. Zwraca 0, je- 


śli są identyczne. Zwraca wartość 
dodatnią, jeśli napl jest większy 
(w porządku leksykograficznym) 
niż nap2. 

int strnemp(char* napl, char: nap2, int n); Porównuje co najwyżej n pierw- 
szych znaków napisów. 

charx strchr (char: nap, char znak); Zwraca podnapis rozpoczynający 
się od pierwszego wystąpienia da- 
nego znaku. 

char: strrchr (char nap, char znak); Zwraca podnapis rozpoczynający 
się od ostatniego wystąpienia da- 
nego znaku. 

char: strstr (char napl, char: nap2); Zwraca podnapis rozpoczynający 
się od pierwszego wystąpienia pod- 
napisu nap2 w napisie napl. 
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3 Cwiczenia 


Zadanie 8.1. Zaimplementuj samodzielnie funkcje z biblioteki <cstring>: strlen (), strcpy 
O, strncpy (O), strcat (), strncat (), stremp(), strstr (), strchr (), strrchr (). 


Zadanie 8.2. Napisz funkcję, która w danym łańcuchu znaków zamieni wszystkie małe litery 
alfabetu łacińskiego na wielkie. 


Zadanie 8.3. Napisz funkcję, która usunie wszystkie znaki odstępów (spacje) z końca danego 
łańcucha znaków. 


Zadanie 8.4. Napisz funkcję, która usunie wszystkie znaki odstępów (spacje) z początku da- 
nego łańcucha znaków. 


Zadanie 8.5. Napisz funkcję, która usunie z danego napisu wszystkie znaki niebędące cyframi 
bądź kropką. 


Zadanie 8.6. Napisz funkcję, która odwróci kolejność znaków w danym napisie. 


Zadanie 8.7. Napisz funkcję, która jako parametr przyjmuje dwa łańcuchy znaków i zwraca 
nowy, dynamicznie alokowany napis będący ich połączeniem, np. dla "ala" i "ola" wynikiem 


2 m 


powinno być " alaola ". 
Zadanie 8.8. Napisz funkcję, która oblicza, ile razy w danym napisie występuje dany znak. 


Zadanie 8.9. Napisz funkcję, która oblicza, ile razy w danym napisie występuje dany inny 
łańcuch znaków, np. w "ababbababa" łańcuch "aba" występuje 3 razy. 


Zadanie 8.10. Napisz funkcję, która dla danej liczby int zwróci dynamicznie alokowany łań- 
cuch znaków, składający się z symboli O lub 1, przechowujący binarną reprezentację argumentu. 


Zadanie 8.11. Napisz funkcję, która dla danego łańcucha znaków, składającego się z symboli O 
lub 1, reprezentującego pewną liczbę w postaci binarnej, zwróci jej wartość jako zmienną typu 
int. 


Zadanie 8.12. Napisz funkcję, która dla danego łańcucha znaków, składającego się z symboli O 
lub 1, reprezentującego pewną liczbę w postaci binarnej, zwróci dynamicznie alokowany napis 
przechowujący jej szesnastkową reprezentację. 


Zadanie 8.13. Palindronf] to ciąg liter, które są takie same niezależnie od tego, czy czytamy 
je od przodu czy od tyłu, np. kobyłamamałybok, możejutrotadamasamadatortujeżom, ikarła- 
patraki. Napisz funkcję sprawdzającą czy dany napis jest palindromem. Zwróć wartość typu 
bool. 


2Zob. http://www.palindromy.pl. 
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1 Reprezentacja macierzy 


Dana jest macierz określona nad pewnym zbiorem 


011 012 *** dim 

Q21 Q22 *** dam 
A= 

Ani dn2 *** Anm 


gdzie n — liczba wierszy, m — liczba kolumn. 


Macierze mogą być wygodnie reprezentowane w języku C++ na dwa sposoby: 


— jako tablice tablic, 


— jako tablice jednowymiarowe. 


Preferujemy tutaj sposób pierwszy. Pomimo nieco zawiłego tworzenia i usuwania tego typu 
obiektów, zapewnia on bardzo wygodny dostęp do poszczególnych elementów. 


int n= ...; // liczba wierszy 
int m = — ...; // liczba kolumn 


sjtyp*x* A; // tablica o elementach typu "typx" 

sA = new typ*[n]; 

s for (int i=0; i<n; ++i) 

7 A[i] = new typim]; 

o // A to n=elementowa tablica tablic m-elementowych 

1|// teraz np. A[OJ[3] to element w I wierszu i IV kolumnie .... 
sl for (int i=0; i<n; ++i) 


4 delete [] A[i|; 
s delete |] A; 


Z, drugiej strony, macierze możemy reprezentować za pomocą jednowymiarowych tablic. 
Łatwo się je tworzy, jednak odwoływanie się do elementów jest dość skomplikowane. 


int n ...; // liczba wierszy 
int m= ...; // liczba kolumn 


y 


s|typ* A; // jednowymiarowa tablica 


sA = new typ[n*m]; 


1|// teraz np. A[lxn+3] to element w II wierszu i IV kolumnie .... 


o| delete [] A; 
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2 Przykładowe algorytmy z wykorzystaniem macierzy 


Omówimy teraz następujące przykładowe algorytmy: 


a) dodawanie macierzy, 
b) mnożenie macierzy, 
c) rozwiązywanie układów równań liniowych metodą eliminacji Gaussa. 


Będziemy rozpatrywać macierze o wartościach rzeczywistych (reprezentowanych przez typ 
double). 


2.1 Dodawanie macierzy 


Dane są dwie macierze A, B typu n x m. Wynikiem dodawania A + B jest macierz: 


011 012 *** dim by ba +e bim 
Q21 Q22 *': dam boj bog + bom 
|| + . 

Ani dn2 *** Anm bni bno RR bnm 
aiy dr aiz +bi2 *** Qim + bim 
as + b21 a22 +b22 *** Q2m + bəm 
An1 + bni An2 + bno *** Anm + brm 


Oto kod funkcji służącej do dodawania macierzy. 


void dodaj(doublexx A, doublexx B, int n, int m, 
doublexx C) 


/x A, B — macierze wejściowe typu n*m 
C — macierz wynikowa typu nxm (pamięć już przydzielona) 
x / 
assert(n > 0 && m > 0); 
for (int i=0; i<n; ++i) 
for (int j=0; j<m; ++j) 
C[i][j] = A[EI[J] + B[EJ[J I: 
} 
Uwaga 


Zamiana kolejności pętli wewnętrznej i zewnętrznej 


for (int j=0; j<m; ++j) 
for (int i=0; i<n; ++i) 
Cli1[71 = A[+117] + BIEJLJI: 
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w przypadku macierzy dużych rozmiarów powoduje, że program może wykonywać się wolniej. 
Na komputerze skromnego autora niniejszego skryptu dla macierzy 10000 x 10000 czas wyko- 
nywania rośnie z 2.7 aż do 24.1 sekundy. Drugie rozwiązanie nie wykorzystuje bowiem w pełni 
szybkiej pamięci podręcznej komputera. 


2.2 Mnożenie macierzy 


Niech A — macierz typu n x m oraz B — macierz typu m x r. Wynikiem mnożenia macierzy 
A - B jest macierz C typu n x r, dla której 


m 
Cij = 5 Qikbkj, 
k=1 


gdzie 1 < i < n, 1 < j < r (por. rys.[T). 


o. 


m 


Rysunek 1: Ilustracja algorytmu mnożenia macierzy. 


A oto kod funkcji służącej do wyznaczenia iloczynu dwóch macierzy. 


void mnoz(doublexx A, doublexx B, int n, int m, int r, 
doublex C) 


/x A — macierz wejściowa typu nxm 
B — macierz wejściowa typu mxr 
C — macierz wynikowa typu nxr (pamięć już przydzielona) 
*/ 
assert(n > 0 && m > 0 && r > 0); 
for (int i=0; i<n; ++i) 
for (int j=0; j<r; ++) 
{ 
C[i][j] = 0; 
for (int k=0; k<m; ++k) 
C[E][5] += A[E][K] * B[k][j]; 
} 
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2.3 Rozwiązywanie układów równań liniowych 
Dany jest oznaczony układ równań postaci 
Ax=b, 


gdzie A jest macierzą typu n x n, b jest n-elementowym wektorem wyrazów wolnych oraz x 
jest n-elementowym wektorem niewiadomych. 
Rozpatrywany układ równań można zapisać w następującej postaci. 


011 Q12 *** din Tı bı 
Q21 Q22 *** dn T2 bą 
Ani dn2 *** Ann Tn bn 


Uproszczona metoda eliminacji Gaussa polega na sprowadzeniu macierzy rozszerzonej 
[A]b] do postaci schodkowej [4'|b']: 


1 1 t / ri $ 
411 Ga Q13 >?  din-1 Ain bi 
zł al ke d / / 
22 023 2,n—1 42m 2 
1 1 t U 
0O 0 a33 >> lana Gn bs 
F + t 

0 0 0 © dn_ln-=1 În-in n—1 
1 1 
0 0 0 -.. 0 Śr b, 


Dokonuje się tego za pomocą następujących tzw. operacji elementarnych: 


a) pomnożenie dowolnego wiersza macierzy rozszerzonej przez niezerową stałą, 
b) dodanie do dowolnego wiersza kombinacji liniowej pozostałych wierszy. 


Następnie uzyskuje się wartości wektora wynikowego x korzystając z eliminacji wstecznej: 


dla kolejnych i = n,n — 1,...,1. 


Rozpatrzmy przykładową implementację tej metody w języku CH. 


N 


void schodkowa(doublexx A, doublex b, int n); 
void eliminwst(doublexx A, doublex b, int n, doublex x); 


void gauss(doublexx A, doublex b, int n, doublex x) 
5| { 
assert(n > 0); 
schodkowa(A, b, n); 
eliminwst(A, b, n, x); 
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Funkcja sprowadzająca macierz do postaci schodkowej: 


void schodkowa(doublex* A, doublex b, int n) 
{ 
for (int i=0; i<n; ++i) 


( /x dla każdego wiersza */ 


assert(A[i][i]!=0.0);// jeśli nie=>macierz osobliwa 


// assert(fabs(A[i]j[i]) > le—10); // lepiej 


for (int j=i+1; j<n; ++j) 
{ // wiersz j:=(wiersz j)—A[j][i]/A[i][i]*»(wiersz 


for (int k=i; k<n; ++k) 
A[jl[k] —=A[JI[E]/A[iE][E]*A[E][Kk]; 


BLI] —= ALĻJ]Li]/ALi][i]*b[i]; 


i) 


Eliminacja wsteczna: 


void eliminwst(doublexx A, doublex b, int n, doublex x) 


( 


for (int i=n—1; i>=0; —i) 

{ 
x[i] = b[i]; 
for (int k=i+l; k<n; ++k) 
{ 


x[i] —= A[i][k] * x[k]; 


) 
x[i] /= A[LJ[E]; 
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3 Cwiczenia 


Zadanie 9.1. Dana jest macierz A typu n x m o wartościach rzeczywistych oraz liczba k € R. 
Napisz funkcję, która wyznaczy wartość kA, czyli implementującą mnożenie macierzy przez 
skalar. 


Zadanie 9.2. Dana jest macierz A typu n x m o wartościach rzeczywistych oraz wektor b € 
R”. Napisz funkcję, która zwróci macierz [A|b), czyli A rozszerzoną o nową kolumnę, której 
wartości pobrane są z b. 


Zadanie 9.3. Dana jest macierz A typu n x m o wartościach rzeczywistych oraz wektor 
b e R”. Napisz funkcję, która zwróci macierz A rozszerzoną o nowy wiersz, którego war- 
tości pobrane są z b. 


Zadanie 9.4. Dana jest macierz A typu 2 x 2 o wartościach rzeczywistych. Napisz funkcję, 
która zwróci wyznacznik danej macierzy. 


Zadanie 9.5. Dana jest macierz A typu 3 x 3 o wartościach rzeczywistych. Napisz funkcję, 
która zwróci wyznacznik danej macierzy. 


x Zadanie 9.6. Dana jest macierz kwadratowa A o 4 wierszach i 4 kolumnach zawierająca 
wartości rzeczywiste. Napisz rekurencyjną funkcję, która zwróci wyznacznik danej macierzy. 
Skorzystaj wprost z definicji wyznacznika. Uwaga: taka metoda jest zbyt wolna, by korzystać 
z niej w praktyce. 


Zadanie 9.7. Napisz funkcję, która rozwiązuje układ 2 równań liniowych korzystając z metody 
Cramera. Poprawnie identyfikuj przypadki, w których dany układ nie jest oznaczony. 


Zadanie 9.8. Napisz funkcję, która rozwiązuje układ 3 równań liniowych korzystając z metody 
Cramera. Poprawnie identyfikuj przypadki, w których dany układ nie jest oznaczony. 


Zadanie 9.9. Dana jest macierz A o wartościach całkowitych. Napisz funkcję, która zwróci jej 
transpozycję. 


Zadanie 9.10. Dana jest macierz A typu n x m o wartościach całkowitych. Napisz funkcję, 
która dla danego 0 < i < ni0 < j < m zwróci podmacierz powstałą przez usunięcie z A 
i-tego wiersza i j-tej kolumny. 


Zadanie 9.11. Dana jest kwadratowa macierz A o wartościach całkowitych. Napisz funkcję, 
która sprawdzi, czy macierz jest symetryczna. Zwróć wynik typu bool. 


Zadanie 9.12. Dana jest macierz kwadratowa A o wartościach rzeczywistych typu nxn. Napisz 
funkcję, która zwróci jej ślad, określony jako 


tr(4) = an + a22 +*** + ann = AC 
i 


Zadanie 9.13. Dla danej macierzy kwadratowej A napisz funkcję, która zwróci jej diagonalę 
w postaci tablicy jednowymiarowej. 
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Zadanie 9.14. Dla danej macierzy kwadratowej A napisz funkcję, która zwróci jej macierz 
diagonalną, czyli macierz z wyzerowanymi wszystkimi elementami poza przekątną. 


Zadanie 9.15. Kwadratem łacińskim stopnia n nazywamy macierz kwadratową typu n x n 
o elementach ze zbioru (1,2,...,n)j taką, że żaden wiersz ani żadna kolumna nie zawierają 
dwóch takich samych wartości. Napisz funkcję, która sprawdza, czy dana macierz jest kwadra- 
tem łacińskim. Zwróć wynik typu bool. 


Zadanie 9.16. Kwadratem magicznym stopnia n nazywamy macierz kwadratową typu n x n 
o elementach ze zbioru liczb naturalnych taką, że sumy elementów w każdym wierszu, w każdej 
kolumnie i na każdej z dwóch przekątnych są takie same. Napisz funkcję, która sprawdza, czy 
dana macierz jest kwadratem magicznym. Zwróć wynik typu bool. 
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4 Wskazówki do ćwiczeń 


Wskazówka do zadania 9.6! Wyznacznik macierzy 4 x 4 można policzyć ze wzoru 
4 . . 
det A = ) (Day det Aj; 
i=1 
gdzie j jest dowolną liczbą ze zbioru (1,2, 3,4), a A; jest podmacierza 3 x 3 powstałą przez 


opuszczenie ¿-tego wiersza i j-tej kolumny. Do wyznaczenia det A, ; można skorzystać z nieco 
zmodyfikowanej funkcji z poprzedniego zadania, której należy przekazać A, i oraz j. 
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1 Struktury w języku C++ 


Do tej pory omawialiśmy następujące ogólne typy zmiennych: zmienne skalarne, tablicowe 


i wskaźnikowe. 


W języku C++ możemy tworzyć własne typy złożone będące reprezentacją iloczynu ska- 


larnego różnych innych zbiorów (typów). Są to tzw. struktury. Oto składnia ich definicji: 


struct NazwaNowegoTypu 
( 
typl nazwaPolal; 
typ2 nazwaPola2; 
Lata 
typN nazwaPolaN; 


|: // konieczny średnik! 


Zmienne typu złożonego deklarujemy w standardowy sposób, czyli np. 


nazwaStruktury identyfikator; 


Do poszczególnych pól struktury możemy odwołać się za pomocą kropki, np. 


identyfikator .nazwaPola 


Pola zmiennej typu złożonego traktujemy jak zwykłe zmienne odpowiednich typów. 


Przykład: struktura reprezentująca zbiór N x R. 


struct NR 
{ 
int polen, 
double poler; 
E 


int main() 


( 
NR x; // tzn. xENxR 
x.polen = 1; 
x.poler = 3.14; 


cout << x; // BŁĄD! operacja niezdefiniowana 
cout << x.polen; // OK 


return 0; 


Oczywiście można też tworzyć tablice i wskaźniki do obiektów tego typu: 
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int main () 

( 
NR x[3]; // tablica 

100; 

100.0; 


x[0].polen 
x[0].poler 
EAS tes 
NRx wska = «x[1]; 
EL sica AED 
return 0; 


Rozważmy teraz następujące funkcje: 


void fI(NR a) 


// a przekazana przez wartość — kopiowana 
cout << a.polen << " " << a.poler; 
void f2(NR£ a) 
// a przekazana przez referencję — niekopiowana 
cout << a. polen << " " << 8.poler; 


Jak już wiemy, przekazanie zmiennej przez referencję jest wydajniejsze niż podanie jej przez 
wartość. W pierwszym przypadku nie zostanie utworzona jej kopia. Jednakże zmiany wprowa- 


dzone w przekazanych obiektach będą widoczne na zewnątrz funkcji. 
Idąc dalej tym śladem, rozpatrzmy 2 kolejne funkcje. 


void f3(const NR& a) 


// a przekazana przez referencję — niekopiowana 
// + zabezpieczona przed zmianą 
cout << a.polen << " | << a.poler ,; 
void f4(NRx a) 
// a przekazana przez wskaźnik — niekopiowana 
cout << (xa).polen << " * << (xa), poler; 
// równoważny zapis (!) 
cout << a—>polen << " " << a—>poler; 
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W funkcji /3() dodatkowo zabezpieczyliśmy obiekt przed przypadkową zmianą za pomocą 
modyfikatora const. 

Jak widzimy, przekazanie zmiennej przez wskaźnik w tym kontekście jest równoważne 
przekazaniu jej przez referencję. Dostęp do pól struktury przekazanej przez wskaźnik możemy 
uzyskać za pomocą operatora —>, który jest bardzo wygodny. 


2 Podstawowe abstrakcyjne struktury danych 


W niniejszym paragrafie omówimy następujące dynamiczne abstrakcyjne struktury danych: 


a) listy jednokierunkowe, 
b) stosy, 

c) kolejki, 

d) kolejki priorytetowe, 

e) listy dwukierunkowe, 
f) drzewa binarne. 


Służą one do przechowywania różnego rodzaju danych. Każda z nich charakteryzuje się innymi 
właściwościami. Np. lista jednokierunkowa pozwala bardzo szybko dodawać nowe elementy 
do zbioru, a w kolejce priorytetowej mamy łatwy dostęp do danych uporządkowanych. 

Dla uproszczenia w omawianych strukturach danych będziemy przechowywali jeden ele- 
ment typu int. 


2.1 Lista jednokierunkowa 


Lista kierunkowa (ang. linked list), podobnie jak tablica, służy do przechowywania elementów 
tego samego typu. W odróżnieniu jednak od tablicy nie ma ona z góry ustalonego rozmiaru. 
Pozwala na efektywne wstawianie 1 usuwanie elementów. Odbywa się to kosztem czasu ich 
wyszukiwania. 

Dane przechowywane są w węzłach następującej postaci: 


struct wezel 
int elem; // element(y), który przechowujemy w węźle 
wezelx nast; // wskaźnik na następny element 


E 


Węzły umieszczone są w różnych (dowolnych) miejscach w pamięci. Każdy z nich jest osobno 
przydzielany dynamicznie. 

Lista składa się z węzłów połączonych ze sobą w kierunku od pierwszego do ostatniego 
elementu. Zatem dodatkowo należy zapamiętać, gdzie leży pierwszy element: 


wezel* pocz; // tzw. głowa listy 


Schemat przykładowej listy jednokierunkowej przechowującej elementy 3, 5 oraz 7 przed- 
stawia rys.|l| Zwróćmy uwagę na to, w jaki sposób węzły mogą być rozlokowane w pamięci. 
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Y 
(int) 5 
(wezel*) Oxff0001cd--- 
ue»  Oxffabc508 - - - = . : 
: (int) 3 
: (wezel*) Oxffabc508 : 
s::ł:>  Oxffabc500 - - - -> : : 
A TETTEN A NAE AN A ; 
w a 
: (int) 7 
A (wezel*) 0 mo . 
=> OXff000lcd — — — : 
NAAS : 
VETAS : 
ocz * : 
NA = = (wezel JOR HABeSd0 E 


Rysunek 1: Przykładowa lista jednokierunkowa przechowująca elementy 3,5,7. Schemat gra- 
ficzny i organizacja pamięci. 


Przyjrzymy się implementacji następujących operacji: 


— sprawdzenie, czy dany element występuje na liście, 
wstawienie elementu na początek listy, 
wstawienie elementu na koniec listy, 


— usunięcie elementu z początku listy, 
— usunięcie elementu z końca listy. 


Uwaga 


Należy zwrócić szczególną uwagę na przypadek, gdy lista początkowo jest pusta (gdy pocz 
== NULL)! 


2.1.1 Wyszukiwanie elementu 


Na tys. [2| i B| znajdziemy ilustrację, krok po kroku, w jaki sposób przebiega wyszukiwanie 
elementów 5 oraz 2 na liście zawierającej (3,5,7). 
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Rysunek 2: Wyszukiwanie elementu 5 w liście jednokierunkowej. 
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wsk e NULL 


Rysunek 3: Wyszukiwanie elementu 2 w liście jednokierunkowej. 
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W pierwszym przypadku element zostaje odnaleziony już w drugim kroku. W kolejnym 
stwierdzamy, że nie ma go wcale na danej liście. 

Zauważmy, że korzystamy tutaj z dodatkowego wskaźnika, za pomocą którego poruszamy 
się po odwiedzanych węzłach. Wskaźnik ten rozpoczyna swoje poszukiwania od głowy listy. 
Jest to bowiem jedyny dostępny bezpośrednio element. Do każdego kolejnego dostajemy się za 


pomocą pola nast. 


Oto iteracyjna wersja funkcji formalizującej powyższe kroki. 


N 


bool szukaj(wezelx pocz, int x) 


( 
wezelx* wsk = pocz; 
while (wsk != NULL) 
( 
if (wsk—>elem == x) 
return true; // znaleziony 
else 
wsk = wsk-—>nast; // przejdź do następnego węzła 
} 
return false; // doszliśmy do końca listy 
} 


Wywołanie: 


N 


szukaj(pocz, x); 


Podany wyżej algorytm można również zapisać równoważnie w formie rekurencyjnej. 


bool szukaj2(wezelx wsk, int x) 


( 
if (wsk == NULL) 
return false; // doszliśmy do końca listy 
else if (wsk—>elem == x) 
return true; // znaleziony 
else 
return szukaj2(wsk—>nast, x); // szukaj dalej 
} 
Wywołanie: 


szukaj? (pocz, x); 


2.1.2 Wstawianie elementu 


Najpierw zajmijmy się wstawianiem elementu na początek listy. Rys. [4] ilustruje kolejne kroki 
potrzebne do utworzenia listy (6,3,5,7) z listy (3,5,7). Zwróćmy uwagę, że po dokonaniu tej 
operacji zmienia się głowa listy. 
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Rysunek 4: Wstawianie elementu 6 na początek listy jednokierunkowej. 
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Oto kod stosownej funkcji w języku C++. Zauważmy, że pierwszym parametrem jest refe- 
rencja na zmienną wskaźnikową. Przekazujemy tutaj głowę listy, która będzie zmieniona przez 
tę funkcję. 


¡void wstawpocz(wezelx*8 pocz, int x) 

2 { 

3 wezelx* nowy = new wezel, 

4 nowy—>elem = x; 

3 

6 nowy->nast = pocz; /* wskazuj na to, na co wskazuje pocz, 
może być NULL x/ 

7 

8 pocz = nowy; /x teraz lista zaczyna się 

9 od nowego węzła */ 

10) ) 


Rozważmy teraz, w jaki sposób wstawić element na koniec listy. Rys. 5]i[ójilustrują, w jaki 
sposób utworzyć listę (3,5,7,1), mając daną listę (3,5,7). Zauważmy, że pierwszym krokiem, 
jaki należy wykonać, jest przejście na koniec listy. Tutaj korzystamy z dwóch dodatkowych 
wskaźników. 


Poniższa funkcja implementuje interesującą nas operację. Zauważmy, że pierwszym pa- 
rametrem jest referencja do wskaźnika. Jeżeli lista jest pusta, przekazana głowa może zostać 
zmieniona. 


| void wstawkon(wezelxóż pocz, int x) 


| { 
3 wezelx nowy = new wezel; 
4 nowy->elem = x; 


5 nowy->nast = NULL; 


7 if (pocz == NULL) // czy lista pusta? 


8 pocz = nowy; 

9 else 

o { 

1 wezelx* wsk = pocz; 

2 while (wsk-—>nast != NULL) 

3 wsk = wsk->nast; // przejdź na ostatni element 
4 

5 wsk=>nast = nowy; // nowy ostatni element 

6 ) 

1) 


A oto równowazna wersja rekurencyjna. 
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Rysunek 5: Wstawianie elementu 1 na koniec listy jednokierunkowej cz. I. 
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Rysunek 6: Wstawianie elementu 1 na koniec listy jednokierunkowej cz. II. 


void wstawkon2(wezelxż wsk, int x) 


( 
if (wsk != NULL) 
wstawkon2(wsk-—>nast, x); // idź dalej 
else 
( 
wsk = new wezel, 
wsk—>elem = x; 
wsk->nast = NULL; 
) 
) 


2.1.3 Usuwanie elementu 


Podobnie jak wyżej, i tutaj rozpatrzymy dwa przypadki, w których usuwany element znajduje 


się na początku bądź na końcu listy. 


Rys. [7] przedstawia możliwy sposób usuwania elementów z początku listy jednokierunko- 
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Rysunek 7: Usunięcie elementu z początku listy jednokierunkowej. 


AiPP-X, s. 13. 


wej. Implementacje tej operacji w C++ przedstawia poniższy kod. Zauważmy, że po jej do- 
konaniu możliwe jest osiągnięcie stanu glowa == NULL, co oznacza, że usunięty element był 
jedynym. 


void usunpocz(wezelxK pocz) 


( 
if (pocz != NULL) 
( 
wezelx* wsk = pocz; 
pocz = pocz—>nast; 
delete wsk; 
} 
) 


Zastanówmy się teraz, w jaki sposób dokonać wykasowania ostatniego elementu. 
Rys. [8Jprzedstawia krok po kroku kolejne działania. Formalizuje je następujący kod w C++. 


void usunkon(wezelx8 pocz) 


( 
if (pocz == NULL) 
return ; // lista pusta 
else if (pocz—>nast == NULL) 
{ // tylko jeden element 
delete pocz; 
pocz = NULL; 
} 
else 
{ // > 1 element 
wezelx* wsk = pocz; 
// przejdź na przedostatni element 
while (wsk—>nast—>nast != NULL) 
wsk = wsk—>nast , 
delete wsk—>nast ; 
wsk—>nast = NULL; 
} 
} 
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Rysunek 8: Usunięcie elementu z końca listy jednokierunkowej. 
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Omawianą procedurę można zapisać również w postaci rekurencyjnej. 


void usunkon2(wezelx8£ wsk) 
{ 
if (wsk == NULL) return; // lista była pusta 
else if (wsk—>nast == NULL) // jesteśmy na końcu 
{ 
delete wsk; 
wsk = NULL; 
) 
else // idziemy dalej 
usunkon2(wsk->nast); 


2.1.4 Uwaga na temat wydajności 


Porównajmy na koniec dwie struktury danych: listę jednokierunkową oraz zwykłą tablicę jed- 
nowymiarową, przechowujące n elementów. Poniższa tabela zestawia liczbę elementów, które 
należy rozpatrzeć, aby móc zrealizować podstawowe operacje. 


Operacja Tablica Lista 
Dostęp do i-tego elementu 1 i 
Wyszukiwanie <n <n 
Wstawianie na początek n 1 
Wstawianie na koniec n n (*) 
Kasowanie z początku n 1 
Kasowanie z końca n n 


(*) Liczbę operacji potrzebnych do wstawienia elementu na koniec listy można zredukować 
do 1 jeśli dodatkowo będziemy przechowywać wskaźnik na koniec listy. Wymaga to jednak 
przerobienia wszystkich funkcji. 


Zauważmy, że dostęp do poszczególnych elementów w tablicy jest natychmiastowy. Z kolei 
lista ma tę zaletę, iż szybko pozwala zwiększać bądź zmniejszać swój rozmiar. Nie wymaga 
ona bowiem kopiowania wszystkich elementów za każdym razem. 


2.2 Stos 


Stos (ang. stack) jest strukturą danych typu LIFO (ang. last-in-first-out), tzn. elementy są zdej- 
mowane z niej w kolejności odwrotnej niż były na niej umieszczane. 
Stos udostępnia tylko dwie operacje: 


— umieść (ang. push) — wstawia element na początek stosu, 
— wyjmij (ang. pop) — usuwa 1 zwraca element z początku stosu. 


Bardzo łatwo jest zaimplementować stos z użyciem już poznanych algorytmów dla listy 
jednokierunkowej. Dlatego pozostawiamy go jako ćwiczenie. 
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2.3 Kolejka 


Kolejka (ang. queue) jest strukturą danych typu FIFO (ang. first-in-first-out), czyli elementy są 
udostępniane w tej samej kolejności, w jakiej były w niej umieszczane. Udostępnia ona dwie 
następujące operacje: 


— umieść (ang. enqueue) — wstawia element na koniec kolejki, 
— wyjmij (ang. dequeue) — usuwa i zwraca element z początku kolejki. 


Ze względów wydajnościowych kolejkę dobrze jest implementować jako listę jednokierun- 
kową, w której przechowuje się dodatkowo wskaźnik na ostatni element. Jej implementację 
pozostawiamy jako ćwiczenie. 


2.4 Kolejka priorytetowa 


Jako kolejne ćwiczenie pozostawiamy także implementację tzw. kolejki priorytetowej (ang. 
priority queue), która jest modyfikacją listy jednokierunkowej. Przechowuje ona elementy wg 
porządku <= (bez duplikatów). Udostępnia następujące operacje: 


— wstawienie elementu wg porządku <=, 
— pobranie elementu najmniejszego, 
— usunięcie elementu najmniejszego. 


2.5 Lista dwukierunkowa 


Lista dwukierunkowa pozwala na bardzo szybkie wstawianie i usuwanie elementów zarówno 
z końca, jak i z początku listy. Dodatkowo pozwala na przeglądanie elementów w kierunku 
przeciwnym. Dane przechowywane są w węzłach następującej postaci: 


struct wezel 
int elem; // element(y), który przechowujemy w węźle 
wezelx nast; // wskaźnik na następny element 
wezelx poprz; // wskaźnik na poprzedni element 


| 


Należy rzecz jasna zapamiętać wskaźnik na pierwszy (pocz) i wskaźnik na ostatni (kon) 
element listy dwukierunkowej. 


Przykładowa lista dwukierunkowa (3,5,7) przedstawiona jest na rys. g 
Jako że jej implementacja jest dość podobna do omawianych już konstrukcji, pozostawiamy 
ją jako ćwiczenie. 


2.6 Drzewo binarne 


Drzewo poszukiwań binarnych (BST, ang. binary search tree) jest dynamiczną abstrakcyjną 
strukturą danych, która przechowuje dane uporządkowane względem relacji <=. Umożliwia ono 
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Rysunek 9: Przykładowa lista dwukierunkowa. 


często dość szybkie wyszukiwanie elementów. Dla danych losowych oczekiwana liczba opera- 
cji koniecznych do znalezienia danej wartości jest rzędu log, n (logarytm jest funkcją, która 
dość wolno zwiększa swe wartości wraz ze wzrostem wartości argumentu, np. log, 1000 = 10, 
a log, 10000 = 13). Niestety, w przypadku pesymistycznym liczba ta rośnie do n, gdzie n 
to rozmiar przechowywanego zbioru danych. Na zajęciach w III semestrze dowiemy się, jak 
równoważyć drzewa w taki sposób, aby wyszukiwanie było zawsze efektywne. 


Dane w drzewie binarnym przechowywane są w węzłach zdefiniowanych następująco: 


struct wezel 
int elem; // element(y), który przechowujemy w węźle 
wezelx lewy; // wskaźnik na lewe poddrzewo (lewego potomka) 
wezelx prawy; // wskaźnik na prawe poddrzewo (prawego 
potomka ) 


E 


Zatem każdy węzeł ma co najwyżej dwóch potomków. 

Drzewo jest strukturą acykliczną, tzn. wychodząc z dowolnego węzła za pomocą wskaźni- 
ków, nie da się do niego wrócić. Ponadto, jeśli istnieje Ścieżka pomiędzy węzłami v i w, jest 
ona określona jednoznacznie. 

Początek drzewa określa węzeł zwany korzeniem (ang. root). Nie da się do niego dojść 
z żadnego innego węzła (inaczej: nie ma on rodzica). Ponadto, drzewo jest spójne, tzn. z korze- 
nia można dojść do każdego innego węzła. Węzły, które nie mają potomków, zwane są inaczej 
liśćmi. 

Rozpatrzmy dowolny węzeł v. Ważną cechą drzewa jest to, że 

a) wszystkie elementy w jego lewym poddrzewie są nie większe niż element przechowy- 
wany w węźle v, 

b) wszystkie elementy w jego prawym poddrzewie są większe niż element przechowywany 
wy. 


Często rozpatruje się drzewa (i tak czynimy w tym przypadku), w których zabronione jest prze- 
chowywanie duplikatów elementów. 


Rozważymy następujące operacje na drzewie binarnym. 
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N 


— wyszukiwanie zadanego elementu, 
wypisywanie elementów wg porządku <=, 
— wstawianie, 


— usuwanie. 


2.6.1 Wyszukiwanie elementu 


Rys. [TOJilustruje krok po kroku wyszukiwanie elementu 2 w drzewie składającym się z elemen- 
tów 41, 2, 3,5, 6,8). 


Uwaga 

Zwróćmy uwagę, że istnieje wiele postaci drzew binarnych, które mogłyby służyć do prze- 
chowywania powyższego zbioru wartości. W szczególności drzewo takie mogłoby się zredu- 
kować do listy, np. jeśli żaden z węzłów nie miałby lewego bądź prawego poddrzewa. W tym 
przypadku wyszukiwanie może wymagać przejrzenia wszystkich elementów. 


Oto wersja iteracyjna funkcji służącej do wyszukiwania elementów o zadanej wartości. 


bool szukaj(wezelx korzen, int x) 


( 


wezelx wsk = korzen, 
while (wsk != NULL) 
{ 


if (x == wsk—>elem) return true; 
else if (x < wsk—>elem) wsk = wsk—>lewy ; 
else wsk = wsk—>prawy ; 
) 
return false; // nie znaleziono 


Poniżej znajduje się równoważna funkcja rekurencyjna. 


bool szukaj2(wezel*x wsk, int x) 


( 


if (wsk == NULL) return false; 

else if (x == wsk—>elem) return true; 

else if (x < wsk—>elem) return szukaj2(wsk—>lewy ,x); 
else return szukaj2(wsk=>prawy ,x); 
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Rysunek 10: Wyszukiwanie elementu 2 w przykładowym drzewie binarnym. 
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2.6.2 Wypisywanie wszystkich elementów wg porządku 


Procedura wypisująca wszystkie elementy przechowywane w drzewie uwzględniająca porządek 
<= jest ze swej natury rekurencyjna. Skoro wartości w lewym poddrzewie danego węzła są 
mniejsze (przypominany, omijamy duplikaty) od wartości w danym węźle, dajmy na to, v, przed 
wypisaniem wartości w v należy wypisać wartości w lewym poddrzewie v. 

Implementacja wyszukiwania jest dość prosta. 


void wypisz(wezelx* wsk) 


( 


3 if (wsk == NULL) return; // tu nie ma nic do roboty 


N 


5 Wypisz(wsk—>lewy); // najpierw lewe poddrzewo 
6 cout << wsk—>elem; // teraz wypisywanie 
7 wypisz(wsk—>prawy); // a potem prawe poddrzewo 


2.6.3 Wstawianie elementu 


Wstawianie elementu 4 do przykładowego drzewa binarnego ilustruje rys. 


Najprościej wstawić element jako liść. Zatem najpierw należy znaleźć węzeł, do którego 
nowy element będziemy doczepiać jako potomka. Oto wersja rekurencyjna stosownej funkcji: 


¡void wstaw(wezelx8£ wsk, int x) 
| 
3 if (wsk == NULL) 

4 { // tutaj będzie nowy liść 


5 wsk = new wezel, 

6 wsk—>elem = x; 

7 wsk->lewy = NULL; 

8 wsk=>prawy = NULL; 

9 ) 

10 else if (x == wsk—>elem) return; // nic nie rób 

1 else if (x < wsk—>elem) return wstaw(wsk—>lewy, x); 
12 else return wstaw(wsk-—>prawy, x); 


2.6.4 Usuwanie elementu 


Usuwanie elementu z drzewa jest dość skomplikowane. Musimy rozpatrzeć trzy przypadki: 
a) liść — po prostu usunąć, 


b) węzeł z jednym potomkiem — zastąpić potomkiem, 
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korzen 


Rysunek 11: Wstawianie elementu do przykładowego drzewa binarnego. 


c) węzeł v z dwoma potomkami — zastąpić lewym (prawym) dzieckiem (ozn. w). Lewe 
(prawe) poddrzewo w pozostaje bez zmian. Prawe (lewe) poddrzewo v staje się prawym 
(lewym) poddrzewem w. Prawe (lewe) poddrzewo w zostaje doczepione za elementem 
najmniejszym (największym) w prawym (lewym) poddrzewie v. 


Zainteresowanym przedstawiamy przykładową implementację. 


void usun(wezelx8 wsk, 


( 
if (wsk == NULL) 
else if (x 
else if (x 
else { 

if (wsk-—>lewy 


wezelx usuwany 
wsk = wsk=>prawy; 


< wsk— 
> wsk—>elem) usun(wsk=>prawy, x); 


int x) // rekurencyjna 


return; // elementu brak 


>elem) usun(wsk—>lewy, x); 


danego 


NULL) { 
wsk; 


1/ 


zastąp prawym, nawet jeśli NULL 


delete usuwany; 
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else if (wsk-—>prawy == NULL) { 
wezelx* usuwany = wsk; 
wsk = wsk=>lewy; // lewy nie jest NULL; 
delete usuwany; 

) 

else 

( // ma obydwu potomków 


"t 


// zastąpimy "sprytnie usuwany węzeł jego prawym 


potomkiem 


// znajdź największy element mniejszy niż usuwany—> 
elem 

wezelx wsk2 = wsk—>lewy; 

while (wsk2-—>prawy != NULL) 
wsk2 = wsk2—>prawy; 


// elementy większe od wsk2->elem, ale mniejsze niż 

// wsk=>prawy—>elem będą "adoptowane" przez wsk2, 

// które będzie się znajdować gdzieś w lewym podrzewie 
wsk 

wsk2—>prawy = wsk—>prawy—>lewy; 


// powiększone lewe wsk podrzewo przenosimy jako lewe 
poddrzewo 

// prawego potomka wsk, żeby własności drzewa 

// binarnego zostały zachowane 


wsk->prawy-—>lewy = wsk->lewy; 

// zapamiętaj element usuwany 

wezel* usuwany = wsk; 

// prawy potomek usuwanego wskakuje na miejsce swego 
rodzica 


wsk = wsk=>prawy; 


delete usuwany; 
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Podziękowania. Serdecznie dziękuję Jowicie Hąci, Kasi Fokow i Michałowi Dębskiemu za 
pomoc w udoskonaleniu niniejszych materiałów. 
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3 Cwiczenia 


Zadanie 10.1. Napisz samodzielnie pełny program w języku C++, który implementuje i testuje 
(w funkcji main()) następujące operacje na liście jednokierunkowej przechowującej wartości 
typu double: 


a) Wyznaczenie sumy wartości wszystkich elementów. 

b) Wyznaczenie sumy wartości co drugiego elementu. 

c) Wyszukiwanie danego elementu. 

d) Wstawienie elementu na początek listy. 

e) Wstawienie elementu na koniec listy. 

f) Wstawienie elementu na 2-ta pozycję listy. 

g) Usuwanie elementu z początku listy. Usuwany element jest zwracany przez funkcję. 

h) Usuwanie elementu z końca listy. Usuwany element jest zwracany przez funkcję. 

i) Usuwanie elementu o zadanej wartości. Zwracana jest wartość logiczna w zależności od 
tego, czy element znajdował się na liście, czy nie. 

j) Usuwanie i-tego w kolejności elementu. Usuwany element jest zwracany przez funkcję. 


Zadanie 10.2. Rozwiąż powyższe zadanie, implementując listę jednokierunkową, która dodat- 
kowo przechowuje wskaźnik na ostatni element. 


Zadanie 10.3. Rozwiąż powyższe zadanie, implementując listę dwukierunkową. 


Zadanie 10.4. Dla danych dwóch list jednokierunkowych napisz funkcję, która je połączy, np. 
dla (1,2,5,4) oraz (3,2,5) sprawi, że I lista będzie postaci (1,2,5,4,3,2,5), a II zostanie skasowana. 


Zadanie 10.5. Napisz samodzielnie pełny program w języku C++, który implementuje i testuje 
(w funkcji main()) stos (LIFO) zawierający dane typu int. 


Zadanie 10.6. Napisz samodzielnie pełny program w języku C++, który implementuje i testuje 
(w funkcji main()) zwykłą kolejkę (FIFO) zawierającą dane typu char (napisy). 


x Zadanie 10.7. Zaimplementuj operacje enqueue() i dequeue() zwykłej kolejki (FIFO) typu 
int korzystając tylko z dwóch gotowych stosów. 


Zadanie 10.8. Napisz samodzielnie pełny program w języku C++, który implementuje i testuje 
(w funkcji main()) kolejkę priorytetową zawierającą dane typu int. 


Zadanie 10.9. Napisz funkcję, która wykorzysta kolejkę priorytetową do posortowania danej 
tablicy o elementach typu int. 


Zadanie 10.10. Napisz samodzielnie pełny program w języku C++, który implementuje i te- 
stuje (w funkcji main()) następujące operacje na drzewie binarnym przechowującym wartości 
typu double: 


a) Wyszukiwanie danego elementu. 

b) Zwrócenie elementu najmniejszego. 

c) Zwrócenie elementu największego. 

d) Wypisanie wszystkich elementów w kolejności od najmniejszego do największego. 

e) Wstawianie danego elementu. Jeśli wstawiany element znajduje się już w drzewie nie 
należy wstawiać jego duplikatu. 

f) Usuwanie danego elementu. Usuwany element jest zwracany przez funkcję. 
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4 Wskazówki do ćwiczeń 


Wskazówka do zadania [10.7| Dane są stosy A oraz B. 

Operacja enqueue(): odłóż element na stos A. 

Operacja dequeue(): zdejmij element ze stosu B. Jeśli stos jest pusty, przerzuć na niego 
wszystkie elementy ze stosu A korzystając z metod pop() i push(). 
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